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Abstract 

Past  research  has  indicated  that  implementation  of  a  pulsed  detonation  combustor 
(PDC)  into  a  high-bypass  turbofan  engine  yields  a  more  efficient  engine  at  design 
conditions.  It  is  proposed  that  performance  gains  can  be  made  utilizing  this  hybrid  engine 
off-design.  A  hybrid  high-bypass  turbofan  engine  with  a  PDC  model  was  evaluated  for  a 
range  of  Mach  numbers,  altitudes,  and  fill  fractions  in  the  Numerical  Propulsion  System 
Simulation  (NPSS).  Results  were  compared  to  a  conventional  baseline  high-bypass 
turbofan  engine  that  shares  the  same  architecture  with  the  hybrid.  The  NPSS  baseline 
engine  was  validated  using  the  Aircraft  Engine  Design  System  (AEDsys)  program  and 
the  net  thrust  and  specific  fuel  consumption  agreed  to  within  one  percent.  The  effect  of 
detonation  on  the  core  air  flow  is  calculated  using  a  closed  fonn  solution  for  the 
Chapman- Jouguet  Mach  number  with  a  total  energy  correction  applied.  Results  indicate 
that  fill  fraction  can  be  adjusted  to  reduce  the  TSFC  to  that  of  the  baseline  engine  and 
lower  at  some  thrust  levels.  With  careful  selection  of  design  parameters,  results  suggest  a 
pulsed  detonation  combustor  may  be  an  appropriate  candidate  for  inclusion  in  a  hybrid 
turbofan  engine. 
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OFF-DESIGN  ANAYSIS  OF  A  HIGH  BYPASS  TURBOFAN  USING  A  PULSED 


DETONATION  COMBUSTOR 

I.  Introduction 

Since  the  early  1940’s,  pulsed  detonation  engines  (PDE)  have  been  studied  as  a 
means  of  increasing  bum  efficiency  in  an  engine  as  a  result  of  its  supersonic  detonative 
mode  of  combustion  over  conventional  subsonic  deflagration.  Detonations  provide  a 
much  more  efficient  means  of  combusting  a  fuel-oxidizer  mixture  due  to  increased 
thennodynamic  efficiency  as  a  result  of  the  pressure-rise  associated  with  detonation. 
Additionally,  with  its  potential  for  a  cycle  time  of  more  than  ten  times  that  of  a  traditional 
pulsejet  engine  and  fewer  moving  parts  to  maintain,  PDEs  hold  the  promise  for 
applications  across  the  flight  envelope  spanning  subsonic,  supersonic  and  hypersonic 
flight. 

More  advanced  concepts  such  as  a  hybrid-PDE  have  been  studied  in  which  a 
pulsed  detonation  combustor  (PDC)  is  incorporated  into  a  gas  turbine  engine  as  the 
primary  combustion  system  with  the  intention  of  increasing  efficiency  by  utilizing  the 
strengths  of  both  engines.  In  this  type  of  system  the  exhaust  from  the  detonation  chamber 
drives  the  downstream  turbine  which  provides  power  to  the  compressor,  which,  in  turn, 
provides  the  air  flow  to  fill  and  purge  the  detonation  chamber.  Although  a  novel  idea,  the 
hybrid  system  is  not  without  its  challenges.  Low-vapor  pressure  hydrocarbon  fuels  must 
be  used  efficiently  as  key  PDE  cycle  parameters  such  as  ignition  time  and  deflagration  to 
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detonation  time  depend  on  the  properties  of  the  fuel.  Noise  is  also  a  substantial  issue  as 
detonations  are  significantly  louder  than  deflagration  combustion.  The  periodic,  high- 
pressure  pulses  must  be  assessed  on  turbine  performance  and  the  life  of  the  engine.  The 
first  flight  of  an  aircraft  powered  by  a  PDE  took  place  on  January  31,  2008  operating 
under  its  own  power  for  10  seconds  at  an  altitude  of  approximately  100  feet.  With  this 
demonstration  proving  that  a  PDE  can  be  integrated  into  an  aircraft  frame  without 
experiencing  structural  problems,  PDEs  are  increasingly  recognized  as  a  realizable 
technology  for  future  aerospace  propulsion. 

Purpose 

A  substantial  amount  of  work  on  PDEs  and  hybrid  PDEs  has  been  accomplished, 
with  significant  developments  being  made  in  the  last  fifteen  years.  The  theoretical 
analysis  of  Petters  and  Felding  (Petters  and  Felding,  2002:6)  indicated  that  a  PDE-hybrid 
with  the  same  inlet  airflow  as  a  baseline  turbofan  engine  produced  a  2%  higher  thrust  and 
an  1 1%  reduction  in  thrust  specific  fuel  consumption  (TSFC).  Similar  studies  by  Andrus 
(Andrus,  2007:81)  showed  that  an  optimal  hybrid  engine  operating  at  design  conditions 
could  yield  an  8%  decrease  in  TSFC  while  maintaining  thrust. 

The  experimental  work  of  Schauer  et  al.  (Schauer  et  ah,  2003),  Deng  et  al.  (Deng 
et  al.,  2008),  and  Rasheed  et  al.  (Rasheed  et  al.,  2005)  all  investigated  a  detonation  driven 
turbine  at  design  conditions;  however,  the  only  experimental  work  done  on  PDEs  at  off- 
design  conditions  to  date  is  that  of  Glaser  et  al.  (Glaser  et  al.,  2004).  Glaser’s  work 
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suggests  performance  gains  can  be  made  by  varying  key  parameters  such  as  equivalence 
ratio  and  fill  fraction.  The  objective  of  this  thesis  is  to  build  on  Andrus’s  (Andrus,  2007) 
work  by  perfonning  a  simulated  off-design  analysis  of  a  hybrid-PDE  design  to  evaluate 
the  effects  on  thrust  and  TSFC. 

Procedure 

The  procedure  for  performing  the  off-design  analysis  closely  mirrors  the  steps 
used  by  Andrus  in  performing  his  comparative  analysis  of  a  high  bypass  turbofan  using  a 
pulsed  detonation  combustor  with  a  conventional  baseline  turbofan.  A  baseline  high- 
bypass  turbofan  is  modeled  in  both  the  Numerical  Propulsion  System  Simulation  (NPSS) 
and  in  the  Aircraft  Engine  Design  system  (AEDsys)  programs.  The  comparison  is  made 
to  ensure  identical  engine  configuration  between  the  two  programs.  The  NPSS  program  is 
the  primary  software  program  used  in  this  thesis  to  evaluate  the  off-design  perfonnance 
of  a  hybrid-PDE.  It  was  developed  via  a  cooperative  effort  between  industry  and  NASA 
to  predict  and  analyze  the  aero  thermodynamic  behavior  of  commercial  jet  aircraft, 
military  applications,  and  space  transportation  with  the  goal  of  reducing  development 
time  and  cost  of  a  new  engine  by  half.  AEDsys  was  developed  by  Mattingly  (Mattingly  et 
ah,  2006)  for  educational  use  in  the  field  of  gas  turbine  engine  design  and  allows  the  user 
to  perfonn  design  point  and  parametric  cycle  analysis  for  various  engines. 

After  validating  the  accuracy  of  the  baseline  engine  in  both  AEDsys  and  NPSS,  a 
hybrid-PDE  model  with  the  identical  configuration  as  the  baseline  engine,  with  the 
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exception  of  a  pulsed  detonation  combustor  replacing  the  conventional  combustor,  is 
developed  in  NPSS  and  run  at  off-design  conditions.  An  analysis  of  the  hybrid  engine 
performance  is  evaluated  for  a  range  of  Mach  numbers,  altitudes,  and  fill  fractions.  The 
effects  on  thrust  and  TFSC  are  compared  to  that  of  the  NPSS  baseline  engine  running  at 
the  same  off-design  conditions. 


Significance  of  Research 

Much  research  has  been  done  on  turbofan  engines  with  a  pulsed  detonation 
combustor  at  design  conditions,  but  there  is  a  very  limited  amount  of  literature  on  the 
performance  of  these  engines  off-design.  Because  the  majority  of  an  engine’s  operation 
are  at  off-design  conditions,  significant  cost  savings  could  be  realized  if  hybrid  turbofan 
engines  are  more  efficient  than  conventional  engines.  In  addition,  a  hybrid  engine  may  be 
cheaper  to  build  and  less  expensive  to  maintain  than  a  conventional  engine,  offering 
additional  long  tenn  savings. 

Organization 

This  thesis  compares  the  perfonnance  of  a  hybrid-PDE  to  that  of  a  conventional 
turbofan  using  the  NPSS  program.  Chapter  two  contains  a  thorough  discussion  of  pulsed 
detonation  thermodynamics  as  well  as  prior  work  on  PDEs.  Chapter  three  describes  the 
baseline  turbofan  model  used  in  AEDsys  and  NPSS,  as  well  as  the  hybrid-PDE  engine 
and  its  combustor  section.  Chapter  four  is  a  comprehensive  analysis  of  the  hybrid-PDE 
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performance  as  compared  to  that  of  the  conventional  baseline  engine.  Chapter  five 
contains  the  conclusions  of  this  research  as  well  as  recommendations  for  future  work. 
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II.  Literature  Review 


Introduction 

This  section  presents  a  thorough  discussion  of  the  underlying  thermodynamics  of 
the  hybrid-PDE  engine  as  the  results  of  this  thesis  are  essentially  governed  by  the  basic 
models  of  detonation  to  include  the  Chapman- Jouguet  and  Zeldovich-von  Nuemann- 
Doring  (ZND)  theories. 

The  Chapman- Jouguet  theory  allows  for  the  calculation  of  the  detonation  velocity 
of  a  detonation  wave  with  known  pressures  and  densities  of  the  unburned  gases  for  a 
given  q.  The  steady  state  solution  of  the  detonation  wave  requires  knowledge  of  the 
equilibrium  thennodynamic  calculations.  Experimental  results  have  shown  to  agree  well 
with  the  detonation  velocities  resulting  from  this  theory. 

Zeldovich,  von  Nuemann,  and  Doring  (Kuo,  1986)  present  a  model  for  detonation 
wave  structure  in  which  parameters  such  as  detonation  limits,  initiation  energy,  tube 
diameter,  etc.  are  known.  Unlike  the  Chapman- Jouguet  theory,  experimental 
measurements  do  not  agree  with  the  model  calculations,  mainly  because  the  ZND 
structure  is  unstable  and  only  observed  experimentally  under  transient  conditions. 
Experimental  observations  show  that  the  self-sustained  detonations  have  a  three- 
dimensional  cell  structure;  however,  there  are  currently  no  acceptable  theories  that  define 
this  cell  structure. 
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Combustion  Waves 


In  order  to  understand  how  a  pulsed  detonation  combustor  can  be  more  efficient 
than  a  conventional  combustor,  it  is  first  necessary  to  understand  the  differences  between 
the  detonations  and  deflagrations  of  the  two  burners  that  produce  engine  thrust.  A 
detonation  is  a  supersonic  shock  wave  that  propagates  through  a  fluid  due  to  an  energy 
release  in  a  reaction  zone.  A  deflagration  is  a  wave  that  propagates  at  a  subsonic  rate  by 
heat  transfer.  Detonations  generate  higher  pressures  and  have  increased  wave  speeds, 
thus  producing  greater  thrust  than  deflagrations.  Figure  2.1  shows  a  schematic  of  a 
stationary  one-dimensional  combustion  wave  in  which  subscript  one  and  two  denote 
conditions  of  the  unburned  gases  ahead  of  the  wave  and  burned  gases  behind  the  wave, 
respectively.  Deflagration  and  detonation  wave  properties  are  compared  in  Table  2.1. 


SlQtl^Raf y  Ca-mbuilion  Wave 


{Unburnfld? 

u. 


k" 


i  H  u  r  i' id  I 


Pi.  t2.  pi 


Figure  2.1  Stationary  one-dimensional  combustion  wave  (Kuo,  1986:233) 
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Table  2.1.  Qualitative  differences  between  detonation  and  deflagration  in  gases 
(Kuo,  1986:234) 


Detonation 

Deflagration 

iq/at 

5-10 

0.0001-0.03 

u2/ut 

0.4-0. 7 

4-6 

P2/P1 

13-55 

-0.98 

T/Ti 

8-21 

4-16 

P2/P1 

1. 7-2.6 

0.06-0.25 

A  combustion  wave  is  fonned  in  a  tube  when  a  combustible  gas  mixture  is  ignited 
at  the  closed  end  of  a  tube.  The  properties  in  Table  2.1  show  that  the  burned  gases  are 
higher  in  temperature  and  density  than  the  unburned  gases.  This  increase  in  density 
initiates  a  compression  wave  that  travels  towards  the  deflagration  wave  front,  causing  the 
wave  to  accelerate.  Density  increases  as  the  deflagration  wave  continues,  causing  more 
and  more  compression  waves  to  form.  The  waves  accelerate  as  pressure  and  density 
increase,  thus  causing  them  to  amalgamate  at  the  deflagration  wave  front.  If  the  tube  is 
sufficiently  long,  a  shock  wave  will  form  that  is  strong  enough  to  ignite  the  mixture 
ahead  of  the  wave  front.  A  detonation  is  obtained  as  the  continuous  compression  waves 
in  the  reaction  zone  keep  the  shock  from  decaying.  The  detonation  is  inherently  self- 
sustaining  in  that  the  detonation  front  initiates  a  chemical  reaction  by  compression  by 
diffusing  heat. 


Chapman-Jouguet  Theory 

To  solve  for  the  Hugoniot  curve  on  which  the  Chapman-Jouguet  points  are  found, 
we  must  first  start  with  the  conservation  equations  (Glassman,  1996:226): 
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Pllll  =  p^u2 


(2.1) 


P\  +  PlUl  =  Pi  +  PlUl 


(2.2) 


CpT\  +  ^  ul  +  Q  Cpl 2  ^  W2 


(2.3) 


Pi  =  PiRTi 
Pi  =  PiRTi 


(2.4) 


where  Eqs.  2.1,  2.2,  and  2.3  are  the  mass,  momentum,  and  energy  respectively.  Equation 
2.4  is  simply  the  equation  of  state.  The  four  equations  can  be  reduced  to  one  equation 
with  two  unknowns,  p2  and  po,  by  combining  Eqs.  2.1  and  2.2  to  yield  (Kuo,  1986:236): 


AX  = 


Pz  ~P\ 

1  1 


A  A 


(2.5) 


or  in  terms  of  Mach  number: 


rMt  = 


Pi_ 
2  _  Pi 


-1 


1- 


Pl_ 

Pi 


(2.6) 


where  Eqs.  2.5  and  2.6  are  known  as  the  Raleigh-line  relation. 

The  Hugoniot  relation  can  be  found  by  combining  Eqs.  2.1  -  2.4  to  yield 
(Glassman,  1996:228): 


Y 


7-1 


' P±_  Pi 

V  Pi  Ay 


-\{Pl-Pl) 


X__j_ 

V  A  Pi ) 


(2.7) 
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where  q  is  the  heat  flux.  The  Hugoniot  curve  as  shown  in  Fig.  2.2  is  a  plot  of  the  specific 
volume  (1/p)  and  pressure  of  the  burned  gases  for  given  values  of  heat  flux,  specific 
volume  and  pressure  of  the  unbumed  gases. 


Figure  2.2  Hugoniot  curve  with  Rayleigh  lines  (Kuo,  1986:238) 

The  plot  is  broken  up  into  five  regions  constructed  by  drawing  tangents  and 

vertical  and  horizontal  lines  from  the  origin  (l/pi,pi)  to  the  curve.  The  two  Chapman- 

Jouguet  points  are  at  the  tangents  to  the  curve  and  are  referred  to  as  the  upper  and  lower 

C-J  points  at  the  upper  and  lower  Raleigh  lines,  respectively.  Of  the  five  regions,  only 

regions  I,  II,  and  III  are  physically  possible.  Region  V  does  not  bound  a  valid  solution  as 

P2  and  l/p2  are  greater  than  pi  and  1/pi  and  thus  would  require  a  compression  wave  to 

move  in  the  negative  direction.  Region  IV  is  also  ruled  out  as  the  heat  addition  stipulates 

supersonic  flow;  however,  it  is  not  possible  to  have  heat  addition  and  advance  past  the 

sonic  condition  in  a  constant  area  duct.  Regions  I  and  II  are  the  detonation  regions  of  the 
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curve  and  represent  strong  and  weak  detonations,  respectively.  These  regions  are 
eliminated  due  to  the  structure  of  the  detonation  wave  discussed  in  the  next  section. 
Region  III  is  the  weak  deflagration  region  and  is  often  observed  in  most  experimental 
conditions;  however,  since  deflagration  is  not  of  interest  in  this  thesis,  region  III  is 
irrelevant  to  this  work  as  well.  The  upper  C-J  point  is  of  importance  to  this  research  in 
that  the  wave  speed  at  this  point  corresponds  to  a  minimum  detonation  wave  speed  and 
implies  that  the  Mach  number  of  the  burned  gases  must  be  equal  to  one.  The  method  used 
in  this  work  to  calculate  the  velocity  of  the  wave  at  this  point  will  be  discussed  later  in 
this  chapter. 

Zeldovich-von  Nuemann-Doring  Theory 

The  Zeldovich-von  Nuemann-Doring  (ZND)  model  is  a  one-dimensional  model 
of  the  structure  of  a  detonation  wave.  The  model  assumes  a  one-dimensional,  steady  flow 
with  limited  reactions  in  the  shock  wave  region.  As  shown  in  Fig.  2.3  the  detonation 
wave  consists  of  a  thin  shock  wave  region  followed  by  a  thick  deflagration  region 
consisting  of  an  induction  and  reaction  zone.  The  reactants  are  initially  heated  by  the 
shock  wave  to  a  temperature  which  ensures  a  high  enough  reaction  rate  in  which  the 
deflagration  can  propagate  at  the  same  speed  as  the  shock  wave.  The  thin  shock  layer 
results  in  a  sharp  spike  in  temperature,  pressure,  and  density.  The  peak  pressure  reached 
in  this  region  is  referred  to  as  the  von  Neumann  spike.  This  is  followed  by  relatively 
steady  profiles  through  the  induction  zone  due  to  a  slowly  increasing  rate  of  reaction 
immediately  behind  the  shock  front. 
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Figure  2.3  Thermodynamic  property  variations  across  a  ZND  detonation  wave 
(Kuo,  1986:261) 

The  properties  change  drastically  again  as  reaction  rate  increases  and  then  reach  their 
equilibrium  values  once  the  reaction  has  completed.  The  ZND  detonation  structure  may 
also  be  shown  on  Fig.  2.4  beginning  at  the  Hugoniot  origin  and  moving  up  along  the  left 
of  the  curve  until  it  reaches  the  von  Neumann  spike.  At  this  point  the  pressure  decreases 
and  the  path  merges  with  the  Hugoniot  curve  to  the  upper  C-J  point. 

Although  these  models  assume  a  detonation  to  be  one-dimensional,  it  should  be 
acknowledged  that  detonation  waves  moving  in  tubes  are  actually  three-dimensional  and 
nonsteady  in  nature  in  which  the  flow  proceeds  in  a  cyclic  manner  with  shock  velocity 
fluctuations  about  the  equilibrium  C-J  value. 
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Figure  2.4  Detonation  structure  (Williams,  1985:193) 


Thermodynamic  Cycle  Analysis 

The  pulsed  detonation  engine  thermodynamic  cycle  is  described  by  Heiser  and 
Pratt  (Heiser  and  Pratt,  2002:2)  as  being  identical  to  that  of  an  ideal  Humphrey  cycle  used 
in  turbojets  and  ramjets  with  the  exception  of  heat  addition  during  the  combustion 
process.  The  ideal  Humphrey  cycle  is  sometimes  used  to  estimate  the  thermal  efficiency 
of  the  PDE  cycle  as  it  replaces  the  Brayton  cycle’s  constant-pressure  heat  addition 
process  with  a  constant-volume  heat  addition  process.  The  T-s  diagram  of  these  three 
processes  is  shown  in  Fig.  2.5. 
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Figure  2.5  Ideal  PDE,  Humphrey,  and  Brayton  cycle  temperature-entropy  diagrams 
(Heiser,  2002:4) 

The  thermal  efficiency  of  the  PDE  cycle  is  slightly  greater  than  that  of  the 
Humphrey  cycle  and  much  greater  than  that  of  the  Brayton  cycle.  The  thermal  efficiency 
of  the  PDE  cycle  as  proposed  by  Heiser  is  identical  to  that  of  the  Fickett-Jacobs  cycle  as 
described  by  Wintenberger  (Wintenberger  and  Sheperd,  2004: 12)  in  which  an  upper  limit 
is  computed  to  be  the  amount  of  mechanical  work  in  a  cycle  produced  by  an  unsteady 
detonation  process. 

In  Fig.  2.5,  from  state  0  to  3,  an  isentropic,  adiabatic  compression  takes  place  in 
all  three  cycles,  raising  the  temperature  to  T3.  It  is  the  process  from  state  3  to  4  in  which 
the  PDE,  Humphrey,  and  Brayton  cycles  differ.  In  the  Brayton  cycle,  a  constant  pressure 
heat  addition  takes  place  and  increases  the  temperature  T3  of  the  combustor  inlet  to  T4  at 
the  combustor  outlet.  From  state  3  to  4  in  the  PDE  cycle,  the  ZND  detonation  wave  is 


14 


seen  in  which  the  process  is  constrained  by  the  Chapman- Jouguet  condition,  requiring  the 
Mach  number  at  state  4  to  be  sonic.  The  path  from  state  3  to  4  differs  slightly  from  that  of 
the  Brayton  and  Humphrey  cycles  in  that  from  state  3  to  3  a  the  heat  addition  process 
generates  entropy  via  the  adiabatic  normal  shock  wave,  and  from  state  3a  to  4  entropy  is 
generated  via  a  constant  area  heat  addition  process.  The  process  from  state  4  to  10  and 
state  10  to  0  of  the  three  cycles  are  identical  in  that  an  isentropic  expansion  process  takes 
place  followed  by  a  heat  rejection  to  close  the  cycle. 

The  derivation  for  the  solutions  for  the  Chapman- Jouguet  Mach  number,  the 
entropy  difference  from  states  3  to  4,  and  cycle  thennal  efficiency  are  shown  in  Appendix 
A.  They  are  given  here  by  Eqs.  2.8,  2.9,  and  2.10  respectively: 
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The  Chapman- Jouguet  Mach  number  is  calculated  using  the  non-dimensional  heat 
addition  q  and  y/  ,  the  ratio  of  static  temperature  at  state  3  to  the  free-stream  static 
temperature  at  state  0.  Mcj  is  then  used  to  calculate  the  entropy  rise  from  state  3  to  4  and 
the  thermal  efficiency  of  the  cycle.  According  to  Heiser  and  Pratt  (Heiser  and  Pratt, 

2002),  all  the  fluid  properties  at  the  detonation  tube  exit  can  be  solved  for  using  Eqs.  2.9 
and  2. 13  to  solve  for  the  entropy  and  pressure  at  state  4. 

?±=  (2.13) 

Po  r + 1  Po 

The  PDE  thennodynamics  as  described  by  Heiser  and  Pratt  is  used  in  this  thesis 
with  the  addition  of  a  correction  factor  by  Dyer  and  Kaemming  (Dyer  and  Kaemming, 
2002:5).  They  note  that  using  the  pressure  and  entropy  at  the  detonation  tube  exit  to  solve 
for  all  the  properties  of  the  fluid  at  this  point  is  inaccurate  because  it  ignores  the  eventual 
pressure  loss  that  the  gas  will  go  through  due  to  expansion  waves.  They  propose  using 
entropy  and  the  change  in  enthalpy  liberated  by  the  combustion  process  to  solve  for  the 
properties  at  the  detonation  tube  exit.  Available  energy  is  calculated  using  the  known  CJ 
entropy  originally  calculated  by  Heiser  and  Pratt  (Heiser  and  Pratt,  2002:3),  Eq.  2.9,  with 
the  known  system  enthalpy  of  (h0+qadd),  with  qadd  being  the  heat  flux  into  the  system,  to 
ensure  that  energy  is  conserved. 

Pulsed  Detonation  Engine  Cycle 

A  basic  understanding  of  the  PDE  cycle  is  necessary  to  understand  how  the 
combustor  section  of  the  hybrid  PDE  performs  in  this  research.  The  cycle  consists  of  four 
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distinct  processes:  fill  phase,  detonation  initiation,  detonation  wave  propagation,  and 
purge  phase.  The  combustor  is  then  recharged  with  another  fuel/air  mixture  and  the  cycle 
repeats. 


Figure  2.6  illustrates  the  various  stages  in  a  pulsed  detonation  engine  cycle.  At 


©Purging  stage 

(chamber  is  at  ambient  condition) 


PDE  Cycle 


Figure  2.6  PDE  cycle  stages 

station  one  in  the  diagram,  the  detonation  chamber  is  at  ambient  conditions.  The  fill 
phase  is  shown  at  station  2  as  a  valve  that  seals  one  end  of  the  detonation  chamber  opens, 
permitting  the  fuel/air  mixture  into  the  chamber.  The  volume  of  the  fuel/air  mixture  at 
ambient  conditions  to  the  tube  volume  is  the  fill  fraction  iff).  This  is  one  of  the  variables 
analyzed  in  this  research  to  determine  engine  performance  over  a  range  of  values. 

After  the  fuel/air  mixture  enters  the  chamber,  the  valve  closes  and  a  detonation 
wave  is  initiated  near  the  closed  end  of  the  chamber  as  shown  at  station  3  of  the  figure. 

At  the  onset  of  this  stage,  a  spark  plug  deposits  a  spark  that  causes  a  deflagration  wave  to 
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form  at  the  end  of  the  tube.  The  deflagration  wave  propagates  through  the  tube  and 
transitions  to  a  detonation  wave  before  reaching  the  open  end  of  the  tube.  This  transition 
is  known  as  the  deflagration  to  detonation  transition  (DDT),  and  the  time  elapsed 
between  the  fonnation  of  the  deflagration  wave  and  detonation  wave  is  known  as  the 
DDT  transition  time.  The  detonation  wave  then  propagates  to  the  tube  exit  at  the 
Chapman- Jouguet  condition.  As  shown  by  station  4  in  Fig.  2.6,  the  region  ahead  of  the 
detonation  wave  contains  the  unbumed  gas  at  state  one.  The  burned  gas  at  state  2  just 
behind  the  wave  is  at  a  significantly  higher  pressure  and  temperature  than  state  one; 
however,  the  burned  gas  at  state  3  near  the  near  the  closed  end  of  the  chamber  will  have  a 
lower  temperature  and  pressure  than  the  gas  at  state  2  with  an  intennediate  condition 
existing  between  states  2  and  3. 

Upon  reaching  the  end  of  the  tube,  the  detonation  wave  exits,  producing  the  thrust 
of  the  engine.  The  purge  phase  begins  as  a  pressure  differential  in  the  tube  creates 
rarefaction  waves  which  propagate  into  the  tube  and  expel  the  burned  gases.  Pressure  and 
temperature  in  the  chamber  eventually  decay  to  ambient  levels  and  the  exhaust  velocity 
goes  to  zero.  The  detonation  tube  can  then  be  filled  with  a  new  fuel/air  mixture  to  begin 
the  cycle  once  again. 

Prior  Work  on  Hybrid-Pulse  Detonation  Engines 

A  significant  amount  of  research,  both  experimental  and  analytical,  has  been  done 
on  integrating  a  pulsed  detonation  combustor  into  a  turbine  system  with  the  hopes  of 
increasing  thrust  and  decreasing  fuel  consumption  of  an  aircraft  engine.  Petters  and 
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Felder  (Petters  and  Felder,  2002:4)  and  Andrus  (Andrus,  2007:81)  demonstrated  that,  in 
theory,  a  pulsed  detonation  hybrid  engine  can  reduce  TSFC  by  8  to  10%  at  design 
conditions.  GE  Global  Research  (Rasheed  et  ah,  2005)  built  and  evaluated  an  eight  tube 
PDC  integrated  with  a  single-stage  axial  turbine.  Results  indicated  the  ability  to  produce 
detonations  at  10  and  20  Hz  conditions  showing  promise  for  operability  over  a  wide 
range  of  conditions.  Noise  signatures  and  internal  structural  damage  due  to  the  cyclic 
pulsations  of  the  detonations  are  a  cause  for  concern  in  implementing  a  PDC  into  a 
turbine  system.  Caldwell  and  Gutmark  (Caldwell  and  Gutmark,  2008:1)  perfonned 
experimental  studies  to  ascertain  the  flow  field  and  suggest  that  shock  reflection  and 
blowdown  jet  interaction  length  and  time  scales  could  minimize  noise  and  structural 
damage.  During  Schauer  et  al.  (Schauer  et  ah,  2003: 1)  testing  of  a  PDE  into  a  radial 
turbine,  the  turbine  withstood  all  detonations  into  the  inlet,  as  well  as  significantly 
weakened  the  strength  of  the  detonation  shocks  in  the  exhaust  nozzle.  This  experimental 
work  (Rasheed  et  al.,  2005)  (Caldwell  and  Gutmark,  2008)  (Schauer  et  al.,  2003)  among 
others  prove  that  after  overcoming  a  few  hurdles,  these  engines  can  become  a  reality. 

Though  much  work  has  been  performed  on  hybrid-PDEs  at  design  conditions, 
hardly  any  off-design  analysis,  either  experimentally  or  analytically,  has  been 
accomplished.  Off-design  analysis  detennines  the  performance  of  an  engine  at  a  given  set 
of  conditions  for  a  fixed  geometry  determined  from  a  design  operating  point.  Glaser  et  al. 
(Glaser  et  al.,  2004:1)  experimentally  investigated  the  off-design  perfonnance  of  a  pulsed 
detonation  engine  by  varying  the  equivalence  ratio  and  fill  fraction  parameters.  Their 
PDE  system  utilized  a  single  stainless  steel  PDE  tube  1”  in  diameter  and  24”  in  length.  A 
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fuel/oxidizer  mixture  of  ethylene  and  oxygen  was  injected  into  the  tube  via  a  controlled 
solenoid  valve.  A  relationship  between  the  wave  speed  and  the  equivalence  ratio  was 
found  and  is  shown  in  Fig.  2.7. 


?roo 


Figure  2.7  Relationship  between  equivalence  ratio  and  detonation  wave  speed  (Glaser  et  al.,  2004:6) 


As  can  be  seen  from  the  plot,  the  wave  speed  increases  with  equivalence  ratio 
before  leveling  out  at  a  maximum  equivalence  ratio  of  1.7  and  wave  speed  of  2583  m/s. 
The  effect  of  fill  fraction  on  wave  speed  was  also  detennined.  These  results  can  be  seen 
in  Fig.  2.8. 
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Figure  2.8  Relationship  between  fill  fraction  and  detonation  wave  speed  (Glaser  et 
al.,  2004:7) 


Two  different  air/fuel  ratios  were  investigated,  with  the  two  mixtures  diverging  at 
a  fill  fraction  of  approximately  0.6  before  leveling  off  at  a  fill  fraction  of  about  1.0.  These 
results  indicate  that  perfonnance  gains  may  be  made  at  an  equivalence  ratio  greater  than 
one  and  a  completely  filled  detonation  tube.  Although  these  tests  were  not  performed  on 
a  hybrid-PDE,  they  indicate  favorable  results  for  pulsed  detonation  off-design  studies, 
thus  furthering  the  need  for  hybrid-PDE  off-design  research. 

Chapter  three  will  describe  the  baseline  and  hybrid-PDE  models  in  AEDsys  and 
NPSS  and  detail  how  the  pulsed  detonation  thennodynamics  are  incorporated  into  the 
hybrid  combustor  section  to  solve  for  the  hybrid  engine  perfonnance. 
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III.  Baseline  and  Hybrid  Models 


Introduction 

To  ensure  the  accuracy  of  the  baseline  model  in  NPSS,  it  was  compared  to  the 
same  baseline  model  in  AEDsys  (Mattingly  et  al.,  2002)  at  the  design  point.  The  baseline 
turbofan  models  were  run  in  off-design  mode  in  both  the  AEDsys  and  NPSS  programs 
and  the  results  compared.  The  hybrid  turbofan  model  was  developed  and  run  off-design 
and  its  perfonnance  compared  to  that  of  the  baseline  model  off-design.  This  chapter 
describes  the  baseline  and  hybrid  models,  including  the  changes  made  to  the  hybrid 
model  to  perform  an  off-design  analysis. 

Baseline  High  Bypass  Turbofan  Engine  in  AEDsys  and  NPSS 

The  modeled  engine  is  based  on  the  parameters  of  the  TF-39-GE-1C  engine  used 
on  the  C-5  Galaxy,  as  this  engine  has  known  operating  parameters  and  is  relevant  to  the 
Air  Force.  The  component  efficiencies  were  unknown,  however,  and  were  selected  to 
correspond  to  a  technology  level  projected  ten  to  twenty  years  in  the  future.  The 
efficiencies  can  be  found  in  Mattingly’s  Table  4.4  (Mattingly  et  al.,  2002: 107)  under 
level  4  technology.  Table  3.1  shows  the  parameters  for  the  TF-39-GE-1C  engine  as 
compared  with  the  notional  baseline  of  the  engine  modeled. 
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Table  3.1  Parameters  for  baseline  engine 


Engine 

TF39-GE-1C 

Baseline  Model 

Aircraft 

C-5 

N/A 

Mach  number 

0.0 

0.0 

Altitude  (ft) 

0.0 

0.0 

Fn(lbf) 

40805 

41500 

Weight  (lb) 

7186 

N/A 

Length  (in) 

100 

N/A 

Max  Diam  (in) 

203 

N/A 

BPR 

8 

8 

FPR 

1.56 

1.56 

OPR 

26 

26 

TSFC (l/h) 

0.315 

0.325 

Mass  flow  ( Ibm/s) 

1549 

1500 

Max  TIT  (R) 

2810 

2900 

Some  of  the  baseline  parameters  are  slightly  different  than  those  for  the  TF39-GE-1C 
engine  and  were  chosen  for  simplicity. 

The  baseline  engine  coded  in  NPSS  utilizes  the  architecture  of  a  high  bypass  split 
stream  turbofan  as  described  by  Mattingly  (Mattingly  et  ah,  2002:569-587).  The  model 
reference  stations  and  NPSS  configuration  are  shown  in  Fig.  3.1. 


NozS 
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Figure  3.1  Baseline  NPSS  high  bypass  turbofan  engine  configuration  with  reference  stations 
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In  Fig.  3.1,  the  core  and  bypass  flows  split  at  the  fan  entry  and  two  mixers  are 


employed,  one  at  the  burner  exit  to  combine  the  fluid  exiting  the  burner  with  bleed  flow 


(MIX40)  and  another  at  the  high  pressure  turbine  exit  to  combine  the  fluid  exiting  the 


turbine  with  bleed  flow  (MIX44).  The  bleeds  are  5.0%,  which  is  the  default  value  in 


AEDsys.  The  model  file  that  defines  the  baseline  NPSS  engine  is  found  in  Appendix  B. 

The  engine  was  also  modeled  in  AEDsys  in  order  to  compare  results  to  the  NPSS 
model  and  ensure  its  accuracy.  Table  3.2  shows  the  input  variables  for  the  AEDsys 
baseline  turbofan  engine. 

Table  3.2  AEDsys  baseline  turbofan  engine  input  parameters 


Input  Variables 

Value 

Mach  Number 

0.0 

Altitude 

0.0  ft 

Temperature 

518.67  R 

Pressure 

14.696  psia 

Mass  Flow 

1500  Ibm/s 

Tt4 

2900  R 

Compr  Press  Ratio 

26 

LPC  PR 

1.56 

Fan  PR 

1.56 

Bypass  Ratio 

8 

Fuel  Heating  Value 

18400  BTU/lbm 

Cp  Compressor 

0.2415  BTU/lbm  R 

gamma  c 

1.3986 

Cp  Turbine 

1.2957 

Bleed  Air 

1% 

Coolant  Air 

5% 

Coolant  Air  2 

5% 

Diffuser  Max  Press  Ratio 

0.995 

Burner  Press  Ratio 

0.96 

Nozzle  Press  Ratio 

0.985 

Fan  Press  Ratio 

0.98 

Fan  Polytropic  Eff 

0.89 

LPC  Polytropic  Eff 

0.89 

HPC  Polytropic  Eff 

0.90 

HPT  Polytropic  Eff 

0.89 

LPT  Polytropic  Eff 

0.89 

Burner  Eff 

0.995 

Mech  Shaft  LP  Spool  Eff 

0.99 

Mech  Shaft  HP  Spool  Eff 

0.99 

Max  Tt4 

3200  R 

Max  Comp  Press  Ratio 

32 

Max  Press  at  station  3 

650  psia 

Max  Temp  at  station  3 

1860  R 

Max  %  Ref  RPM  -  LP  Spool 

110 

Max  %  Ref  RPM  -  HP  Spool 

110 
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All  parameters  in  Table  3.2  are  also  input  into  NPSS  with  the  exception  of  the 


polytropic  efficiencies.  NPSS  requires  adiabatic  efficiency  inputs  for  turbines  and 
compressors  as  compared  to  the  AEDsys  requirement  for  polytropic  efficiency  inputs. 
The  two  efficiencies  are  related  by  Eqs.  3.1  and  3.2  (Oates,  1997:  214  &  222): 

n  (n~1)/r<  -1 


>h  = 


n. 


(.Yc^VY^c 


-1 


(3.1) 


n, 


1  ~n. 


e,(Y,-Wr, 


I  -71 


0,-iVr, 


(3.2) 


These  relationships  are  used  to  calculate  the  adiabatic  efficiencies  used  in  the  baseline 
NPSS  model  listed  in  Table  3.3. 

Table  3.3  Baseline  model  NPSS  adiabatic  efficiency  inputs 


NPSS  efficiency  inputs 

Value 

Fan  adiabatic  eff 

0.8827 

LPC  adiabatic  eff 

0.8827 

HPC  adiabatic  eff 

0.8573 

HPT  adiabatic  eff 

0.9057 

LPT  adiabatic  eff 

0.9084 

The  other  difference  between  the  AEDsys  and  NPSS  programs  is  the 
thermodynamic  model.  The  AEDsys  thennodynamic  package  is  a  subroutine  termed 
FAIR.  FAIR  is  an  8th  order  polynomial  fit  to  JANAF  specific  heat  data  for  pure  air,  and 
CEA  data  for  vitiated  air  (Mattingly  et  ah,  2002:  89-91).  NPSS  however,  gives  the  user 
control  over  the  thennodynamic  quantities  by  offering  a  choice  of  six  different 
thennodynamic  packages.  Table  3.4  lists  the  available  models  and  provides  a  description 
of  each.  The  NPSS  data  in  this  thesis  was  generated  using  the  GasTbl  thermodynamics 
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package  as  it  is  the  simplest  to  implement  in  NPSS  and  a  close  match  to  the  AEDsys 
thennodynamic  package.  However,  the  GasTbl  package  is  limited  to  results  whose 
equivalence  ratio  is  less  than  one;  therefore,  solutions  with  equivalence  ratios  greater  than 
one  could  not  be  investigated. 


Table  3.4  NPSS  Thermodynamics  Packages 


Package 

Description 

CEA 

Implementation  of  NASA  chemical 
equilibrium  code 

Janaf 

Implementation  of  NIST  gas  properties 
prepared  by  Honeywell 

GasTbl 

Package  developed  by  Pratt  &  Whitney 
based  on  Therm,  but  adding  humidity 
calculations  and  some  chemical 

AIIFuel 

Package  developed  by  General  Electric  that 
contains  gas  properties  and  fuel  properties 

INGTherm 

Package  used  to  calculate  the  properties  of 
ideal  gases,  contains  only  properties  for 
Helium  and  Zenon  at  present 

FPT 

Package  used  to  define  NPSS  table  and/or 
functions  that  describe  the  thermodynamic 
properties  of  the  fluid 

The  differences  in  the  thennodynamics  packages  account  for  a  1 .0% 
difference  in  net  thrust  and  0.8%  difference  in  TSFC  at  SLS  of  the  baseline  engine  in 
NPSS  and  AEDsys  at  the  design  point.  Andrus  (Andrus,  2007:23)  provides  a  complete 
explanation  on  how  the  thermodynamics  packages  differ  and  how  they  contribute  to  these 
differences. 
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Hybrid  Turbofan  Engine  and  Pulsed  Detonation  Combustor  in  NPSS 

The  hybrid-PDE  model  contains  the  identical  architecture  of  the  baseline  turbofan 
model  shown  in  Fig.  3.1,  with  the  exception  of  the  burner  section  which  is  replaced  by  a 
pulsed  detonation  combustor  (PDC).  The  NPSS  Model  code  for  the  hybrid  turbofan 
engine  can  be  found  in  Appendix  C. 

The  hybrid  engine  and  PDC  inputs  are  listed  in  Table  3.5. 

Table  3.5  Hybrid  engine  on-design  configuration 


Input  Variables 

Value 

Mass  Flow 

1500  lb  m/s 

Tt4 

2900  R 

Compr  Press  Ratio 

26 

LPC  PR 

1.56 

Fan  PR 

1.56 

Bypass  Ratio 

8 

Tube  Inner  Diameter 

2.1  in 

Tube  Length 

36  in 

Number  of  Tubes 

24  in 

ARvalve 

0.7 

Mvalve 

0.8 

Equivalence  Ratio 

0.9012 

Frequency 

64.906  Hz 

Purge  Fraction 

0.5 

Fuel  Heating  Value 

18400  BTU/lbm 

Cp  Compressor 

0.2415  BTU/lbm  R 

gamma  c 

1.3986 

Cp  Turbine 

1.2957 

Bleed  Air 

1% 

Coolant  Air 

5% 

Coolant  Air  2 

5% 

Diffuser  Max  Press  Ratio 

0.995 

Burner  Press  Ratio 

0.96 

Nozzle  Press  Ratio 

0.985 

Fan  Press  Ratio 

0.98 

Fan  Adiabatic  Eff 

0.8827 

LPC  Adiabatic  Eff 

0.8827 

HPC  Adiabatic  Eff 

0.8573 

HPT  Adiabatic  Eff 

0.9057 

LPT  Adiabatic  Eff 

0.9084 

Burner  Eff 

0.995 

Mech  Shaft  LP  Spool  Eff 

0.99 

Mech  Shaft  HP  Spool  Eff 

0.99 

Max  Tt4 

3200  R 

Max  Comp  Press  Ratio 

32 

Max  Press  at  station  3 

650  psia 

Max  Temp  at  station  3 

1860  R 

Max  %  Ref  RPM  -  LP  Spool 

110 

Max  %  Ref  RPM  •  HP  Spool 

110 
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The  inputs  in  Table  3.5  are  identical  to  those  of  the  baseline  engine  in  NPSS,  with 
the  exception  of  the  new  parameters  added  for  the  PDC  section:  tube  dimensions,  number 
of  tubes,  ARyaive,  M valve,  equivalence  ratio,  frequency,  and  purge  fraction.  ARvaiVe  is  the 
ratio  of  the  inlet  valve  cross-sectional  area  to  the  detonation  tube  cross-sectional  area  and 
M valve  is  the  Mach  number  into  the  valve.  The  values  of  these  parameters  are  the  results 
of  a  parametric  study  performed  for  optimal  engine  performance  at  SLS.  Altitude,  Mach 
number,  and  fill  fraction  are  the  parameters  varied  to  analyze  the  hybrid  engine 
performance  off-design.  The  NPSS  PDC  burner  element  (BRN36)  configuration  is  shown 
in  Fig.  3.2. 


3.1 


•MR40 


Figure  3.2  Pulsed  Detonation  Combustor  (BRN36)  configuration  with  station 
numbers 
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Flow  leaves  the  high  pressure  compressor  and  enters  the  pulsed  detonation 
combustor.  The  flow  control  at  the  detonation  tube  inlet  is  modeled  as  a  pressure  loss 
AP  IP  term  between  the  inlet  and  detonation  tubes.  This  pressure  loss  matches  the  dry- 
duct  pressure  loss  experienced  by  the  conventional  combustor.  For  the  flow  going  into 
the  tubes,  this  dry  duct  pressure  loss  is  intended  to  represent  pressure  loss  through  a 
valve.  Since  the  detonation  engine  operates  at  a  higher  equivalence  ratio  than 
conventional  engines,  it  requires  less  air  to  mix  with  the  fuel  for  a  similar  enthalpy 
generation.  Balancing  the  mass  flow  through  the  tubes  necessitates  shunting  some  of  the 
air  around  the  detonation  tubes  through  an  internal  bypass.  The  mass  flow  rate  through 

the  internal  bypass  1  ni[BF  1  is  defined  as: 

mOP  ~  mccn_ar  ~ m»ibe _ar  (3 .2) 

f72 

where  con  _ar  is  the  combined  mass  flow  rate  entering  the  PDC  before  it  is  split.  The 
internal  bypass  ratio  (iBPR)  equals: 

iBPR=  m:BF  (3.3) 

In  order  to  determine  the  mass  flow  rate  into  the  detonation  tubes  '  )=  the 

fill  fraction  iff)  and  purge  fraction  (pf)  must  first  be  defined.  The  fill  air  is  the  air  mixed 
with  the  fuel,  while  the  purge  air  is  the  unmixed  portion.  The  purge  air  is  used  to  expel 
the  burned  gases  and  also  serves  to  cool  the  detonation  tubes  between  cycles.  The  fill  and 
purge  fractions  are  defined  in  terms  of  volume  of  air  for  their  respective  portions  of  the 
tube  filling  process: 
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(3.4) 


ff  = 


Pf  = 


V , 


fuel -air  mix 

A  */ 

^ tube  ’'tube 


V 

purge  air 

A  *  / 

iiibe  L tube 


(3.5) 


Since  the  air  is  stopped  in  the  tube  when  the  valve  is  closed,  total  density  ( pt )  can  be 

found.  The  purge  and  fuel-air  masses  are  calculated  by  multiplying  the  total  density  by 
the  tube  volume  and  the  purge  fraction  and  fill  fractions,  respectively: 


m 


ir  Pf  '  ^ tube  ’  Pt 


purge  air 


^ fuel— air  mix  ff  ^ 'tube  P t 


(3.6) 


(3.7) 


where  Vtube  Afube  *  ltube- 

Equations  3.6  and  3.7  represent  the  amount  of  purge  air  and  fuel-air  mixture  that 
flow  into  one  tube  during  each  cycle.  The  amount  of  air  to  send  through  the  valve  at  the 
opening  of  the  detonation  tube  is  calculated  in  Eqs.  3.8  -  3.12: 


m fuel-air  mix  ~  m fill  air  +  171  fuel 


171  fuel-air  mix  _  m  fill  air  m fuel 
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fill  air  fill  air  fill  air 
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ff*Kube*P, 
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(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 
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where  men  air  is  the  mass  of  the  air  that  is  detonated  and  mtuci  is  the  fuel  used  during  one 
cycle.  Once  the  mass  of  air  flowing  into  the  tubes  during  one  cycle  is  known,  the  time 

averaged  steady  state  mass  flow  rate  into  the  detonation  tubes  1  ;  is  calculated  by 

multiplying  the  total  mass  of  air  into  the  tubes  by  the  user  inputs  of  frequency  if)  and  the 
number  of  tubes  (n tubes): 


K +  ™purp  _*r)  *  / 

(3  - 1 3> 

The  thermodynamics  within  the  pulsed  detonation  combustor  are  modeled  after 
the  work  of  Heiser  and  Pratt  (Heiser  and  Pratt,  2002:1)  with  a  Dyer  and  Kaemming 
correction  (Dyer  and  Kaemming,  2002:1)  to  more  accurately  conserve  system  energy.  To 
calculate  the  detonation  properties  at  the  tube  exit,  implementation  into  NPSS  required  a 


few  modifications  to  Equations  2.7  and  2.8.  The  quantity  -  as  defined  by  Eqs.  2.11 
and  2. 12  is  rearranged  as: 


0  =  V^ 

J  C/3 


(3.14) 


This  allows  for  burner  inlet  and  exit  enthalpies,  specific  heat,  and  inlet  engine 
temperature  parameters  to  be  used  to  solve  for  the  Chapman- Jouguet  Mach  number,  Eq. 
2.8,  and  entropy  gain  across  the  burner,  Eq.  2.9.  The  pressure  rise  across  the  shock  is  then 
calculated  as: 


Pa  _  1  +  Y^cj  Pi  >  , 

p0  r+ 1  p0 


(3.15) 
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Having  solved  for  the  pressure  and  entropy  at  the  detonation  tube  exit,  all  properties  of 
the  fluid  at  this  point  (station  4)  are  solved  for  through  thermodynamic  relationships. 
However,  as  mentioned  in  the  previous  chapter,  instead  of  using  the  pressure  and  entropy 
at  the  detonation  tube  exit  to  solve  for  all  the  properties  of  the  fluid  at  this  point,  which 
ignores  the  eventual  pressure  loss  that  the  gas  will  go  through  due  to  expansion  waves, 
entropy  and  the  change  in  enthalpy  liberated  by  the  combustion  process  are  used  to  solve 
for  the  properties  at  the  detonation  tube  exit.  Available  energy  is  calculated  using  the 
known  CJ  entropy  originally  calculated  by  Heiser  and  Pratt  (Heiser  and  Pratt,  2002: 1), 
Eq.  2.9,  with  the  known  system  enthalpy  of  (h0+qadd)  in  order  to  ensure  that  energy  is 
conserved  within  the  system.  The  NPSS  PDC  burner  element  code  is  found  in  Appendix 
D. 

Rasheed  et  al.  (Rasheed  et  ah,  2006)  showed  that  exhausting  a  pulsed  detonation 
combustor  directly  into  a  turbine  lowers  the  turbine  efficiency  and  has  structural 
ramifications  affecting  the  engine  life.  The  hybrid  model  is  based  on  the  assumption  that 
the  flow  into  the  turbine  is  steady  flow.  A  subelement  to  the  PDC  was  created  that  allows 
for  the  application  of  a  pressure  drop  and  enthalpy  loss,  however,  no  such  loss  is  applied 
in  this  model. 


Hybrid  Turbofan  Engine  Off-Design 

Fill  fraction  is  the  primary  method  of  thrust  control  as  prior  work  on  PDEs 
indicated  that  performance  gains  may  be  made  at  fill  fractions  other  than  one  (Glaser  et 


al.,  2004:1).  Frequency  and  equivalence  ratio  can  also  be  used  as  variables  to  throttle  the 
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engine  as  shown  by  the  work  of  Schauer  et  al.(Schauer  et  ah,  2001 :6)  and  Hoke  et  al. 
(Hoke  et  al.,  2005:2).  Frequency  is  chosen  to  be  a  user  input  and  constant  throughout 
flight,  whereas  for  the  hybrid  engine  to  operate  at  altitudes  above  13,000  ft,  the 
equivalence  ratio  has  to  be  varied  to  gain  the  maximum  thrust  at  a  given  Mach  number 
and  altitude.  Therefore,  the  equivalence  ratio  is  adjusted  at  each  operating  condition  to 
yield  maximum  thrust.  The  engine  constraints  are  shown  at  the  bottom  of  Table  3.2  and 
controls  were  implemented  into  the  solver  to  ensure  the  model  stayed  within  these 
constraints. 

When  running  the  NPSS  solver  off-design,  an  error  was  given  for  the  constant 
area  mixer,  MIX39,  which  combines  the  flow  coming  out  of  the  detonation  tubes  with  the 
internal  bypass  flow  of  the  PDC.  In  design  mode,  the  user  provides  a  Mach  number  for 
the  tube  flow  into  the  mixer,  which  determines  the  primary  entrance  area  of  the  mixer. 
The  area  of  the  internal  bypass  flow  entering  the  mixer  is  determined  by  varying  the  area 
until  the  static  pressure  of  the  two  streams  equal.  This  conserves  energy,  continuity,  and 
momentum  when  mixing  the  flows  exiting  the  tubes  and  the  internal  bypass  flow  into 
one. 

Running  the  model  in  off-design  mode  at  various  fill  fractions,  Mach  numbers, 
and  altitudes,  however,  yielded  a  static  pressure  difference  between  the  two  flows  when 
entering  MIX39.  As  the  fill  fraction  decreases  from  the  design  fill  fraction,  the  mass  of 
the  fill  air  decreases,  thus  there  is  less  mass  flow  entering  the  tubes  and  more  flow 
entering  the  internal  bypass  as  shown  in  Eqs.  3.13  and  3.2,  respectively: 

+"W_<xr)*  f* 
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(3.13) 


(3.2) 


m-J3>=mc*e  xr 


As  the  mass  flow  of  the  internal  bypass  increases,  the  static  pressure  (p)  decreases  for  a 
fixed  area.  The  internal  bypass  flow  behaves  as  pipe  flow,  and  as  mass  flow  increases, 
velocity  increases  (assuming  incompressible  flow),  Eq.  3.16,  and  thus  the  static  pressure 
must  decrease  to  maintain  a  constant  total  pressure  (pt)  shown  in  Eq.  3.17: 


m  =  pAV 


(3.16) 


— pV2  +  p  =  p,  =  const 


(3.17) 


This  yields  a  lower  static  pressure  entering  the  mixer  from  the  internal  bypass  than  from 
the  detonation  tubes.  The  NPSS  MIX39  requires  that  the  incoming  streams  have  equal 
static  pressures  in  order  to  mix.  To  converge  to  a  solution  off-design,  the  static  pressure 
of  the  internal  bypass  flow  entering  the  mixer  must  be  increased  to  match  the  static 
pressure  of  the  flow  coming  out  of  the  detonation  tubes  into  the  mixer.  This  is 
accomplished  by  decreasing  the  mass  flow  of  the  internal  bypass,  as  mass  flow  and  static 
pressure  have  an  inverse  relationship  as  seen  from  Eqs.  3.16  and  3.17.  To  decrease  the 
mass  flow  of  the  internal  bypass,  a  bleed  was  implemented,  shown  in  Fig.  3.2  as  BLD4, 
in  which  the  mass  flow  is  bled  from  the  internal  bypass  until  the  static  pressure  of  the  two 
streams  entering  MIX39  match.  The  bleed  flow  then  enters  back  into  MIX40  and  MIX44 
on  either  side  of  the  low  pressure  turbine.  The  flow  is  equally  divided  into  M1X40  and 
MIX44  and  does  not  present  a  static  pressure  error.  MIX39  caused  an  NPSS  convergence 
error  because  there  is  an  unequal  amount  of  independent  (9)  and  dependent  variables 
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(10).  The  internal  bleed  is  added  as  an  independent  variable  to  be  varied  until  the  static 
pressure  of  the  two  flows  entering  the  mixer  equal. 

MIX39  allows  for  the  flow  exiting  the  PDC  to  cool  to  below  the  temperature 
constraint  before  entering  the  turbine,  however,  it  has  the  potential  to  be  the  source  of  a 
significant  pressure  loss.  No  pressure  loss  term  is  applied  to  the  mixer  in  this  model,  but 
the  possibility  of  such  a  loss  is  recognized. 

The  next  chapter  utilizes  this  PDC  combustor  configuration  to  analyze  the  hybrid 
turbofan  performance  at  off-design  conditions  as  compared  to  the  conventional  baseline 
model. 
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IV.  Analysis  and  Results 


Introduction 

This  chapter  contains  the  performance  of  the  baseline  turbofan  evaluated  at 
off-design  conditions  in  AEDsys  and  NPSS  and  an  off-design  analysis  of  the  hybrid 
turbofan  modeled  in  NPSS.  The  baseline  turbofan  is  evaluated  in  two  programs  to 
establish  a  foundation  for  comparing  the  engine  off-design. 

Baseline  Turbofan  Off-Design  Performance 

The  on-design  baseline  model  yields  a  thrust  variation  of  approximately  1.0%  and 
thrust  specific  fuel  consumption  (TSFC)  variation  of  approximately  0.8%  between 
AEDsys  and  NPSS.  This  deviation  is  due  to  differences  in  the  thennodynamic  models  of 
the  two  programs  in  specific  heat  and  enthalpy.  The  baseline  turbofan  model  is  run  off- 
design  in  both  AEDsys  and  NPSS.  Throttle  hooks  are  shown  in  Fig.  4. 1  across  six 
different  Mach  number  and  altitude  levels.  The  throttle  hooks  are  generated  by  varying 
the  fuel-air  ratio  (FAR)  to  match  a  selected  a  thrust  value,  which  is  then  plotted  against 
the  corresponding  TSFC.  As  seen  in  the  figure,  the  two  programs  display  agreement  at 
the  design  point;  however,  the  solutions  diverge  at  higher  Mach  numbers  and  altitudes. 
This  divergence  is  due  to  differences  in  off-design  component  efficiencies.  The  AEDsys 
program  component  adiabatic  efficiencies  do  not  change  from  their  on-design  values, 
whereas  the  NPSS  program  utilizes  component  maps  for  off-design  performance. 
Variation  in  efficiencies  for  NPSS  at  off-design  conditions  is  shown  in  Figs.  4.2  and  4.3. 
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The  baseline  turbofan  is  run  at  maximum  thrust  at  maximum  Tt4  at  SLS,  and  at  maximum 


kc  at  all  other  altitudes.  The  AEDsys  low  pressure  and  high  pressure  spool  efficiencies 

remain  constant  at  0.99  as  expected.  The  NPSS  low  pressure  compressor  and  turbine  and 
high  pressure  compressor  efficiencies  change  as  NC/(NC  design)  increases.  The  high 
pressure  turbine  efficiency  remains  constant;  however,  the  fan  efficiency  experiences  a 
severe  drop  as  it  moves  away  from  the  design  speed.  The  variations  for  the  NPSS 
efficiencies  are  expected  and  may  account  for  the  differences  in  the  throttle  hook  results 
between  AEDsys  and  NPSS. 


Baseline  Throttle  Hook  Comparison  NPSS  and  AEDsys 


Figure  4.1  Throttle  hook  baseline  engine  comparison  using  NPSS  and  AEDsys 
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Baseline  Turbofan  LP  Spool  Adiabatic  Efficiencies  NPSS  and  AEDsys 


Nc/(Nc  design) 


-X—  eta  f  NPSS:  M  =  0.0,  Alt  =  0.0ft 

-X—  eta  f  NPSS:  M  =  0.4,  Alt  =  16  kft 

-X—  eta  f  NPSS:  M  =  0.8,  Alt  =  32  kft 

eta  cL  NPSS:  M  =  0.0,  Alt  =  0.0  ft 
eta  cL  NPSS:  M  =  0.4,  Alt  =  16  kft 
eta  cL  NPSS:  M  =  0.8,  Alt  =  32  kft 
eta  tL  NPSS:  M  =  0.0,  Alt  =  0.0  ft 
eta  tL  NPSS:  M  =  0.4,  Alt  =  16  kft 
eta  tL  NPSS:  M  =  0.8,  Alt  =  32  kft 
eta  LP  Spool  AEDsys 


Figure  4.2  Low  pressure  spool  adiabatic  efficiencies  for  NPSS  and  AEDsys  off-design 


Baseline  Turbofan  HP  Spool  Adiabatic  Efficiencies  NPSS  and  AEDsys 


-X—  eta  cH  NPSS:  M  =  0.0,  Alt  =  0.0  ft 
-X—  eta  cH  NPSS:  M  =  0.4,  Alt  =  16  kft 
-X—  eta  cH  NPSS:  M  =  0.8,  Alt  =  32  kft 
eta  tH  NPSS:  M  =  0.0,  Alt  =  0.0  ft 
eta  tH  NPSS:  M  =  0.4,  Alt  =  16  kft 
eta  tH  NPSS:  M  =  0.8,  Alt  =  32  kft 
eta  HP  Spool  AEDsys 


Nc/(Nc  design) 


Figure  4.3  High  pressure  spool  adiabatic  efficiencies  for  NPSS  and  AEDsys  off-design 
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Hybrid  Turbofan  Off-Design  Results 

The  hybrid  turbofan  model  is  evaluated  off-design  over  a  range  of  Mach  numbers, 
altitudes,  and  fill  fractions  to  detennine  the  effects  on  engine  performance.  Engine 
constraints  are  placed  in  NPSS  to  ensure  that  the  model  does  not  violate  the  maximum 
engine  control  values  of  Tt4,  nc ,  Pt3  and  Tt3,  and  NC/(NC  design)  as  listed  as  the  bottom  of 
Table  3.2. 

Code  Verification  and  Operating  Limit 

The  hybrid  engine  is  run  in  NPSS  at  the  design  point  in  off-design  mode  to 
validate  the  off-design  code  for  the  engine.  Table  4. 1  shows  the  engine  data  at  the  design 
point  in  both  design  and  off-design  mode.  The  data  is  very  similar  with  a  TSFC  and 
thrust  variation  of  0.26%  and  0.01%,  respectively. 

The  engine  operating  envelope  is  found  for  both  the  baseline  and  hybrid  engines 
in  NPSS  and  shown  in  Fig.  4.4.  The  envelope  is  attained  by  running  the  model  at  sea 
level  at  increasing  flight  Mach  numbers  until  NPSS  no  longer  converges  to  a  solution. 
The  far  right  boundary  is  obtained  by  gradually  increasing  the  flight  altitude  at  the 
maximum  Mach  number.  The  top  boundary  is  found  via  a  similar  method.  The  models 
are  run  off-design  at  maximum  thrust  and  a  fill  fraction  of  1 .0  using  the  on-design 
parameters  shown  in  Table  3.5  with  the  exception  of  the  equivalence  ratio.  The  design 
equivalence  ratio  only  allows  the  hybrid  engine  to  operate  at  a  maximum  of  13,000  ft.  To 
obtain  maximum  perfonnance,  the  equivalence  ratio  is  varied  to  yield  the  maximum 
thrust  at  the  flight  Mach  number  and  altitude.  This  is  accomplished  in  NPSS  by  varying 
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the  FAR,  which  is  a  user  input,  to  yield  the  maximum  thrust  at  a  particular  flight 
condition.  Thus,  the  equivalence  ratio  varies  along  the  operating  limit  line  ranging  from  a 
minimum  of  0.79  at  M  =  1.7  and  Alt  =  75,000  ft  to  a  maximum  of  0.93  at  M  =  1.8  and 
Alt  =  0.0  ft.  The  baseline  engine  has  a  higher  altitude  limit  of  1 13,000  ft  as  compared  to 
that  of  the  hybrid  engine  with  a  maximum  of  80,000  ft.  Both  models  have  the  same  Mach 
number  limit  of  2.2. 


Table  4.1  Hybrid  engine  test  data  showing  the  design  point  and  the  design  point  run  at  off-design 
using  NPSS 

^icp  npss_1.  6. 5  -  Rev:  —  model:  pdc  hbtf 

run  by:  caitlin  sol ut i onMode=  steady_state  converge*  1 

case:  2  time:  0.000  timestep: 0.0500  therm_package:GasTbl 

iter/pas/Jac/Broy*  4/  4/  0/  2  run:  12/08/09  13:10:32 


Design  values 

Pidmax*  0.99500  Pi  b  =  1.46943  Eta 
Eta  f  =  0.88270  Eta  cL=  0.88270  Eta  cH= 
e  f  =  0.89210  e  CL  =  0.89210  e  CH= 
gam 2  =  1.39978  gam25=  1.39963  gam3= 
Eta  mL=  0.99000  Eta  mH=  0.9900 

PTO  L  =  0. OkW  PTO  I 

Bleed  =  1.00%  Cool  1=  1.00%  Cool 

**  Thrust  scale  does  not  apply... 

CASE 

converged?  (0  -  no,  1  -  yes) 

Mach  Number  at  0 
Temperature  at  0 
pressure  at  0 
Altitude  at  0 
Total  Temp  at  4 
Pi  r 

/  Tau  r 

Pi  d 
Pi  f 

/  Tau  f 

Pi  CL 

/  Tau  CL 
Pi  CH 

Tau  CH 

Tau  ml 
Pi  tH 

Tau  tH 

Tau  m2 
Pi  tL 

/  Tau  tL 

LP  spool  rpm  (%  reference  pt) 

HP  spool  rpm  (%  reference  pt) 

Alpha 

Ptl9/P19 

P0/P19 

Mach  Number  at  19 

Pt9/P9 

P0/P9 

Mach  Number  at  9 
Mass  Flow  Rate  at  0 
Corr  Mass  Flow  at  0 
Flow  Area  at  0  (ft2) 

Flow  Area*  at  0  (ft2) 

Flow  Area  at  8  and  18  (ft2) 

MB  -  Fuel /Air  Ratio  (f) 

Overall  Fuel /Air  Ratio  (fo) 
specific  Thrust  (F/mO) 

Tnrust  spec  Fuel  consumption  (s) 

Thrust  (Fn) 

Fuel  Flow  Rate  (Ibm/hr) 
propulsive  Efficiency  (%) 

Thermal  Efficiency  (%) 

overall  Efficiency  (%) 


b  =  0.99500 

Pi 

n  =  0.98500 

0.85730  Eta 

t  H= 

0.90570  Eta  tL= 

0.90840 

0.89982  e 

t  H= 

0.89259  e  tL= 

0.89357 

1.35506  gam4= 

1.27739  gam4  5= 

gam 5=1. 31906 

1.29762 

1  =  106. 8kW 

hPr 

18400 

2=  1.00% 

Pi 

nf  =  0.98000 

0:  DESIGN  CASE 

1:  OFFDESIGN 

1 

1 

0.0100 

0.0100 

518.6700 

518.6700 

14.6960 

14.6960 

0.0000 

0.0000 

3648.0 

3648. 3 

1.0001 

1.0001 

1.0000 

1.0000 

0.9950 

0. 9950 

1. 5600 

1. 5600 

1.1534 

1.1534 

1. 5600 

1. 5600 

1.1534 

1.1534 

16.6667 

16. 7164 

2. 3730 

2. 3750 

0.9953 

0.9968 

0.2842 

0. 2838 

0.7783 

0.7780 

0.9964 

0.9980 

0.2553 

0.  2553 

0.7503 

0.7503 

100.000 

99. 987 

100.000 

100.04  5 

8.000 

8.002 

1. 5212 

1. 5212 

1.0000 

1.0000 

0.8003 

0.8003 

1. 8524 

1. 8524 

0.8664 

0.8666 

1.0000 

1.0000 

1500.0000 

1499. 9210 

1499.9263 

1499.8473 

1757. 3935 

1757. 3010 

30.3652 

30. 3636 

22.7472 

22.7472 

0.06161 

0.06161 

0.00301 

0.00300 

32.1463 

32.1439 

0. 3410 

0. 3401 

47699.0467 

47692.9194 

16264. 8674 

16220.1770 

2.0466 

2.0465 

40.2805 

40. 3875 

0. 8244 

0. 8265 
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Operating  Envelope 


—  —  —  Baseline  Engine 
Hybrid  Engine 


Flight  Mach  Number 


figure  4.4a  Mach  number  and  altitude  operating  envelope  at  maximum  thrust, 
baseline  and  hybrid  engine  using  NPSS 


Figure  4.4b  shows  the  aircraft  flight  operating  envelope  as  compared  to  the 
operating  envelope  of  the  engine.  The  flight  operating  envelope  is  estimated  using  a  lift 
to  drag  ratio  at  cruise  conditions  on  the  order  of  10  with  each  engine  supporting  100,000 
lbs  of  weight.  This  comes  from  an  assumed  thrust  to  weight  ratio  of  0.4  for  a  40,000  lb 
engine.  The  flight  envelope  is  much  smaller  than  that  of  the  engine,  with  a  maximum 
altitude  of  38,000  ft  and  a  maximum  Mach  number  of  1.83.  The  hybrid  engine’s  flight 
altitude  is  constrained  to  lower  than  that  of  a  conventional  aircraft;  however,  it  is  still 
acceptable  for  flight.  The  lower  altitude  limit  of  the  hybrid  engine  is  due  to  the  limitations 
of  the  internal  bypass.  The  internal  bypass  ratio  is  on  the  order  of  0.3  at  a  fill  fraction  on 
1.0.  This  ratio  increases  as  the  fill  fraction  is  throttled  to  lower  values  as  seen  in  Fig.  4.4c. 
Figure  4.4c  is  generated  by  selecting  the  fill  fraction  and  varying  the  FAR  upward  to 
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increase  thrust.  At  the  FAR  for  maximum  thrust,  any  more  increase  in  FAR  results  in  non 
convergence  in  NPSS  due  to  the  internal  bleed  air  equaling  zero  and  the  static  pressure  of 
the  internal  bypass  flow  entering  MIX39  no  longer  equaling  the  static  pressure  of  the 
flow  exiting  the  detonation  tubes. 


Operating  Envelope  for  the  Hybrid  Engine 


Flight  Mach  Number 


Engine  operating 
envelope 

—  —  —  Flight  operating 
envelope 


figure  4.4b  Mach  number  and  altitude  engine  and  flight  operating  envelope  for  the 
hybrid  engine  at  maximum  thrust  using  NPSS 
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Figure  4.4c  Internal  bypass  ratio  variation  with  fill  fraction  at  SLS  and  cruise 
conditions  at  maximum  thrust 

Component  Data 

Component  interface  data  for  the  baseline  and  hybrid  engines  at  sea  level  static 
and  cruise  conditions  are  shown  in  Tables  4.2  and  4.3.  The  models  are  run  at  the 
configuration  shown  in  Table  3.5  at  maximum  thrust  and  a  fill  fraction  of  1.0.  At  SLS, 
the  hybrid  engine  has  a  12.7%  greater  thrust  than  the  baseline  engine,  but  at  a  cost  of  a 
5.3%  increase  in  TSFC.  At  cruise  conditions  of  M  =  0.8  and  Alt  =  30,000  ft,  the  hybrid 
engine  has  a  thrust  gain  of  14.7%  over  the  baseline  and  it  also  has  a  lower  TSFC  by 
3.1%.  These  results  are  summarized  in  Table  4.3c  and  indicate  that  the  hybrid  engine 
could  yield  better  performance  than  the  baseline  at  cruise  conditions. 
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Table  4.2a  NPSS  component  interface  data  for  the  baseline  engine  at  SLS,  maximum  thrust 


ft ****************************************************************** ************ 

ncp  npss_1. 6. 5  -  Rev:  —  model : Baseline  hbtf  >  run  by:  caitlin 

time:  0.000  timestep: 0.0500  therm_package:GasTbl  Mode:  DESIGN 

Summary  output  Data 

mn  alt  dTs  w  Fg  Fn  tsfc  wfuel  war 
0.010  0.0  0.00  1500.0  42189.2  41668.8  0.3238  13491.12  0.0000 


solutionMode=  steady_state  converge=  1  case:  0 
iter/pas/Jac/Broy=  5/  7/  1/  3  run:  01/07/10  11:28:56 


OPR 

25.870 


INPUT  FLOW 


w 

gamt 

Pt 

Tt 

PS 

TS 

MN 

V 

Aphy/144. 

FAR 

RS 

imp 

FLO 

INLET. Fl_I 

1500.00 

1.39978 

14.697 

518.68 

14.696 

518.67 

0. 0100 

11.16 

1757.394 

0.0000 

0.06856 

3719559.18 

FLl 

SPLIT. Fl_I 

1500.00 

1.39978 

14.623 

518.68 

12.329 

493.97 

0. 5000 

544.91 

40.864 

0.0000 

0.06856 

97950.38 

FLb2 

Fan21.Fl_I 

1333.33 

1.39978 

14.623 

518.68 

12.329 

493.97 

0.  5000 

544.91 

36.324 

0.0000 

0.06856 

87067.00 

FLb3 

Bypassl3.  Fl_I 

1333.33 

1.39963 

22.812 

598.26 

19.286 

569.80 

0. 5000 

585.02 

24.947 

0.0000 

0.06856 

93528.29 

FLb7 

NOZ18. Fl_I 

1333.33 

1.39963 

22.812 

598.26 

19.286 

569.80 

0. 5000 

585.02 

24.947 

0.0000 

0.06856 

93528.29 

FLb8 

Nozsinkl9. Fl_l 

1333.33 

1.39963 

22.356 

598.26 

14.696 

530.38 

0.8003 

903.41 

19.734 

0.0000 

0.06856 

79200.91 

FL2 

LPC20.Fl_I 

166.67 

1.39978 

14.623 

518.68 

12.329 

493.97 

0. 5000 

544.91 

4.540 

0.0000 

0.06856 

10883.38 

FL25 

HPC25. Fl_I 

166.67 

1.39963 

22.812 

598.26 

19.286 

569.80 

0. 5000 

585.02 

3.118 

0.0000 

0.06856 

11691.04 

FL3 

BLD3. Fl_I 

166.67 

1.35513 

380.206 

1419.66 

357.634 

1397.08 

0.3015 

543.78 

0.444 

0.0000 

0.06856 

25662.10 

FL31 

BRN36. Fl_I 

148.33 

1.35513 

380.206 

1419.66 

357.634 

1397.08 

0.3015 

543.78 

0.395 

0.0000 

0.06856 

22839.27 

FL4 

MIX40.  Fl_I 

152.08 

1.28210 

364.998 

2900.00 

0.000 

0.00 

0. 0000 

0.00 

0.000 

0.0253 

0.00000 

0.00 

FL41 

HPT41. Fl_I 

160.41 

1.28490 

364.998 

2830.35 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0239 

0.00000 

0.00 

FL44 

MIX44. Fl_I 

160.41 

1.30753 

86.673 

2112.70 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0239 

0.00000 

0.00 

FL45 

LPT45.Fl_I 

168.75 

1.30942 

86.673 

2080.75 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0227 

0.00000 

0.00 

FL7 

NOZ8. Fl_I 

168.75 

1.33650 

17.608 

1472.58 

13.929 

1387.80 

0.6000 

1072.95 

5.810 

0.0227 

0.06861 

17280.63 

BL  3 

sink39.Fl_I 

1.67 

1.35513 

380.206 

1419. 66 

0.000 

0.00 

0. 0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FL8 

Nozsink9.Fl_l 

168.75 

1.33650 

17.344 

1472.58 

14.696 

1412.26 

0. 5024 

905.79 

6.638 

0.0227 

0.06861 

18797.76 

BURNERS 

TtOUt 

eff 

dPqP 

LHV 

Wfuel 

FAR 

BRN36 

2900.00 

0.9950 

0.0400 

18400 

3.74753 

0.02526 

Table  4.2b  NPSS  component  interface  data  for  the  hybrid  engine  at  SLS,  maximum  thrust  (ff  =  1) 


NCP  NPSS_1. 6. 5  -  Rev:  —  model:  PDC  HBTF 

time:  0.000  timestep :0. 0500  therm_package:GasTbl 


run  by:  caitlin  solutionwode=  steady_state  converge=  1  case:  0 

Mode:  DESIGN  iter/pas/Jac/Broy=  5/  7/  1/  3  run:  03/06/10  09:14:22 


Summary  output  Data 

mn  alt  dTs  w  Fg  Fn  tsfc  wfuel  war  opr 

0.010  0.0  0.00  1500.0  48246.4  47725.9  0.3408  16264.87  0.0000  38.020 


INPUT  FLOW 


w 

gamt 

Pt 

Tt 

PS 

TS 

MN 

V 

Aphy/144. 

FAR 

RS 

imp 

FLO 

INLET. Fl_I 

1500.00 

1.39978 

14.697 

518.68 

14.696 

518.67 

0.0100 

11.16 

1757.394 

0.0000 

0.06856 

3719559.18 

FLl 

SPLIT. Fl_I 

1500.00 

1.39978 

14.623 

518.68 

13.078 

502.60 

0.4000 

439.70 

48.575 

0.0000 

0.06856 

111974.73 

FLb2 

Fan21. Fl_I 

1333. 33 

1.39978 

14.623 

518.68 

13.078 

502.60 

0.4000 

439.70 

43.177 

0.0000 

0.06856 

99533.09 

FLb3 

Bypassl3. Fl_l 

1333.33 

1.39963 

22.812 

598.26 

20.477 

579.73 

0.4000 

472.07 

29.626 

0.0000 

0.06856 

106920.73 

FLb7 

NOZl8.Fl_I 

1333.33 

1.39963 

22.812 

598.26 

20.477 

579.73 

0.4000 

472.07 

29.626 

0.0000 

0.06856 

106920.73 

FLb8 

NozSinkl9.Fl_I 

1333.33 

1.39963 

22.356 

598.26 

14.696 

530.38 

0.8003 

903.41 

19.734 

0.0000 

0.06856 

79200.91 

FL2 

LPC20. Fl_I 

166.67 

1.39978 

14.623 

518.68 

13.078 

502.60 

0.4000 

439.70 

5.397 

0.0000 

0.06856 

12441.64 

FL25 

HPC25. Fl_I 

166.67 

1.39963 

22.812 

598.26 

20.477 

579.73 

0.4000 

472.07 

3.703 

0.0000 

0.06856 

13365.09 

FL3 

BLD3.Fl_I 

166.67 

1.  35513 

380.206 

1419.66 

341.620 

1380.32 

0.4000 

717.35 

0.348 

0.0000 

0.06856 

20825.70 

FL31 

BRN36. Fl_I 

148.33 

1.35513 

380.206 

1419.66 

357.849 

1397.30 

0. 3000 

541.11 

0.397 

0.0000 

0.06856 

22930.36 

F 1 393 

BLD4 . Fl_I 

0.00 

1.39978 

14.696 

518.67 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

Fl  39 

MIX39. Fl_Il 

115.15 

1.25230 

558.776 

3648.04 

362.681 

3336.85 

0.8500 

2286.33 

0.172 

0.0408 

0.06865 

17160.29 

FT  392  MIX39. Fl_I2 

37.70 

1.35513 

364.998 

1419.66 

362.684 

1417.30 

0.0969 

175.89 

0.310 

0.0000 

0.06856 

16413.70 

FL4 

MIX40. Fl_I 

152.85 

1.27177 

454.794 

3161.21 

424.475 

3114.86 

0.3305 

862.59 

0.482 

0.0305 

0.06863 

33574.07 

FL41 

HPT41. Fl_I 

161.18 

1.27494 

454.794 

3081.06 

363.352 

2934.54 

0.6000 

1523.69 

0.317 

0.0288 

0.06862 

24210.48 

FL44 

MIX44. Fl_I 

161.18 

1.29532 

125.458 

2387.41 

96.133 

2246.04 

0.6500 

1455.77 

0.959 

0.0288 

0.06862 

20572.88 

FL45 

LPT45.  Fl_I 

169. 52 

1.29759 

125.458 

2343.79 

96.089 

2204.00 

0.6500 

1443.34 

0.999 

0.0274 

0.06862 

21426.97 

FL7 

NOZ8.Fl_I 

169.52 

1. 31902 

32.033 

1758.62 

28.857 

1714.71 

0.4000 

789.10 

4.734 

0.0274 

0.06862 

23827.24 

BL  3 

Sink39. Fl_I 

1.67 

1.35513 

380.206 

1419.66 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FL8 

NozSink9.Fl_l 

169.52 

1.31902 

31.552 

1758.62 

17.033 

1511.92 

1.0000 

1859.66 

3.000 

0.0274 

0.06862 

17157.42 

Fl  394 

Nozsinkl.Fl_i 

0.00 

1. 50617 

14.696 

0.00 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

BURNERS 

TtOUt 

eff 

dPqP 

LHV 

wfuel 

FAR 

BRN36 

4443.75 

0.9950 

0. 0400 

18400 

4.51802 

0.06161 
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Table  4.3a  NPSS  component  interface  data  for  the  baseline  engine  at  cruise,  maximum  thrust 


ncp  NPSSJ..6.5  -  Rev:  —  model : Baseline  hbtf  > 

time:  0.000  timestep:0.0500  therm_package:GasTbl 


run  by:  Caitlin  solutionMode=  5TEADY_STATE  converge=  1  case:  1 

Mode:  offdesign  iter/pas/:ac/Broy=  19/  29/  1/17  run:  01/21/10  15:05:18 


mn  alt 
0.800  30000.0 


summary  output  Data 

dTs  w  Fg  Fn  tsfc 

0.00  779.5  30754.6  11483.9  0.6705 


Wfuel  WAR  OPR 

7700.22  0.0000  31.841 


INPUT  FLOW 


w 

gamt 

Pt 

Tt 

Ps 

Ts 

MN 

V 

Aphy/144. 

FAR 

RS 

imp 

FLO 

INLET. Fl_I 

779.47 

1.40010 

6.679 

464.31 

4.364 

411.69 

0.8000 

795.44 

34.249 

0.0000 

0.06856 

40793.83 

FL1 

SPLIT. Fl_I 

779.47 

1.40010 

6.646 

464.31 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FLb2 

Fan21. Fl_l 

685.97 

1.40010 

6.646 

464.31 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FLb3 

Bypassl3. Fl_l 

685.97 

1. 39969 

10.978 

564.78 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FLb7 

N0Z18. Fl_I 

685.97 

1.39969 

10.978 

564.78 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FLb8 

Nozsinkl9.Fl_i 

685.97 

1. 39969 

10.758 

564.78 

5.699 

470.64 

1.0000 

1063. 58 

19. 734 

0.0000 

0.06856 

38871.01 

FL2 

LPC20.Fl_I 

93.49 

1.40010 

6.646 

464.31 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FL25 

HPC25. Fl_I 

93.49 

1.39969 

12.369 

568.03 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FL3 

BLD3.  Fl_I 

93.49 

1.35799 

212.674 

1370.63 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FL31 

BRN36.  Fl_I 

83.21 

1.35799 

212.674 

1370.63 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FL4 

MIX40.  Fl_I 

85.35 

1.28227 

204.167 

2882.66 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0257 

0.00000 

0.00 

FL41 

HPT41.Fl_I 

90.02 

1.28514 

204.167 

2811.73 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0243 

0.00000 

0.00 

FL44 

MIX44.  Fl_I 

90.02 

1.30729 

50.031 

2111.99 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0243 

0.00000 

0.00 

FL45 

LPT45.  Fl_I 

94.70 

1. 30927 

50.031 

2077.93 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0231 

0.00000 

0.00 

FL7 

Noz8. Fl_I 

94.70 

1.34283 

7.178 

1360.36 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0231 

0.00000 

0.00 

BL  3 

sink39.Fl_l 

0.93 

1.35799 

212.674 

1370.63 

0.000 

0.00 

0.0000 

0.00 

0.000 

0.0000 

0.00000 

0.00 

FL8 

Nozsink9.Fl_i 

94.70 

1.34283 

7.070 

1360.36 

4.364 

1201.26 

0.8714 

1456.01 

6.638 

0.0231 

0.06861 

8456.78 

BRN36 


BURNERS 

TtOut  eff  dPqP  LHV 
2882.66  0.9950  0.0400  18400 


Wfuel  FAR 

2.13895  0.02571 


Table  4.3b  NPSS  component  interface  data  for  the  hybrid  engine  at  cruise,  maximum  thrust  (ff  =  1) 


NCP  NPS5_1. 6. 5  -  Rev:  —  model:  POC  HBTF 

time:  0.000  timestep: 0.0500  therm_package:GasTbl 


run  by:  caitlin  solutionMode=  STEADY_STATE  converge=  1  case:  1 

Mode:  offdesign  iter/pas/Jac/Broy=  14/  25/  1/12  run:  03/06/10  09:14:23 


Summary  Output  Data 

mn  alt  dTs  w  Fg  Fn  TSFC  wfuel  war  opr 

0.800  30000.0  0.00  754.2  31818.0  13173.0  0.6495  8556.12  0.0000  44.695 


INPUT  FLOW 


w 

gamt 

Pt 

Tt 

PS 

TS 

MN 

V 

Aphy/144 . 

FLO 

INLET. Fl_I 

754.16 

1.40010 

6.679 

464.31 

4.364 

411.69 

0. 8000 

795.44 

33.137 

FLl 

SPLIT. Fl_I 

754.16 

1.40010 

6.646 

464.31 

5.877 

448.29 

0.4228 

438.76 

48.575 

FLb2 

Fan21.Fl_i 

665.31 

1.40010 

6.646 

464.31 

5.891 

448.58 

0.4188 

434.74 

43.177 

FLb3 

Bypassl3.Fl_I 

665.31 

1.39972 

10.473 

547.22 

9.272 

528. 51 

0.4209 

474.28 

29.626 

FLb7 

NOZ18. Fl_l 

665.31 

1.39972 

10.473 

547.22 

9.272 

528.51 

0.4209 

474.28 

29.626 

FLb8 

Nozsinkl9. Fl _l 

665.31 

1.39972 

10.263 

547.22 

5.441 

456.01 

1.0000 

1046.81 

19.734 

FL2 

LPC20. F 1 _ 1 

88.85 

1.40010 

6.646 

464.31 

5.763 

445.79 

0.4558 

471.77 

5.397 

FL25 

HPC25. Fl_I 

88.85 

1.39971 

11.408 

552.82 

10.142 

534.56 

0.4134 

468.50 

3.703 

FL3 

BLD3. Fl_I 

88.85 

1.35941 

197.081 

1346.97 

177.010 

1309.17 

0.4002 

699.99 

0.348 

FL31 

BRN36. F 1 _ I 

79.07 

1.35941 

197.081 

1346.97 

185.469 

1325.49 

0.3001 

527.96 

0.397 

Fl  393 

BLD4. Fl_I 

0.00 

1.39978 

14.696 

518.67 

0.000 

0.00 

0.0000 

0.00 

0.000 

Fl  39 

MIX39. F 1 _ 11 

62.83 

1.25587 

298.535 

3524.52 

188.177 

3203.63 

0.8744 

2308.19 

0.172 

Fl  392 

MIX39. Fl_I2 

18.62 

1.35941 

189.197 

1346.97 

188.167 

1345.03 

0.0897 

158.93 

0.310 

FL4 

MIX40. Fl_I 

81.45 

1.27401 

239.297 

3083.78 

223.337 

3038.25 

0.3303 

852.12 

0.482 

FL41 

HPT41. F 1 _ I 

85.89 

1.27721 

239.297 

3003.82 

191.170 

2860.12 

0.5994 

1504.12 

0.317 

FL44 

MIX44.Fl_I 

85.89 

1.29731 

67.055 

2331.29 

52.100 

2199. 58 

0.6320 

1401.79 

0.959 

FL45 

LPT45. Fl_I 

90.34 

1.29954 

67.055 

2286.95 

52.095 

2156.95 

0.6316 

1388.43 

0.999 

FL7 

NOZ8.Fl_I 

90.34 

1.32172 

16.788 

1704.22 

15.123 

1661.33 

0. 3999 

777.40 

4.734 

BL  3 

sink39.Fl_i 

0.89 

1.35941 

197.081 

1346.97 

0.000 

0.00 

0.0000 

0.00 

0.000 

FL8 

Nozsink9.Fl_I 

90.34 

1.32172 

16.537 

1704.22 

8.920 

1463.35 

1.0000 

1831.54 

3.000 

Fl  394 

Nozsinkl. Fl _l 

0.00 

1.50617 

14.696 

0.00 

0.000 

0.00 

0. 0000 

0.00 

0.000 

BURNERS 

TtOut 

eff 

dPqP 

LHV 

wfuel 

FAR 

BRN36 

4304.85 

0.9950 

0.0400 

18400  2. 

37670  0. 

05918 

FAR 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0. 0000 
0.0000 
0.0000 
0.0393 
0. 0000 
0.0301 
0.0285 
0.0285 
0.0270 
0.0270 
0.0000 
0.0270 
0.0000 


RS 

0.06856 

0.06856 

0.06856 

0.06856 

0.06856 

0.06856 

0.06856 

0.06856 

0.06856 

0.06856 

0.00000 

0.06864 

0.06856 

0.06863 

0.06862 

0.06862 

0.06862 

0.06862 

0.00000 

0.06862 

0.00000 


imp 
39469.29 
51394.06 
45615.04 
49361.74 
49361.74 
37108.89 
5781.91 
6702.35 
10798. 37 
11889.13 
0.00 
9165.29 
8500.78 
17666.07 
12737.20 
10939.36 
11392.06 
12491.16 
0.00 
8996.48 
0.00 
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Table  4.3c  Maximum  power  output  comparison  baseline  and  hybrid  engines 


Baseline  engine 

Hybrid  engine 

Alt  =  0.0  ft,  M  =0.0 

Thrust  (Ibf) 

41668.8 

47725.9 

TSFC (1/h) 

0.3238 

0.3408 

Mass  flow  rate  (Ibm/s) 

1500 

1500 

Tt4(R) 

2900.0 

3161.2 

Alt  =  30  kft,  M  =  0.8 

Thrust  (Ibf) 

11483.9 

13173.0 

TSFC (1/h) 

0.6705 

0.6495 

Mass  flow  rate  (Ibm/s) 

779.5 

754.2 

Tt-  ( R) 

2882.7 

3083.8 

Table  4.4  shows  a  comparison  of  the  fluid  properties  for  both  the  baseline  and 
hybrid  engine  at  the  entrance  and  exit  of  the  combustor  at  SLS  and  cruise  conditions  at 
maximum  thrust.  The  properties  of  the  flow  entering  station  3.1  are  fairly  similar  for  the 
baseline  and  hybrid  engines.  The  properties  of  the  flow  exiting  the  PDC  are  taken  at  the 
exit  of  MIX39  (station  4.0  as  shown  in  Fig.  3.4),  which  combines  the  tube  and  internal 
bypass  flows.  With  no  internal  bypass  (and  no  mixer),  the  flow  exiting  the  PDC  tubes 
could  be  much  higher  than  the  3200  R  limit  as  seen  from  Tables  4.2b  and  4.3b  at  station 
3.9.  M1X39  allows  the  flow  to  cool  below  the  temperature  constraint  of  3200R  at  a  cost 
of  a  reduced  stagnation  pressure.  MIX39  is  an  NPSS  element  that,  through 
thennodynamic  analysis,  drops  the  exit  stagnation  pressure  due  to  constant  area  mixing 
(e.g.,  see  Oates  pg.  166).  The  burner  pressure  ratio  of  the  baseline  engine  is  0.96  for  SLS 
or  cruise.  The  burner  pressure  ratio  for  the  PDC  ranges  from  1 .20  at  SLS  to  1 .2 1  at  cruise 
(M  =  0.8,  Alt  =  30,000  ft).  If  the  exit  stagnation  pressure  of  the  flow  is  considered  at  the 
exit  of  the  detonation  tubes  (station  3.9  in  Tables  4.2b  and  4.3b),  the  pressure  ratio  of  the 
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burner  is  1 .47  at  SLS  and  1.51  at  a  cruise.  Thus,  the  internal  mixer  causes  about  a  23.7% 


reduction  in  the  combustor  pressure  ratio  for  SLS  or  cruise. 


Table  4.4  Combustor  Properties  at  SLS  and  Cruise,  Maximum  Power  (ff  =  1) 


Baseline  SLS 

Hybrid  SLS 

Baseline  Cruise 

Hybrid  Cruise 

Altitude 

0.0  ft 

0.0  ft 

30,000  ft 

30,000  ft 

Fit  Mach  number 

0.0 

0.0 

0.8 

0.8 

Eq  Ratio 

0.90 

0.86 

iBPR 

0.342 

0.308 

Mcj 

3.37 

3.42 

Station  3.1  flow  in: 

Tt 

1419.7  R 

1419.7  R 

1370.6  R 

1347.0  R 

Pt 

380.2  psi 

380.2  psi 

212.7  psi 

197.1  psi 

Mass  flow  rate 

148.3  lbm/s 

148.3  lbm/s 

83.2  lbm/s 

79.1  lbm/s 

Mach  number 

0.3 

0.3 

Station  4.0  flow  out: 

Tt 

2900.0  R 

3161.2  R 

2882.7  R 

3038.8  R 

Pt 

365.0  psi 

454.8  psi 

204.2  psi 

239.3  psi 

Mass  flow  rate 

152.1  lb  m/s 

152.9  lbm/s 

85.3  lbm/s 

81.5  lbm/s 

Mach  number 

0.33 

0.33 

The  combustor  properties  in  Table  4.4  are  calculated  assuming  steady  flow 
through  the  PDC.  Applying  a  4%  pressure  loss  in  the  TTSS  element  to  account  for 
unsteady  losses  yields  a  0.5%  reduction  in  thrust  and  a  0.5%  increase  in  TSFC.  An  8% 
loss  would  yield  a  1.2%  reduction  in  thrust  and  a  1.2%  increase  in  TSFC.  These  losses 
may  be  low  due  to  their  application  at  the  exit  of  the  detonation  tubes.  Exit  losses  do  not 
capture  valve  losses.  In  order  to  model  the  pressure  losses  due  to  the  opening  and  closing 
of  the  valves  in  to  the  detonation  tubes,  unsteady  affects  should  be  included  in  the 
thennodynamics  of  the  PDC  cycle.  This  was  not  included  in  this  work. 


47 


Hybrid  Turbofan  Performance  Comparison 

Throttle  hooks  for  the  hybrid  engine  are  perfonned  at  the  configuration  shown  in 
Table  3.5  at  maximum  thrust  at  a  fill  fraction  of  1.0  and  compared  to  that  of  the  baseline 
engine  in  Figs.  4.5a  and  4.5b. 


Figure  4.5a  Throttle  hook  comparison  baseline  and  hybrid  turbofan  engines  in 
NPSS 
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Figure  4.5b  Zoomed  in  view  of  Fig.  4.5a 


Figure  4.5b  shows  a  comparison  of  the  hybrid  and  baseline  engine  throttle  hooks 
at  sea  level  static  and  30,000  ft  at  a  variety  of  Mach  numbers.  The  baseline  engine  yields 
a  lower  TSFC  at  lower  thrust  values  for  a  majority  of  the  thrust  range;  however,  as  the 
thrust  increases,  the  TSFCs  cross.  This  crossing  also  occurs  at  the  higher  Mach  numbers. 
It  is  seen  that  the  baseline  engine  has  a  larger  thrust  range.  The  lower  hybrid  range  of 
thrust  occurs  since  the  internal  bypass  bleed  air  (BLD4)  at  station  3.93  (Fig.  3.2)  is  much 
greater  than  the  air  flowing  through  the  internal  bypass,  due  to  reduced  bypass  flow 
needed  to  balance  the  static  pressure  at  MIX39,  and  the  solver  cannot  converge  on  a 
solution.  The  internal  bypass  bleed  air  variation  with  thrust  is  shown  in  Fig.  4.5c.  As  the 
thrust  increases,  less  air  is  bled  from  the  internal  bypass  to  balance  static  pressures  in 
MIX39.  The  thrust  increases  until  the  bleed  air  is  zero,  after  which  the  static  pressures 
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will  not  balance.  If  it  were  possible  to  provide  additional  air  pumped  into  the  mixer,  one 
could  increase  the  thrust.  No  attempts  were  made  for  this  work.  The  internal  bypass  bleed 
is  the  limiting  factor  for  the  hybrid  engine  performance. 


Figure  4.6  shows  the  effects  of  fill  fraction,  with  throttle  hooks  run  at  sea  level 
static  and  cruise  conditions  at  maximum  thrust  at  the  configuration  shown  in  Table  3.5. 
Figure  4.6  shows  that  fill  fraction  can  be  adjusted  to  reduce  the  TSFC  at  any  thrust  level 
to  roughly  that  of  the  baseline  engine  and  perhaps  slightly  lower.  The  best  results  are  seen 
at  the  maximum  thrust  level  of  each  fill  fraction  shown  in  Figs.  4.6b  and  4.6d. 
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Throttle  Hooks  for  Hybrid  Engine  at  SLS  (various  fill  fractions) 


ff  =  0.4 


-  -ff  =  0.6 


*— ff  =  0.8 
H — ff  =  1.0 
■* —  Baseline 


Thrust  (Ibf) 


Figure  4.6a  Throttle  hooks  at  various  fill  fractions  at  SLS 


Throttle  Hooks  for  Hybrid  Engine  at  SLS  (various  fill  fractions) 


Thrust  (Ibf) 


Hybrid 


Baseline 


Figure  4.6b  Throttle  hooks  at  fill  fractions  from  0.4  to  1.0  at  maximum  thrust  at 
SLS 
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Throttle  Hooks  for  Hybrid  Engine  at  Cruise  (various  fill  fractions) 


Thrust  (Ibf) 


ff  =  0.4 


-  *ff  =  0.6 


ff  =  0.8 


H - ff  =  1.0 

—  Baseline 


Figure  4.6c  Throttle  hooks  at  various  fill  fractions  at  cruise  (M  =  0.8,  Alt  =  30kft) 


Throttle  Hooks  for  Hybrid  Engine  at  Cruise  (various  fill  fractions) 


Thrust  (Ibf) 


Hybrid 


Baseline 


Figure  4.6d  Throttle  hooks  at  fill  fractions  from  0.4  to  1.0  at  maximum  thrust  at 
cruise  (M  =  0.8,  Alt  =  30kft) 
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For  comparison  purposes,  throttle  hooks  were  run  at  various  frequencies  in  Fig. 
4.7  at  the  configuration  shown  in  Table  3.5  at  a  fill  fraction  of  1.0.  The  results  are  very 
similar  and  indicate  that  frequency  may  be  considered  as  an  additional  throttling 
parameter  for  future  research. 


Throttle  Hooks  for  Hybrid  Engine  at  SLS  (various  frequencies) 
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Figure  4.7a  Throttle  hooks  at  various  frequencies  at  SLS 
Throttle  Hooks  for  Hybrid  Engine  at  SLS  (various  frequencies) 
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Figure  4.7b  Throttle  hooks  at  various  frequencies  at  maximum  thrust  at  SLS 
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Throttle  Hooks  for  Hybrid  Engine  at  Cruise  (various  frequencies) 
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Figure  4.7c  Throttle  hooks  at  various  frequencies  at  cruise  (M  =  0.8,  Alt  =  30kft) 


Throttle  Hooks  for  Hybrid  Engine  at  Cruise  (various  frequencies) 


Thrust  (Ibf) 


- Hybrid 

x  Baseline 


Figure  4.7d  Throttle  hooks  at  various  frequencies  at  maximum  thrust  at  cruise 
(M  =  0.8,  Alt  =  30kft) 


54 


Figures  4.8a,  4.8b,  and  4.8c  show  the  effect  of  flight  Mach  number  for  the  hybrid 
engine  as  compared  to  the  baseline  at  various  altitudes  and  fill  fractions  at  sea  level  and 
30,000  ft.  The  model  is  run  at  the  configuration  in  Table  3.5  at  maximum  thrust.  At 
higher  Mach  numbers  the  baseline  engine  thrust  is  slightly  higher,  except  for  a  fill 
fraction  of  one  at  sea  level.  At  cruising  altitude  there  is  a  range  of  fill  fractions  where  the 
hybrid  engine  can  match  the  baseline. 


Thrust  Variation  with  Flight  Mach  Number  (ff  =  1.0,  various  altitudes) 
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Figure  4.8a  Variation  of  thrust  with  flight  Mach  number  comparison  of  the  baseline 
and  hybrid  engines  in  NPSS  at  0.0  ft,  30,000  ft,  and  60,000  ft  (maximum  thrust) 
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Thrust  Variation  with  Flight  Mach  Number  at  Sea  Level  (various  fill  fractions) 
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Figure  4.8b  Variation  of  thrust  with  flight  Mach  number  comparison  of  the  baseline 
and  hybrid  engines  in  NPSS  at  fill  fractions  of  0.6,  0.8,  and  1.0  at  sea  level 
(maximum  thrust) 


Thrust  Variation  with  Flight  Mach  Number  at  30kft  (various  fill  fractions) 
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Figure  4.8c  Variation  of  thrust  with  flight  Mach  number  comparison  of  baseline  and 
hybrid  engines  at  fill  fractions  of  0.6,  0.8,  and  1.0  at  30,000  ft  (maximum  thrust) 
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The  effects  of  fill  fraction  on  thrust  are  plotted  at  SLS  and  cruise  conditions  in 
Fig.  4.8d.  For  the  configuration  listed  in  Table  3.5  and  at  maximum  thrust,  in  general,  a 
lower  fill  fraction  corresponds  to  a  lower  thrust. 


Thrust  Variation  with  Fill  Fraction  for  Hybrid  Engine 


Fill  fraction 


M  =  0.0,  Alt  =  0.0  ft 


M  =  0.8,  Alt  =  30000  ft 


Figure  4.8d  Thrust  variation  with  fill  fraction  for  the  hybrid  engine  at  SLS  and 
cruise  (maximum  thrust) 


The  thrust  is  divided  by  the  freestream  pressure  to  determine  the  effects  of 
altitude  on  the  engine.  As  shown  in  Fig.  4.8e,  the  altitude  has  a  significant  effect  on  thrust 
as  the  curves  at  SLS  and  cruise  conditions  are  nearly  identical. 
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Figure  4.8e  Thrust  divided  by  free  stream  pressure  variation  with  fill  fraction  for 
the  hybrid  engine  as  SLS  and  cruise  (maximum  thrust) 


Fill  fraction  is  also  varied  in  Fig.  4.9  to  detennine  the  effects  on  TSFC  at  SLS  and 
cruise  conditions.  The  model  is  run  at  maximum  thrust  at  the  configuration  in  Table  3.5. 
Results  indicate  that  at  cruise  conditions  the  lowest  TSFC  is  found  at  a  fill  fraction  of  0.7. 
The  two  curves  are  at  different  thrusts.  At  a  thrust  of  13,000  lbs,  TSFC  at  cruise  is  0.6468 
and  TSFC  at  SLS  is  0.2794. 
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Figure  4.9  TSFC  variation  with  fill  fraction  for  the  hybrid  engine  at  SLS  and  cruise 
(maximum  thrust) 


The  mass  flow  rate  variation  with  thrust  is  shown  in  Figs.  4. 10a,  4. 10b,  and  4. 10c. 
The  model  is  run  at  the  configuration  in  Table  3.5  at  maximum  power  at  SLS  and  cruise 
conditions  and  fill  fractions  of  0.6,  0.8,  and  1.0.  The  baseline  engine  has  a  higher  mass 
flow  rate  at  all  thrusts  than  the  hybrid.  The  fill  fraction  affects  the  mass  flow  rate  only  by 
the  range  of  thrust  it  covers. 
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Figure  4.10a  Mass  flow  rate  variation  with  thrust  comparison  of  the  baseline  and 
hybrid  engines  at  SLS  and  cruise  conditions  (maximum  thrust) 


Mass  Flow  Rate  Variation  with  Thrust  Comparison  at  SLS  (various  fill  fractions) 
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Figure  4.10b  Mass  flow  rate  variation  with  thrust  comparison  of  the  baseline  and 
hybrid  engines  at  SLS  at  fill  fractions  of  0.6,  0.8,  and  1.0  (maximum  thrust) 
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Mass  Flow  Rate  Variation  with  Thrust  Comparison  at  Cruise  (various  fill 

fractions) 
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Figure  4.10c  Mass  flow  rate  variation  with  thrust  comparison  of  the  baseline  and 
hybrid  engines  at  30,000  ft  at  fill  fractions  of  0.6,  0.8,  and  1.0  (maximum  thrust) 


Component  Performance 

The  fan,  high  pressure  compressor  (HPC),  low  pressure  turbine  (LPT),  and 
turbofan  bypass  ratio  variation  with  Tt2/To  are  shown  in  Figs.  4.1 1-  4.14.  The  model  is 
run  at  the  configuration  in  Table  3.5  at  maximum  power  at  a  fill  fraction  of  1.0  at  SLS 
and  30,000  ft.  For  the  baseline  model,  the  fan  and  HPC  pressure  ratios  break  at  1.0  at 
30,000  ft.  The  hybrid  model  is  already  past  the  break  point  in  the  plots  shown.  The  LPT 
is  choked  for  both  models.  The  bypass  ratio  is  shown  to  steadily  increase  for  both  the 
baseline  and  hybrid  engines. 
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Figure  4.11  Fan  pressure  ratio  variation  with  TVTo  comparison  of  the  baseline  and 
hybrid  engines  at  SLS  and  cruise  (maximum  thrust) 


HPC  Pressure  Ratio  Variation  with  Tt2/To  Engine  Comparison 
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Figure  4.12  High  pressure  compressor  ratio  variation  with  T^/To  comparison  of  the 
baseline  and  hybrid  engines  at  SLS  and  cruise  (maximum  thrust) 
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LPT  Pressure  Ratio  Variation  with  Tt2/To  Engine  Comparison 
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Figure  4.13  Low  pressure  turbine  pressure  ratio  variation  with  Tt2/To  comparison  of  the 
baseline  and  hybrid  engines  at  SLS  and  cruise  (maximum  thrust) 


Turbofan  Bypass  Ratio  Variation  with  Tt2/To  Engine  Comparison 
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Figure  4.14  Turbofan  bypass  ratio  variation  with  Tt2/To  comparison  of  the  baseline  and 
hybrid  engines  at  SLS  and  cruise  (maximum  thrust) 
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Component  Adiabatic  Efficiencies 

It  was  of  interest  to  determine  how  changing  the  adiabatic  efficiencies  of  various 
components  affected  the  performance  of  the  hybrid  engine.  Changes  in  the  inlet,  fan, 
turbine,  compressor,  and  burner  efficiencies  were  evaluated  via  throttle  hooks  shown  in 
Figs.  4.16  -  4.21,  respectively.  The  dotted  lines  represent  the  on-design  efficiencies.  The 
model  was  run  at  maximum  thrust  at  the  configuration  shown  in  Table  3.5.  With  a  5% 
decrease  in  efficiency,  each  component  resulted  in  an  average  of  1.2%  to  1.5%  increase 
in  TSFC  with  the  exception  of  the  burner.  The  burner  efficiency  is  used  to  calculate  the 
heat  addition  into  the  system  as  shown  in  Eq.  2.11.  Decreasing  the  efficiency  of  the 
burner  in  Fig.  4.21  less  than  97.5%  leads  to  choking,  in  which  case  the  solver  cannot 
converge  to  a  solution.  Decreasing  burner  efficiency  from  its  original  value  of  99.5%  to 


97.5%  resulted  in  a  3.4%  average  decrease  in  TSFC. 


Figure  4.15  Throttle  hooks  for  inlet  efficiency  of  the  hybrid  engine  (maximum 
thrust) 
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Figure  4.16  Throttle  hooks  for  fan  efficiency  of  the  hybrid  engine  (maximum  thrust) 


Figure  4.17  Throttle  hooks  for  high  pressure  turbine  efficiency  of  the  hybrid  engine 
(maximum  thrust) 
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Figure  4.18  Throttle  hooks  for  low  pressure  turbine  efficiency  of  the  hybrid  engine 
(maximum  thrust) 


Figure  4.19  Throttle  hooks  for  low  pressure  compressor  efficiency  of  the  hybrid 
engine  (maximum  thrust) 
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Figure  4.20  Throttle  hooks  for  high  pressure  compressor  efficiency  of  the  hybrid 
engine  (maximum  thrust) 


Figure  4.21  Throttle  hooks  for  burner  efficiency  of  the  hybrid  engine  (maximum 
thrust) 
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Parameter  On-Design  Choices 

The  design  choices  listed  in  Table  3.5  used  for  the  analysis  above  are  based  on  a 
parametric  study  performed  with  the  turbofan  model  at  design  conditions.  Since  the 
current  model  was  modified  for  off-design  performance,  changes  in  these  design 
parameters  may  improve  engine  performance  at  design  and/or  off-design  conditions. 

Throttle  hooks  were  run  at  SLS  and  cruise  conditions  for  various  on-design 
frequencies,  purge  fractions,  and  equivalence  ratios  shown  in  Figs.  4.22  -  4.24, 
respectively,  to  detennine  if  changing  these  parameters  could  improve  performance. 
Figure  4.22  indicates  that  a  lower  design  frequency  may  improve  TFSC,  but  at  the 
expense  of  thrust  range.  A  lower  equivalence  ratio  yields  TSFC  improvements  similar  to 
that  of  frequency;  however,  a  design  equivalence  ratio  below  0.88  chokes  the  burner  inlet 
and  a  solution  cannot  be  converged  in  Fig.  4.23.  The  design  purge  fraction  may  yield  a 
lower  TSFC  as  a  purge  fraction  of  0.75  yields  a  2.2%  decrease  in  TSFC  below  design 
value  of  0.5  with  the  same  thrust  range  shown  in  Fig.  4.24.  As  seen  from  Eq.  3.6,  the 
mass  of  the  purge  air  increases  as  the  purge  fraction  increases.  The  increased  purge 
fraction  increases  the  mass  flow  through  the  tubes  while  decreasing  the  flow  through  the 
internal  bypass.  This  results  in  an  increase  in  thrust  as  well  as  TSFC,  however,  thrust 
increases  more  than  TSFC.  A  thorough  parametric  study  should  be  conducted  to  identify 
such  design  parameters  as  number  of  tubes,  tube  geometry,  etc.,  that  may  yield  better 
perfonnance  at  design  and  off-design  conditions. 
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Figure  4.22  Throttle  hooks  at  various  design  frequencies  at  SLS  and  cruise  (M  =  0.8, 
Alt  =  30,000  ft) 
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Figure  4.23  Throttle  hooks  at  various  design  equivalence  ratios  at  SLS  and  cruise 
(M  =  0.8,  Alt  =  30,000  ft) 
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Throttle  Hooks  Purge  Fraction  Design  Choice  (SLS) 


10000  20000  30000  40000 

Thrust  (Ibf) 


50000 


pf  =  0.25,  M  =0.0,  Alt  =  0.0ft 
pf=0.5,  M  =0.0,  Alt  =  0.0ft 
pf  =  0.75,  M  =0.0,  Alt  =  0.0ft 


60000 


6000 


Throttle  Hooks  Purge  Fraction  Design  Choice  Cruise 


8000 


10000 

Thrust  (Ibf) 


12000 


pf  =  0.25,  M  =  0.8,  Alt  =  30  kft 
pf  =  0.5,  M  =  0.8,  Alt  =30  kft 
pf  =  0.75,  M  =  0.8,  Alt  =  30  kft 


14000 


Figure  4.24  Throttle  hooks  at  various  design  purge  fractions  at  SLS  and  cruise  (M  : 
0.8,  Alt  =  30,000  ft) 
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V.  Conclusions  and  Recommendations 


Introduction 

A  turbofan  engine  with  a  pulsed  detonation  combustor  may  have  performance 
gains  over  a  conventional  turbofan,  specifically  in  the  areas  of  thrust  and  thrust  specific 
fuel  consumption.  Research  (Andrus,  2007)  has  shown  that  at  design  conditions,  the 
hybrid  engine  may  allow  an  8.0%  decrease  in  TSFC  while  maintaining  thrust.  The 
objective  of  this  work  was  to  develop  a  hybrid  engine  model  with  a  pulsed  detonation 
combustor  to  run  off-design  in  NPSS,  and  to  evaluate  the  performance  of  the  hybrid 
engine  at  various  off-design  conditions. 


Hybrid  Turbofan  Engine  Off-Design  Performance 

To  detennine  the  performance  of  the  hybrid  turbofan  engine,  the  model  was  run 
over  a  range  of  off-design  conditions,  at  various  Mach  numbers,  altitudes,  and  fill 
fractions,  and  compared  to  that  of  the  baseline  engine.  Equivalence  ratio,  frequency,  and 
fill  fraction  were  all  potential  parameters  to  be  used  to  throttle  the  hybrid  engine.  After 
perfonning  an  operating  limit  analysis  on  the  hybrid  engine,  it  was  discovered  that  at 
constant  equivalence  ratio,  the  aircraft  has  a  maximum  operating  altitude  of  13,000  ft.  In 
order  for  the  hybrid  engine  to  operate  at  realistic  cruising  altitudes,  the  equivalence  ratio 
was  adjusted  until  the  maximum  thrust  for  the  given  operating  condition  was  reached. 
This  allowed  for  a  maximum  engine  operating  envelope  of  M  =  2.2,  and  Alt  =  80,000  ft. 
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The  estimated  aircraft  flight  envelope  resulted  in  a  maximum  altitude  of  38,000  ft.  This  is 
less  than  that  of  an  aircraft  with  a  conventional  burner,  but  acceptable  for  flight. 

Frequency  or  fill  fraction  could  have  been  chosen  as  the  independent  throttling 
parameter,  but  it  was  decided  that  frequency  would  remain  fixed  and  that  fill  fraction 
would  be  throttled  for  this  model.  Results  indicate  that  adjusting  either  of  these 
parameters  can  reduce  the  TSFC  at  any  thrust  level  to  roughly  that  of  the  baseline  engine 
and  slightly  lower  at  some  thrust  levels,  particularly  at  cruise  conditions.  These  results 
are  significant  as  incorporating  hybrid  PDEs  into  aircraft  may  save  fuel  costs. 


Recommendations 

The  hybrid  model  described  in  this  research  has  shown  to  yield  perfonnance  gains 
over  a  conventional  engine,  but  not  without  limitations.  The  engine  yields  its  lowest 
TSFC  at  maximum  thrust  which  occurs  when  the  internal  bypass  bleed  is  zero.  The  thrust 
range  is  limited  due  to  the  internal  bypass  of  the  PDC.  The  model  could  also  not  be  run  at 
fill  fractions  greater  than  one  due  to  limitations  of  the  internal  bypass.  Modifying  the 
burner  architecture  to  eliminate  the  internal  bypass  to  accommodate  fill  fractions  greater 
than  one  may  be  considered.  This  modification  would  also  eliminate  the  need  for  an 
internal  mixer,  thus  increasing  the  pressure  of  the  flow  exiting  the  PDC.  Should  the 
internal  bypass  remain  in  the  PDC,  the  static  pressure  entering  MIX39  can  be  controlled 
by  decreasing  the  area  of  the  duct  into  the  mixer.  Such  a  valve  could  allow  for  better 
engine  perfonnance  as  all  flow  would  be  maintained  within  the  PDC.  Throttling 
frequency  as  well  as  fill  fraction  may  result  in  gains.  Additionally,  distributing  all  of  the 
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bleed  flow  from  the  internal  bypass  into  MIX  40  immediately  preceding  the  high  pressure 
turbine  may  result  in  better  perfonnance  than  splitting  the  bleed  flow  equally  between 
MIX40  and  MIX44.  A  boost  pump  should  also  be  implemented  to  increase  the  pressure 
of  the  internal  bleed  flow  into  MIX40  due  to  the  pressure  increase  of  the  flow  exiting  the 
PDC.  Design  parameters  such  as  frequency,  purge  fraction,  equivalence  ratio,  and  tube 
geometry  affect  off-design  performance  and  a  complete  parametric  study  should  be 
perfonned  to  obtain  the  design  parameters  which  yield  the  optimal  design  and  off-design 
perfonnance. 
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Appendix  A.  Derivation  for  Ideal  PDE  Cycle  Thermodynamics 


This  is  a  derivation  of  the  entropy  rise,  thennal  efficiency,  and  Chapman- Jouguet 
Mach  number  according  to  the  Heiser  and  Pratt  thermodynamics  used  in  this  thesis 
(Shapiro,  1953:193). 


Consider  the  flow  through  a  control  volume  of  Fig  A.  1 .  The  continuity  equation 
for  constant  area  is 


k  =  K 

A  (A.l) 


The  momentum  equation  is 


pi  -p2  =  fh  A  (V2-V1) 


and  noting  for  a  perfect  gas  pV2  =  ypM 2 ,  this  may  be  arranged  to  give 


Pi  _1  +  /M2 

Pi  l  +  xAf22 


(A.2) 
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From  the  perfect-gas  law 


Pi  _  Pi  ^2 
Pi  Pi  Ti 


(A. 3) 


Inserting  equations  A.l  and  A.2  into  A. 3  yields 


T2  1 +/M{V2 
Tx  1  +  /M22Vl 


(A.4) 


From  the  definition  of  the  Mach  number  and  perfect  gas 

M 2  _V2c1  _  V2  lj[ 
M\  ~  Vx  c2  ~  V,  \  T2 

Using  the  value  of  V2/V 1  from  Eq.  A. 5,  Eq.  A.4  becomes 

r,  Ml  (i  +  r^i2)' 


From  the  energy  equation 
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Elimination  of  T1/T1  from  Eqs.  A. 6  and  A. 7  yields 


(A. 5) 


(A.  6) 


(A.  7) 


76 


(A. 8) 


t02  _  m;  (i  +  /Mi  y  l+r 2 

T0i  Mf  (l  +  yM;  )2  1  +  LZ 1  Mj2 

Setting  Mi=l  and  M2=MCj  (Shapiro,  1953:195) ,  Eq.  A. 6  becomes 


t,  (i  ZLZk 

1 ‘  (l  +  rM2CJf 


Eq.  A. 8  becomes 

^  2q  +  l)Ma,(A^M;. 
T«>  {i+rM'Z 


Similarly, 


v2  _pi  _{y+x)M2cj 

p2  i  yMCJ 

p2 _  r+ 1 

Pi  i+yMcj 

From  the  definition  of  isentropic  stagnation  pressure 

BlZ+L±m>Y' 
p  \  2  J 


(A. 9) 


(A. 10) 


thus, 


77 


P  02  _  Pl_ 
Pox  Px 


V 


2 


2 


(A.  11) 


f 

1  + 

V 


Y~ 1 
2 


M 


J 


7- 1 


J 


Substituting  Eq.  A.  10  and  Mi=l  and  M2=Mcj  into  Eq.  A.l  1  yields 


P  02 
Pox 


y  +  l 

i +yM2CJ 


1+Y _ ^-M2 

1 T  2  CJ 


y- 1 


y  +  l 


Defining  the  change  in  entropy  as 


^  =  ln 


T2/Tx 


r- 1 

{p2  /  Px)  r 


(A.  12) 


and  substituting  Eq.  A. 9  and  A.  10  into  Eq.  A.  12  yields  the  Heiser  and  Pratt  change  in 
entropy: 


=  -In 


M 1 

1V1CJ 


y+ 1 

f  y  +  l  ^  r 


[i+rKjJ 


Deriving  thermal  efficiency 


(  \ 

(T*-T)  =  c„T' 

exp 

fio  ~  *0 

-1 

c 

\  p  J 

=  CpT0 


exp 


f  \ 

*4  -^3 

v  c,  J 


-1 


=  CpT0 


fl  ,  lfl  \ 


CJ 


r+i 


M. 


CJ 


1  +  yM  t 


y  + 1 


and  since 


9ai3 


and  q  = 

cSn 


the  thennal  efficiency  becomes 


n tk 


=i- 


1 

'l+rMtj) 

r - 

Ik... 

l  /+1  J 

v-l 

r 


To  define  Mcjwe  first  rearrange  the  mass  and  momentum  equations 


Pi  +  P\K  =  P2  + 


PLV> 

Pi 


to  yield 


V?=v\ 


f  \ 

Pi  ~  Pi 

V  V1“V2  J 

f  \ 

Pl~  Pi 

V  V1~V2  J 


The  energy  equation  \  +  ^-l/l2  =  h2  +  V;  may  be  rearranged  with  the  momentum 
equation  to  yield 
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>h-\=\(v?-vi) 

Substituting  the  mass  balance  equation  we  get  the  Hugoniot  relation: 

K~K  =-(/?2-a)(vi +v2> 

Assuming  a  perfect  gas  and  hx  =  cpT ,  h2  =  cpT2  -  qadd  we  get 

(  Pl^l  -  P\V\  )  ~  Qadd  =  T  ( P2  ~  Pi  )(V2  +  V1 ) 
y  —  t  / 

The  expression  for  the  Mach  number  of  the  Rayleigh  process  at  station  1  can  be  written 

as 


(A. 13) 


(A.  14) 


M 


2 

1 


(  \ 
vi  P2-P1 

yp\v i-v2  , 


(A.  15) 


Solving  eqs.  A.  13  and  A.  14  to  detennine  the  volume  and  pressure  ratios  in  terms  of 
approaching  Mach  number  yields 

v2  1 ±  -  if  -2(r + i)(r  -  i)m^  /  «,2 

v,  ~  (r+i)M? 

Pl  i+rPl  +  ^K-i)i-2(r+i)(y-i)^r^^i 

Pi  (r  +  1) 


For  any  Mi  the  value  of  qadd  is  found  by  setting  the  quantity  under  the  radical  equal  to 
zero.  This  yields 


y- 1 

2 


2(y  +  l)AfJ, 
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jSolving  for  Mcj  with  q  =  cpT0  and  ^  =  qT$  f  ^yields 


Ma  =(/+l)(9/v/)+l  +  J[(/+l)(?/»  +  l]‘-l 

V  (A.  16) 
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Appendix  B.  Model  File  for  Baseline  Turbofan  Engine 


// 

// - 

//TURBOJET  ENGINE  BUILD  | 

//  I 

//  BUILD  AND  VERIFY  TURBOJET  | 

//  I 

//DESIGN  POINT  ONLY  | 

//  I 

// - 

//  TURBOJET  CONFIGURATION 

// - 

cout «  "\t - \n" 

«  "\t  Baseline  High  Bypass  Turbofan  built  to  match  AEDsys  \n" 
«  "\t - \n\n"; 

//  Set  model  name 

MODELNAME  =  "Baseline  HBTF  CmpareAEDsys.mdl  with  mixers"; 

// - 

//  set  the  thermo  package 

// - 

setThermoPackage("GasTbl"); 

//  setThermoPackage("Janaf'); 

// - 

//  include  the  standard  intepretted  things 

// - 

#include  <lnterplncludes.ncp> 

#include  "ncp.view" 

//#include  "bleed_macros.fnc" 

//#include  "NewDuct.int" 

// - 

//  #include  the  definition  file  for  the  user  defined  engine 
//  performance  component 

// - 

#include  "EngPerf. cmp"  ; 


// - 

//  MODEL  DEFINITION 
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// - 

//  mmmmmmmmmm  flight  conditions  mmmmmmmmmm# 


Element  FlightConditions  AMBIENT  { 


alt  =  0.0; 

MN  =  0.01; 

//  Ps  =  14.696; 
//  Ts  =  59.0; 
W=  1500.00; 


} 


//  Specify  Design  conditions 
//  design  altitude  (ft) 

//  design  Mach  number 
//  ambient  pressure  (psia) 

//  ambient  temperature  (F) 

//  design  mass  flow  (lbm/s) 


//mmmmmmmmmmmmm#  inlet  mmmmmmmmmmmmmm 

Element  Inlet  INLET  { 


eRamBase  =  0.995;  //Ram  Recovery  Factor? 

} 


//  mmmmmmmmmm##  splitter  ############################### 

Element  Splitter  SPLIT  { 

BPR=  8.0;  // Bypass  Ratio 


} 

// ##########################  FAN  ############################### 

//  here  the  fan  represents  the  outer  portion  of  the  Low  pressure 
//  compressor  spool 
Element  Compressor  Fan21  { 

//  //  use  these  lines  if  no  compressor  map  is  imlemented 
//  effDes  =  0.88042; //0.882886; 

//  PRdes  =  1.56; 

//  use  these  lines  if  compressor  map  is  used... 

#include  "fan.map"  ;  //Compressor  sub-element  map 
Smap. effDes  =  0.8827;  //0. 88289; 

Smap. PRdes  =  1.56; 

} 

//  #####################  Bypass  Duct /  Nozzle/  Sink  ################### 

Element  Duct  Bypass  13  { 

//  AEDsys  assumes  flow  in  bypass  duct  is  isentropic 
//  dPqPbase  =  0.015;//  pressure  loss  through  the  bypass  duct 

} 
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Element  Nozzle  Noz  18  { 

//  Cfg  =  0.995; 

dPqP  =  1.0-0.98;  //  pressure  loss  from  nozzle  inlet  to  throat 
PsExhName  =  "AMBIENT.Fl_O.Ps"; 

//  AEDsys  uses  a  fixed  convergent  nozzle  for  bypass  exit 
switchType  =  "CONIC"; 

} 

Element  Flo  wEnd  Noz  Sink  19  { 

} 

//  mmmmmmmm  Low  Pressure  Compressor  ########################### 

Element  Compressor  LPC20  { 

//  //  use  these  lines  if  no  compressor  map  is  implemented 
//  effDes  =  0.88042;  //  set  the  design  point  isentropic  efficiency 
//  PRdes=  1.56; 

//  use  these  lines  if  compressor  map  is  used... 

#include  "lpc.map"; 

Smap. effDes  =  0.8827;//  set  design  point  isentropic  efficiency 
Smap.PRdes  =  1.56; 


} 

//  ######################  High  Pressure  Compressor  ################## 
Element  Compressor  HPC25  { 

//  //  use  these  lines  if  no  compressor  map  is  implemented 
//  effDes  =  0.85755;  //  set  the  design  point  isentropic  efficiency 
//  PRdes=  16.66667; 

//  use  these  lines  if  compressor  map  is  used... 

#include  "hpc.map"  ;  //  Compressor  sub  element  map 
Smap.  effDes  =  0.8573  ;  //  design  point  isentropic  efficiency 
S_map.PRdes  =  16.66667  ;  //  Set  the  pressure  ratio  at  design 


} 

//  mmmmmmmm#  Bleed  starting  point  ######################### 

Element  Bleed  BLD3  { 

//  =========================  BLEEDS  ==================== 
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//  Three  Bleeds  are  taken  off  of  the  back  side  of  the 

//  High  pressure  Compressor 

BleedOutPort  BL_Cool_301  { 

fracW  =  0.05;  //  mass  flow  (5%  for  cooling  turbine) 

} 

BleedOutPort  BL_Cool_302  { 

fracW  =  0.05;  //  mass  flow  (5%  for  cooling  turbine) 

} 

BleedOutPort  BL_Env_303  { 

fracW  =  0.01;  //  mass  flow  fraction  (1%  bleed) 

} 


} 

//  mmmmmmmmmmmmmm  Fuel  mmmmmmmmmmmmmm# 

Element  FuelStart  FUEL32{ 

LHV  =  18400;  //  BTU/lbm  -  Lower  Heating  Value  of  the  fuel  - 

//  default  is  18400  BTU/lbm 


} 

//  mmmmmmmmmmmmmm##  Burner  ######################### 

Element  Burner  BRN36{ 
effBase  =  0.995;  //  component  efficiency 

dPqPBase  =  1.0  -  0.96;  //pi  b  =  1.0-(dP/P)  pressure  drop  across  burner 

//  Change  from  burner  default  of  FAR  to  TEMPERATURE 
switchBurn  =  TEMPERATURE; 

//  Total  temp,  at  exit  (degrees  Rankine)  ||  not  to  be  used  with  FAR 
TtCombOut  =  2900.0; 


} 

//  ##########################  Bleed  Mixer/IGV  ########################### 

Element  Bleed  MIX40  { 

BleedlnPort  Blln40{ 

Pscale  =  0.88; 

} 

} 

//  ##########################  HP  Turbine  ########################### 
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Element  Turbine  HPT41  { 

#include  "hpt.map";  //High  Pressure  Turbine  Map 
Smap.effDes  =  0.9057;//0.90555;0.91075; 

//  InterStageBleedlnPort  Blln41  { 

//  Pfract  =  1 ;  //force  the  bleed  to  come  in  at  enterance 

//  } 

} 

//  mmmmmmmmmmmmm  Bleed  Mixer  ########################### 

Element  Bleed  MIX44  { 

BleedlnPort  Blln44{ 

Pscale  =  0.68; 

} 

//  F1I1.MN  =  .29; 

//  B1  II. Pscale  =  0.92;  //  Scale  pressure  so  that  the  pressure  ratio  across  mixer  =  1 
//  B1_I1.MN  =  0.31; 

} 


//  mmmmmmmmmmmmm  lp  Turbine  mmmmmmmmmmmmm# 

Element  Turbine  LPT45  { 

#include  "lpt.map"  //Low  Pressure  Turbine  Map 
Smap.effDes  =  0.9084;//0.90836;0.90906; 

//  InterStageBleedlnPort  Blln44  { 

//  Pfract  =1.;  //  force  bleed  to  come  in  at  turbine  entrance 

//  } 

} 

//  #########################  Nozzle  ####################### 

Element  Nozzle  Noz8  { 

//Cfg  =  0.995; 

//Cv  =  0.985; 
dPqP  =  1.0-0.985; 

PsExhName  =  "AMBIENT.Fl_O.Ps"; 

switchType  =  "CONIC";  //  AEDsys  uses  a  fixed  convergent  nozzle  for  core  exit 


} 

//  ##########################  Tenninate  Flow  ################ 
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Element  FlowEnd  Sink39  { 

//  BleedlnPort  Blln44{ 

//  Pscale  =  0.96; 

//  } 

//  sink  for  the  environmental  bleed... 

} 

Element  FlowEnd  NozSink9  { 

//  sink  for  the  core  airflow 

} 


// 

%%0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o%%% 

%%%%%%%%% 

//  Put  shafts  in  the  model 

// 

%%0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o0/o%%% 

%%%%%%%%% 


//mm#####################  Low-Pressure  Shaft  ################ 

Element  Shaft  LPShf  { 

ShaftlnputPort  LPC,  FAN,  LPT  ; 

Nmech  =  2000.0; 

inertia  =  1.0;  //  inertia  is  only  needed  for  transient  analysis 
HPX  =  0.0  ;//  +131.;  //+92.30;  //  Horsepower  extracted  from  the  shaft  hp  (  =  325.7 
kW) 

fracLoss  =  1.0  -  0.99;  //  Fractional  loss  on  positive  port  torque  (1.0  -  etam) 

} 

//mmmmmmmmmmmm#  High  Pressure  shaft  mmmmmmmmm# 

Element  Shaft  HPShf  { 

ShaftlnputPort  HPT,  HPC  ; 

Nmech  =  1 1000.0; 
inertia  =1.0; 

HPX  =  143.178  ;//+372;//  +415.;//  +400.0;  //  Horsepower  extracted  from  the  shaft  hp 
(=  105.7  kW)/  etam  (  =  0.99) 

fracLoss  =  1.0  -  0.99;  //  Fractional  loss  on  positive  port  torque  (1.0  -  eta  m) 

//cout « inertia.unitsunits  «endl; 

//quit(); 

} 

//mmmmmmmmmmmm  Engine  Performance  mmmmmmmmmmm 

Element  EngPerf  PERF{ 
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} 


// _ 

//  Flow  Connections  // 

//  // 

//  This  is  where  the  flow  is  defined  for  the  engine  // 

// _ // 

II 


nmmmmmm#  Ambient  to  splitter  mmmmmmmmmmmm# 

linkPorts(  " AMBIENT. F1_0",  " INLET. F1I",  "FLO"  ); 

linkPorts(  "INLET.FI  O",  "SPLIT.F1  I",  "FL1"  ); 


II#############  Bypass  air 

############################# 

linkPorts(  "SPLIT.F102", 

"Fan21.Fl  I", 

"FLb2"  ); 

linkPorts(  "Fan21.Fl  O", 

"Bypass  13. FI  I", 

"FLb3"  ); 

linkPorts(  "Bypass  13. FI  O", 

"Nozl8.Fl  I", 

"FLb7"  ); 

linkPorts(  "Nozl8.Fl  O", 

"NozSinkl9.Fl_I", 

"FLb8"  ); 

limmmmmm#  Core  Air  Flow  ############################# 


linkPorts(  "SPLIT. FI  01", 

"LPC20.F1  I", 

"FL2"  ); 

linkPorts(  "LPC20.F1  O", 

"HPC25.F1  I", 

"FL25"  ); 

linkPorts(  "HPC25.F1  O", 

"BLD3.F1  I", 

"FL3"  ); 

linkPorts(  "BLD3.F1  O", 

"BRN36.F1  I", 

"FL31"  ) ; 

/ /4i4i4i4i4i4i4i4i4i#44i#44i#4  T71  UUUUUUUUUUUUUUUUUUUU-U-I4-I4-I4-U--I4-I4-I4-I4-I4 

linkPorts(  "FUEL32.Fu  O", 

"BRN36.Fu  I", 

"Fu3" ); 

linkPorts(  "BRN36.F1  O", 

"MIX40.F1  I", 

"FL4"); 

linkPorts(  "MIX40.F1  O", 

"HPT41.F1  I", 

"FL41"  ); 

linkPorts(  "HPT41.F1  0", 

"MIX44.F1  I", 

"FL44"); 

linkPorts(  "MIX44.F1  0", 

"LPT45.F1  I", 

"FL45"  ); 

linkPorts(  "LPT45.F1  0", 

"Noz8.Fl  I", 

"FL7"); 

linkPorts(  "Noz8.Fl  0", 

"NozSink9.Fl  I", 

"FL8"  ); 

limmmmmmm  Bleed  port  linkage  ########################## 

//linkBleedCB("BLD3",  "MIX40",  0.05,  1.0,  1.0,  "BL  1"); 
//linkBleedCB("BLD3",  "MIX44",  0.05,  1.0,  1.0,  "BL  2"); 
//linkBleedCB("BLD3",  "Sink39",  0.01,  1.0,  1.0,  "BL  3"); 
linkPorts(  "BLD3.BL_Cool_301",  "MIX40.BlIn40",  "BL  1"); 

linkPorts(  "BLD3.BL_Cool_302",  "MIX44.BlIn44",  "BL  2"); 

linkPorts(  "BLD3.BL_Env_303",  "Sink39.FlJ",  "BL  3"); 

//d>d>d>d>d>4)d>d>4>4)4)d>4>d>d>d>4>d>d>d>4>d>d>d>4>4)d>d>4>4)4)d>4>d>4)d>4>4)4)d>4>d>d>d>4>d>4)d> 

//  Mechanical  (Shaft)  connections 
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//  LPd>4>d>4>d>4>d>4>d>4>d>4>d>4>4>4>d>4>d>4>d>4>4>4>d>4>4>4>d>4>4>4>d>4>4>4>d>4>d>4>d>4>d>4>d>4>4> 


//###############  Low-Pressure  Spool  ####################### 

linkPorts("LPC20.Sh_O",  "LPShf.LPC",  "LP1"); 
linkPorts("LPT45.Sh_0",  "LPShf.LPT",  "LP2"); 
linkPorts("Fan2 1  .Sh_0",  "LPShf.FAN",  "LP3"); 

//##############  High-Pressure  Spool  ####################### 

linkPorts("HPC25.Sh_0",  "HPShf.HPC","HPl"); 
linkPorts("HPT41.Sh_0",  "HPShf.HPT",  "HP2"); 

j I  A AA AA AA AAA AA AAA AA AA AAA AA AA AA AA AA AA AA  AA AAAA AA AA AA AA AA AA AA 

//  Begin  Run  Definition 

// vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 

COUt  «  "AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^" 

«  "  Begin  Run  Input  definitions  \n  " 

« " vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv\n\n" ; 
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Appendix  C.  Model  File  for  Hybrid  Turbofan  Engine 


// 

// - 

//  HYBRID  TURBOFAN  ENGINE  | 

//  I 

// - 

//  CONFIGURATION 

// - 

cout «  "\t - \n" 

«  "\t  Hybrid  Pulsed  Detonation  Combustor  High  Bypass  Turbofan  ...  \n" 

«  "\t - \n\n"; 

//  Set  model  name 

MODELNAME  =  "PDC  HBTF";  //Pulsed  Detonation  Combustor  High  Bypass 
Turbofan"; 

// - 

//  set  the  thenno  package 

// - 

setThennoPackage("GasTbl"); 

//  setThermoPackage("FPT"); 

// - 

//  include  the  standard  intepretted  things 

// - 

#include  <lnterplncludes.ncp> 

#include  "ncp.view" 

// - 

//  #include  the  definition  file  for  the  user  defined  engine 
//  performance  component 

// - 

#include  "EngPerf.cmp"  ; 


// - 

//  MODEL  DEFINITION 

// - 

//  ####################  FLIGHT  CONDITIONS  ##################### 

Element  FlightConditions  AMBIENT  { 

//  Specify  Design  conditions 
alt  =  0.0;  //  design  altitude  (ft) 

MN  =  0.01;  //  design  Mach  number 
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//  Ps  =  14.696;  //  ambient  pressure  (psia) 

//  Ts  =  59.0;  //  ambient  temperature  (F) 

W  =  1500.00;  //  design  mass  flow  (lbm/s) 

} 

//###########################  inlet  mmmmmmmmmmmmmm 

Element  Inlet  INLET  { 

eRamBase  =  .995;  //Ram  Recovery  Factor?  //.995 

} 

//  ######################  Splitter  ############################### 

Element  Splitter  SPLIT  { 

BPR=  8.0;  // Bypass  Ratio 


} 

//  mmmmmmmmmmmmm  fan  mmmmmmmmmmmmmm### 

//  here  the  fan  represents  the  outer  portion  of  the  Low  pressure 
//  compressor  spool 
Element  Compressor  Fan21  { 

//  //  use  these  lines  if  no  compressor  map  is  imlemented 
//  effDes  =  0.88042; //0.882886; 

//  PRdes  =  1.56; 

//  use  these  lines  if  compressor  map  is  used... 

#include  "fan.map"  ;  //Compressor  sub-element  map 
Smap.  effDes  =  0.8827;  //0.88289;//.8827 
Smap.  PRdes  =  1.56; 

} 


//  #####################  Bypass  Duct/  Nozzle/  Sink  ################### 

Element  Duct  Bypass  13  { 

//  AEDsys  assumes  flow  in  bypass  duct  is  isentropic  (pi 09,  #9) 

//  dPqPbase  =  0.015;  //  pressure  loss  through  the  bypass  duct 

} 


Element  Nozzle  Noz  18  { 

//Cfg  =  0.995; 

dPqP  =  1.0-0.98;  //  pressure  loss  from  nozzle  inlet  to  throat 

PsExhName  =  MAMBIENT.Fl_O.Ps"; 
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switchType  =  "CONIC";  //  AEDsys  uses  a  fixed  convergent  nozzle  for  bypass  exit 

) 

Element  Flo wEnd  Noz Sink  19  { 

} 

//  mmmmmmmm  Low  Pressure  Compressor  ########################### 

Element  Compressor  LPC20  { 

//  //  use  these  lines  if  no  compressor  map  is  imlemented 
//  effDes  =  0.88042;  //  set  the  design  point  isentropic  efficiency 
//  PRdes=1.56; 

//  use  these  lines  if  compressor  map  is  used... 

#include  "lpc.map"; 

Smap. effDes  =  0.8827;  //  0.88289;  set  the  design  point  isentropic 

//efficiency//.  8  827 
Smap.PRdes  =  1.56; 

} 

//  mmmmmmmmmmm  High  Pressure  compressor  mmmmmmmmm 

Element  Compressor  HPC25  { 

//  //  use  these  lines  if  no  compressor  map  is  imlemented 

//  effDes  =  0.85755;  //0. 8855338;  //  set  the  design  point  isentropic  efficiency 

//  PRdes=  16.66667; 

//  use  these  lines  if  compressor  map  is  used... 

#include  "hpc.map"  ;  //  Compressor  sub  element  map 
S  map. effDes  =  0.8573  ;  //0. 857535  ;  set  the  maps  design  point 
//isentropic  efficiency//.  8  5  73 

S_map.PRdes  =  16.66667  ;  //  Set  the  pressure  ratio  at  design 

} 

//  mmmmmmmmmmm  Bleed  starting  point  mmmmmmmmmmm### 

Element  Bleed  BLD3  { 

//  ============================  BLEEDS 


//  Three  Bleeds  are  taken  off  of  the  back  side  of  the  High  pressure  Compressor 
BleedOutPort  BL_Cool_301  { 
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//fracBldWork  =  1.0;  //  work  fraction  where  bleed  is  taken 
//fracBldP  =  1.0;  //  Pressure  fraction  where  bleed  is  taken 

fracW  =  0.05;  //  mass  flow  (5%  for  cooling  turbine) 

} 

BleedOutPort  BL_Cool_302  { 

//fracBldWork  =  1.0;  //  work  fraction  (dhb/dh) 

//fracBldP  =  1.0;  //  Pressure  fraction  (dPb/dP) 

fracW  =  0.05;  //  mass  flow  (5%  for  cooling  turbine) 

} 

BleedOutPort  BL_Env_303  { 

//fracBldWork  =  1.0;  //  (dhb/dh)  work  fraction  -  closely  tied  with 
pressure  fraction... 

//fracBldP  =  1 .0;  //  Pressure  Fraction  (dPb/dP) 

fracW  =  0.01;  //  mass  flow  fraction  (1%  bleed) 

} 

} 

//  mmmmmmmmmmmmmm  Fuel  mmmmmmmmmmmmmm# 

Element  FuelStart  FUEL32{ 

LHV  =  18400;  //  BTU/lbm  -  Lower  Heating  Value  of  the  fuel  - 

//  default  is  18400  BTU/lbm 

} 

//  mmmmmmmmmmmmmm##  Burner  ######################### 

//  Element  Burner  BRN36{ 

//  effBase  =  0.995;  //  component  efficiency  ...  ?? 

//  dPqPBase  =  1.0  -  0.96;  //0.04;  //  pi  b  -  pressure  drop  across  burner...  ??  (dP/P) 

// 

//  switchBum  =  TEMPERATURE;  //  Change  from  burner  defauls  using  Fuel-air  Ratio 
(FAR)  to  TEMPERATURE 

//  TtCombOut  =  2900.0;  //  Total  temperature  at  exit  (degrees  Rankine)  ||  not  to  be 

used  with  FAR 

// 

//  //  or  use  the  default  FAR  and  define  what  the  FAR  is... 

//  //  FAR  =  0.02282;  //  Fuel-to-Air  ratio;  not  to  be  used  with  TtCombOut, 

Wfuel,  etc. 

// 

//} 

#include  "PDC  burner  bleed. inf ' 
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Element  PulseDetonationCombustor  BRN36{ 
effBase  =  .995;  //  burning  efficiency//.995 

dPqPBase  =  1.0-0.96;  //  pressure  loss  across  valves/through  bypass 

switchBum  =  FAR;  //  set  fuel-air  ratio  (vs  equivalence  ratio) 

FAR  =  (0.0683  *  1.00);  //approximately  85%  of  stoichiometric  conditions 

purgeFrac  =  0.2;  //  designate  purge  fraction 

iillFrac  =  0.8;  //  designate  fill  fraction 

lTube  =  36;  //  length  of  tube  in  inches 

n_tubes  =  24;  //  number  of  tubes 

dTube  =  2.0;  //  inside  diameter  of  tubes 

tCycle  =  .016776271641;  //  cycle  time 

flowby  =  1 ;  //  percentage  of  internal  bypass  flow  into  mixer39 

} 

//  mmmmmmmmmmmmm  wail  heat  exchange  mmmmmmmmmmmm 

//  ***  not  uses  in  the  current  model  *** 

//Element  Wall  WAFF38{ 

//  Ahxl  =  PI*36;  //  area  of  wall  inside  PDT 

//  Ahx2  =  PI*36*  1 .02;  //  area  that  bypass  flow  sees 

//  ChxDesl  =  0.7;//  heat  transfer  film  coefficient  -  blind  guess... 

//  ChxDes2  =  0.7;// 

//  CpMat  =  0.1481;//specific  heat  of  material  (titanium  @  2160  R) 

//  //  #  tubes  pi/4  length  oD  iD(in)  rho(lbm/ftA3)  Titanium 

//  massMat  =  36.*(PI/4.*(36./12.)*(2.25**2-2.**2)/144.)*280.93;//mass  of  material  in 

lbm 

// 

//} 

//  ######################  Internal  Bypass  Bleed  ######################### 

Element  Bleed  BFD4  { 

//  ============================  BREEDS 


//  Three  Bleeds  are  taken  from  the  internal  bypass  of  the  PDC 

BleedOutPort  BL_Cool_304  { 

//fracBldWork  =  1.0;  //  work  fraction  where  bleed  is  taken 
//fracBldP  =  1 .0;  //  Pressure  fraction  where  bleed  is  taken 

fracW  =  0.499;  //  mass  flow  (50%  for  cooling  turbine) 

} 

BleedOutPort  BL_Cool_305  { 

//fracBldWork  =  1.0;  //  work  fraction  (dhb/dh) 

//fracBldP  =1.0;  //  Pressure  fraction  (dPb/dP) 

fracW  =  0.499;  //  mass  flow  (50%  for  cooling  turbine) 

} 
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} 


//  #############  PDC  bypass  mixer/Transition  to  steady-state  device  ##### 

Element  Mixer  MIX39{ 

Fi  ll .MN  =  .95;  //  Rather  high  MN,  but  it  works  where  lower 
//  values  do  not... 

} 

//  mmmmmmmmmmmmm  Bleed  Mixer/iGv  ########################### 

Element  Bleed  MIX40  { 

BleedlnPort  Blln40{ 

Pscale  =  0.88; 

} 

BleedlnPort  Blln41  { 

Pscale  =  .88; 

} 


//  ##########################  HP  Turbine  ########################### 

Element  Turbine  HPT41  { 

#include  "hpt.map";  //High  Pressure  Turbine  Map 
Smap.effDes  =  0.9057;//0.90555;0.91075;//.9057 


} 

//  mmmmmmmmmmmmm  Bleed  Mixer  mmmmmmmmmmmmm# 

Element  Bleed  MIX44  { 

BleedlnPort  Blln44{ 

Pscale  =  0.68; 

} 

BleedlnPort  Blln45  { 

Pscale  =  .68; 

} 


// ##########################  LP  Turbine  ########################### 

Element  Turbine  LPT45  { 

#include  "lpt.map"  //Low  Pressure  Turbine  Map 
Smap.effDes  =  0.9084;//0.90836;0.90906;//.9084 
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} 


//  #########################  Nozzle  ####################### 

Element  Nozzle  Noz8  { 

//Cfg  =  0.995; 

//Cv  =  0.985; 

dPqP  =  1.0-0.985;  //  pressure  loss  across  the  nozzle 
PsExhName  =  "AMBIENT.Fl_O.Ps"; 

switchType  =  "CONIC";  //  AEDsys  uses  a  fixed  convergent  nozzle  for  core  exit 


} 

//  ##########################  Terminate  Flow  ################ 

Element  FlowEnd  Sink39  { 

//  sink  for  the  environmental  bleed... 

} 

Element  FlowEnd  NozSink9  { 

//  sink  for  the  core  airflow 

} 


Element  FlowEnd  NozSinkl  { 

//  sink  for  the  ibypass  bleed  airflow 

} 

// 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%%%%%%%% 

//  Put  shafts  in  the  model 

// 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%%%%%%%%% 


//#########################  Low-Pressure  Shaft  ################ 

Element  Shaft  LPShf  { 

ShaftlnputPort  LPC,  FAN,  LPT  ; 

Nmech  =  2000.0; 

inertia  =  1.0;  //  inertia  is  only  needed  for  transient  analysis 

HPX  =  0.0  ;  //+92.30;  //  Horsepower  extracted  from  the  shaft  hp  (  =  325.7  kW) 

fracLoss  =  1.0-.99;  //  Fractional  loss  on  positive  port  torque  (1.0 

//-  eta_m)1.0-.99 

} 
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//mmmmmmmmmmmm#  High  Pressure  shaft  mmmmmmmmm# 

Element  Shaft  HPShf  { 

ShaftlnputPort  HPT,  HPC  ; 

Nmech  =  1 1000.0; 

inertia  =  1.0;//  inertia  is  only  needed  for  transient  analysis 

HPX  =  143.178  ;//  +400.0;  //  Horsepower  extracted  from  the  shaft  hp  (  =  105.7  kW)/ 
eta  m  (  =  0.99) 

fracLoss  =  1.0  -  .99;  //  Fractional  loss  on  positive  port  torque  (1.0  -  eta  m) 

//1.0-.99 

} 

//mmmmmmmmmmmm  Engine  Performance  mmmmmmmmmmm 

Element  EngPerf  PERF{ 

} 

// 

// 

// 

// 

// 

// 

//mmmmmm#  Ambient  to  splitter  mmmmmmmmmmmm# 

linkPorts(  "AMBIENT.FI  O",  "  INLET. F1I",  "FLO"  ); 

linkPorts(  "INLET. F1_0",  "SPLIT.F1  I",  "FL1"  ); 

//#############  Bypass  air  ############################# 

linkPorts(  "SPLIT. F102",  "Fan21.Fl  I",  "FLb2"  ); 

linkPorts(  "Fan21.Fl_0",  "Bypass  13. Fl  I",  "FLb3"  ); 

linkPorts(  "Bypass  13. F1_0",  "Nozl8.Fl_I",  "FLb7"  ); 

linkPorts(  "Nozl  8.F1Q",  "NozSinkl9.FlJ",  "FLb8"  ); 

//#############  Core  Air  Flow  ############################# 


linkPorts(  "SPLIT.F101",  "LPC20.F1I",  "FL2" ); 

linkPorts(  "LPC20.F1O",  "HPC25.F1I",  "FL25" ); 

linkPorts(  "HPC25.F10",  "BLD3.F1I",  "FL3" ); 

linkPorts(  "BLD3.F10",  "BRN36.F1I",  "FL31"); 

//##############  Fuel  Flow  ############################## 
linkPorts(  "FUEL32.Fu_0",  "BRN36.Fu_I",  "Fu3"  ); 

//linkPorts(  "BRN36.F1_01",  "WALL38.F1_I1",  "Wal"  ); 

//linkPorts(  "BRN36.F1_02",  "WALL38.F1I2",  "Wa2"  ); 

//linkPorts(  "WALL38.F1_01",  "MIX39.F1_I1",  "F139"); 

//linkPorts(  "WALL38.F1_02",  "MIX39.F1_I2",  "F1392"); 


Flow  Connections  // 

// 

This  is  where  the  flow  is  defined  for  the  engine  // 

// 
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linkPorts(  "BRN36.F101", 
linkPorts(  "BRN36.F1_02M, 
linkPorts(  "BRN36.F103", 
linkPorts(  "BLD4.F10", 
linkPorts(  "MIX39.F10", 
linkPorts(  "MIX40.F1O", 
linkPorts(  "HPT41.F10", 
linkPorts(  "MIX44.F10", 
linkPorts(  "LPT45.F1  O", 
linkPorts(  "Noz8.Fl_0", 


"MIX39.F1I1", 

"MIX39.F1_I2", 

"BLD4.F1I", 

"NozSinkl.FlJ", 

"MIX40.F1I", 

"HPT41.F1I", 

"MIX44.F1J”, 

"LPT45.F1I", 

"Noz8.Fl_I", 

"NozSink9.Fl_I", 


"F139"); 

"F1392"); 

"F1393"); 

"F1394"); 

"FL4"); 
"FL4F'  ); 
"FL44"); 
"FL45"  ); 

"FL7"); 

"FL8" ); 


//##############  Bleed  port  linkage  ########################## 


linkPorts(  "BLD3.BL_Cool_301M, 
linkPorts(  "BLD3.BL_Cool_302", 
linkPorts(  "BLD3.BL_Env_303", 
linkPorts(  "BLD4.BL_Cool_304M, 
linkPorts(  "BLD4.BL_Cool_305", 


"MIX40.BlIn40",  "BL  1"); 

"MIX44.BlIn44",  "BL  2"); 

"Sink39.Fl  I",  "BL  3"); 
"MIX40.BlIn41",  "BL  4"); 

"MIX44.BlIn45",  "BL  5"); 


//4>4>d>4>4>d>4>d>4>4>4>d>4>4)4>d>4>4)4>4)4>d>4>d>4>4)4>4)4)d>4)d>4>4)4>d>4>d>4>d>4>d>4>d>d>d>4)4) 

//  Mechanical  (Shaft)  connections 

//  4>4>d>d>4>d>4)d>4)d>4)d>4)4)4)d>4)4)4)d>4)4)d>d>4)4)4)d>4)d>d>d>4)4)4)d>4)d>4)d>4)d>4)d>4)d>d>d> 


//###############  Low-Pressure  Spool  ####################### 

linkPorts("LPC20.Sh_O",  "LPShf.LPC",  "LP1"); 
linkPorts("LPT45.Sh_0",  "LPShf.LPT",  "LP2"); 
linkPorts("Fan2 1  .Sh_0",  "LPShf.FAN",  "LP3"); 

//##############  High-Pressure  Spool  ####################### 

linkPorts("HPC25.Sh_0",  "HPShf.HPC","HPl"); 
linkPorts("HPT41.Sh_0",  "HPShf.HPT",  "HP2"); 

j I  A AA  AA AA AA AA  AA AAA AA  AAA AA AA AA AA AA AA  AA AA AAAAA AA AA  AA AA AA AA AA 

//  Begin  Run  Definition 

// vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 


COllt  «  ’’AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^II 

«  "  Begin  Run  Input  definitions  \n  " 

« " vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv\n\n" ; 


98 


Appendix  D.  Pulsed  Detonation  Combustor  Code 


#ifndef _ PDC _ 

#define  _PDC_ 

//  *  Air  Force  Institute  of  Technology 
//  *  2950  Hobson  Way,  Bldg  641 
//  *  Wright  Patterson  AFB,  OH  45433 

//* 

//  *  Written  by  Ionio  Q.  Andrus,  Capt.,  USAF 
//  *  Modified  by  Caitlin  R.  Thorn,  Capt.,  USAF 

//  BASED  ON  "Burner.int"  included  in  NPSS,  written  by — 

//  *  NASA  Glenn  Research  Center 
//  *  21000  Brookpark  Rd 
//  *  Cleveland,  OH  44135 

//* 

#include  <lnterplncludes.ncp> 

class  PulseDetonationCombustor  extends  Element  { 

// - 

n  *******  DOCUMENTATION  ******* 

// - 

title  = 

description  =  isA()  +  "  will  calculate  perfonnance  for 
pulsed  detonation  combustor."; 

usageNotes  =  " 

The  burner  element  performs  high  level  burner  perfonnance 
calculations.  This  element  works  with  an  entrance  fluid  and 
fuel  stream.  It  mixes  the  two  flows  together  and  then 
performs  the  burn  calculations.  Please  note  that  the  burner 
has  no  control  over  the  actual  fuel  stream  conditions— fuel  type, 
LHV,  etc.  These  values  are  properties  of  the  fuel  flow  itself 
and  are  usually  set  in  the  FuelStart  element. 
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There  are  two  ways  to  specify  the  burner  exit  conditions.  The 
first  way  is  specify  the  burner  fuel-to-air  ratio.  The  second 
way  is  to  set  equivalence  ratio.  The  type  of  input  used  is 
controlled  by  an  option  switch. 

The  burner  tracks  several  different  pressure  losses.  The  first, 
dPqP,  accounts  for  duct  friction  pressure  drops  and  approximates 
the  pressure  loss  through  valves.  The  second,  dPqPRayleigh, 
accounts  for  the  Rayleigh  pressure  drop.  dPRayleigh  is  input  or 
calculated  -  see  switchHotLoss,  an  iteration  is  necessary  since 
the  pressure  loss  itself  is  a  function  of  the  exit  conditions. 

The  burner  also  allow  two  efficiencies  to  be  input.  The  first 
efficiency,  eff,  refers  to  the  efficiency  based  on  enthalpy 
change.  The  second  efficiency,  effChem,  refers  to  the  efficiency 
based  on  temperature  change.  Both  terms  can  be  input.  However, 
the  enthalpy  efficiency  is  always  applied  first. 

Additionally, 

The  user  can  request  a  pre  burner  pressure  loss  dPqP.  The 
pressure  loss  calculations  are  performed  before  all  the  other 
calculations  are  done.  This  means  that  the  combustion  entrance 
pressure  will  not  match  the  value  indicated  by  the  burner  entrance. 

The  user  can  request  a  heat  transfer  Qhx.  The  heat  transfer 
calculations  are  performed  after  all  the  other  calculations  are 
done.  This  means  that  if  heat  transfer  is  being  used,  the  exit 
temperature  will  not  match  the  value  indicated  by  the  burner 
calculations. 


background  = 


// - 

H  *******  setup  VARIABLES  ******** 

// - 

real  a  dPqP  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Duct  friction  pressure  drop  adder"; 

} 
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real  a  dPqPAud  { 

value  =  0.0;  IOstatus  =  "unset";  units  =  "psia"; 
description  =  "Audit  factor  adder  applied  to  pressure  ratio"; 

} 

real  a_eff  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Adiabatic  efficiency  adder"; 

} 


real  aeffChem  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Chemical  efficiency  adder"; 

} 

real  ARvalve  {  //  Added  15Feb2007  -  IA 

value  =  0.5;  IOstatus  =  "input";  units  =  "none"; 

description  =  "Ratio  of  valve  throat  area  to  tube  cross  section  area"; 

} 

real  deltaS  {  //Added  17Jan2007  -  IA 

value  =  0.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Change  in  entropy  due  to  detonation"; 


real  DDT  {  //Added  17Jan2007  -  IA 

value  =  0.0005;  IOstatus  =  "input";  units  =  "none";  //seconds 
description  =  "Detonation  to  deflaration  time  in  seconds"; 

} 

real  dPqP  { 

value  =  0.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Adjusted  duct  friction  pressure  drop"; 

} 

real  dPqPBase  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Duct  friction  pressure  drop  "; 

} 

real  dPqPRayleigh  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Adjusted  Rayleigh  pressure  drop"; 


} 


real  dTube  {  //Added  17Jan2007  -  IA 

value  =  2.0;  IOstatus  =  "input";  units  =  "none";  //  inches... 
description  =  "Inside  diameter  of  the  detonation  tube"; 


} 

real  eff  { 

value  =  1.0;  IOstatus  =  "output";  units  =  "none"; 
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description  =  "Adjusted  adiabatic  burner  efficiency"; 

} 

real  effBase  { 

value  =1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Adiabatic  burner  efficiency,  from  socket  "; 

} 

real  effChem  { 

value  =1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Adjusted  chemical  efficiency"; 

} 

real  effChemBase  { 

value  =1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Chemical  efficiency,  from  socket"; 

} 

real  eqRatio  { 

value  =1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Equivalence  ratio  for  fuel-air  mixture"; 

} 

real  FAR  { 

value  =  0.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Fuel-to-air  ratio"; 

} 

real  FARDes  { 

value  =  0.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Fuel-to-air  ratio  at  design"; 

} 

real  fillFrac  {  //Added  17Jan2007  -  IA 

value  =  1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Fill  fraction  "; 

} 

real  flowby  {//added  Dec09  -  CT 

value  =  1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Percentage  of  internal  bypass  into  Mixer39"; 

} 

real  fuelFractV  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Fraction  of  the  incoming  flow  velocity  fuel 
enters  the  burner"; 

} 

real  iBPR  {  //added  17Jan2007  -  IA 

value  =  1.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Bypass  ratio  internal  to  the  PDC"; 
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real  iBPRdes  {  //added  lFeb2007  -  IA 

value  =  1.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Bypass  ratio  internal  to  the  PDC  at 
design  conditions"; 

} 

real  lTube  {  //added  17Jan2007  -  IA 

value  =  36;  IOstatus  =  "input";  units  =  "none";  //inches?? 
description  =  "length  of  the  individual  detonation  tubes"; 

} 

real  n_tubes{  //added  17Jan2007  -  IA 

value  =  36;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Total  number  of  detonation  tubes  used 
in  the  PDC"; 

} 

real  MCJ  {  //added  17Jan2007  -  IA 

value  =  3.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Chapman-Jouguet  Mach  number  of  the 
detonation  wave."; 

} 

real  Mvalve  {  //added  15Feb2007  -  IA 

value  =  1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Mach  number  of  flow  passing  through 
the  valve  throat."; 

} 

real  qadd{  //added  17Jan  2007-  IA 

value  =  0.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Heat  addition  due  to  fuel  combustion"; 

} 

real  Qhx  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "Btu/sec"; 
description  =  "Heat  loss  to  thermal  mass  storage"; 

} 

real  PqPRayleigh  { 

value  =1.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Adjusted  Rayleigh  pressure  drop"; 

} 

real  PqPRayleighDelta  { 

value  =  0.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Bounded  Rayleigh  pressure  drop  -  for  loop  only"; 

} 

real  PqPRayleighError  { 

value  =1.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Adjusted  Rayleigh  pressure  drop  error"; 
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} 

real  PqPRayleighMin  { 

value  =  0.05;  IOstatus  =  "input";  units  =  "none"; 

description  =  "Rayleigh  pressure  drop  lower  limit  -  for  loop  only"; 

} 

real  PqPRayleighStep  { 

value  =  0.05;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Maximum  step  for  Rayleigh  pressure  drop 

-  for  loop  only"; 

} 

real  PqPRayleighNew  { 

value  =1.0;  IOstatus  =  "output";  units  =  "none"; 
description  =  "Previous  adjusted  Rayleigh  pressure  drop 

-  for  loop  only"; 


} 

real  purgeFrac  {  //Added  17Jan2007  -  IA 

value  =  0.25;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Purge  fraction  coefficient  for  flow"; 

} 

real  s  dPqP  { 

value  =  1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Duct  friction  pressure  drop  scalar"; 


} 

real  s  dPqPAud  { 

value  =  1.0;  IOstatus  =  "unset";  units  =  "none"; 
description  =  "Audit  factor  scalar  applied  to  pressure  ratio"; 


} 

real  s_eff  { 

value  =  1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Adiabatic  efficiency  scalar"; 


} 

real  s_effChem  { 

value  =  1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Chemical  efficiency  scalar"; 


} 


real  tauBlDn  {  //  Added  17Jan2007  -  IA 

value  =  5.;  IOstatus  =  "input";  units="none"; 
description  =  "Blowdown  time  constant"; 

} 

real  tauValveOpen  {  //  Added  18Jan2007  -  LA. 

value  =  0.33333;  IOstatus  =  "output";  units="none"; 
description  =  "time  valve  open/  time  cycle  -  from  0  to  1"; 

} 
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real  tCycle{  //  Added  17Jan2007  -  IA 

value  =  0.01;  IOstatus  =  "output";  units  =  "none";  //seconds 
description  =  "Detonation  engine  cycle  time  (=  1/frequency)"; 

} 

real  tolRayleigh  { 

value  =  4e-05;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Iteration  tolerance  on  momentum  pressure  drop"; 

} 

real  tolWfuel  { 

value  =  le-05;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Iteration  tolerance  on  temperature  burn"; 

} 

real  TtCombOut  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "R"; 
description  =  "Exit  temperature"; 

} 

real  TtLast  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "R"; 
description  =  "Previous  exit  temperature  -  for  loop  only"; 

} 

real  TTSSeffj  //  Added  17Jan2007  -  IA 

value  =  1.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Efficiency  factor  for  the  transition  device."; 

} 

real  TTSSdPqP{  //  Added  17Jan2007  -  IA 

value  =  0.0;  IOstatus  =  "input";  units  =  "none"; 
description  =  "Change  in  Pressure  divided  by  Pressure 
for  transistion  to  steady  state  calculation."; 

} 

real  tValve{  //  Added  17Jan2007  -  IA 

value  =  0.0002;  IOstatus  =  "input";  units  =  "none";  //seconds 
description  =  "Time  for  valves  to  open/close"; 

} 

real  Wfuel  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "lbm/sec"; 
description  =  "Combustor  fuel  flow"; 

} 

real  WfuelError  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "lbm/sec"; 
description  =  "Combustor  fuel  flow  error"; 


real  WfuelLast  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "lbm/sec"; 


105 


description  =  "Previous  combustor  fuel  flow  -  for  loop  only"; 

} 

real  WfuelNew  { 

value  =  0.0;  IOstatus  =  "input";  units  =  "lbm/sec"; 
description  =  "Next  combustor  fuel  flow  -  for  loop  only"; 

} 

int  countFuel  { 

value  =  0;  IOstatus  =  "output"; 
description  =  "Fuel  loop  counter"; 

} 

int  countFuelMax  { 
value  =  50;  IOstatus  =  "input"; 
description  =  "Fuel  loop  maximum  counter"; 

} 

int  countRayleigh  { 
value  =  0;  IOstatus  =  "output"; 
description  =  "Rayleigh  loop  counter"; 

} 

int  countRayleighMax  { 
value  =  25;  IOstatus  =  "input"; 
description  =  "Rayleigh  loop  maximum  counter"; 

} 

int  HagRayleighLossTooMuch  { 
value  =  0;  IOstatus  =  "output"; 

description  =  "If  true,  Rayleigh  loop  results  in  too  much  loss"; 

} 

int  flagRayleighChoked  { 
value  =  0;  IOstatus  =  "output"; 

description  =  "If  true,  Rayleigh  loop  results  in  supersonic  flow"; 

} 

//  for  backward  compatibly  with  old  "aud" 

FunctVariable  a  dPqPaud  { 
units  =  "none";  IOstatus  =  "input"; 
getFunction  =  "get_aAud";  setFunction  =  "set_aAud"; 

} 

real  get_aAud()  {  return  a_dPqPAud;  } 

void  set_aAud(real  userValue)  {  a_dPqPAud  =  userValue;  } 

FunctVariable  s  dPqPaud  { 
units  =  "none";  IOstatus  =  "input"; 
getFunction  =  "get_sAud";  setFunction  =  "set_sAud"; 

} 
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real  get_sAud()  {  return  s_dPqPAud;  } 

void  set_sAud(real  userValue)  {  s_dPqPAud  =  userValue;  } 

// - 

//  *******  OPTION  VARIABLE  SETUP  ******* 

// - 

Option  switchAud  { 

allowedValues  =  {  "BASE",  "AUDIT"  } 
description  =  "Determines  if  the  audit  factors  are  used"; 
IOstatus  =  "input"; 
trigger=TRUE; 

} 


Option  switchBum  { 

allowedValues  =  {  "FAR",  "EQRATIO"  };//"FUEL",  "WFUEL", 

"TEMPERATURE",  _  mod  18  Dec  2006  -  IA  -  added  "FILLFRACTION" 

description  =  "Switch  detennines  if  burner  is  running  to  fuel  flow,  FAR,  or  T4.  Setting 
option  to  FUEL  will  burn  using  the  burner  value  as  an  input.  Setting  the  option  to 
WFUEL  will  bum  using  the  value  coming  in  from  the  fuel  station."; 
trigger=TRUE; 

} 

Option  switchDes  { 

allowedValues  =  {  "DESIGN",  "OFFDESIGN"  }; 
description  =  "Design  switch"; 
trigger=TRUE; 

} 

//  input  kept  in  for  backward  compatible  (remove  later) 

Option  switchHotLoss  { 

allowedValues  =  {  "INPUT",  "CALCULATE", "input"  }; 

description  =  "Switch  detennines  if  the  hot  pressure  loss  is  input  or  iterated  on"; 

trigger=TRUE; 

} 

// - 

//  ******  SETUP  PORTS,  FLOW  STATIONS,  SOCKETS,  TABLES  ******** 

// - 

//  FLUID  PORTS 

FluidlnputPort  FI _ I  { 

description  =  "Incoming  flow"; 
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} 

FluidOutputPort  F1_0 1  { 
description  =  "Exiting  combustion  flow"; 

} 

FluidOutputPort  F102  { 
description  =  "Exiting  bypass  flow"; 

} 

FluidOutputPort  F103  { 

description  =  "Exiting  bypass  excess  flow"; 

} 

//  FUEL  PORTS 

FuellnputPort  Fu_I  { 
description  =  "Incoming  fuel  flow"; 

} 

//  BLEED  PORTS 
//  THERMAL  PORTS 
//  MECHANICAL  PORTS 
//  FLOW  STATIONS 

// _ flow  stations  modified  18  Dec  2006-  IA 

FlowStation  Fllcomb  { 

description  =  "Inlet  station  to  detonation  tube  section 
of  burner  (after  the  initial  pressure  loss  is  applied)"; 

} 

FlowStation  FlIcombAir  { 

description  =  "Copy  of  the  inlet  station  to  detonation  tube 
section  of  bumer(after  the  initial  pressure  loss  is  applied, 
before  flow  is  split  and  partitioned)"; 

} 

FlowStation  Fl  lprg  { 

description  =  "Station  containing  detonation  tube  purge  fluid"; 

} 
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FlowStation  Fl  Ocomb  { 

description  =  "Exit  station  to  combustion  section  of  burner 
(before  thennal  storage  heat  transfer  is  calculated)"; 

} 

FlowStation  FI  Vit  { 

description  =  "Vitiated  Fluid  flow  station  before  detonation  (cold)"; 

} 

// _ — end  flow  station  modifications 

//  SOCKETS 

Socket  S  dPqP  { 
allowedValues  =  {  "dPqPBase" }; 

description  =  "Dry  duct  and  valve  pressure  loss";  // _ mod  -IA-  18  Dec  2006 

socketType  =  "dPqP"; 

} 

Socket  S_eff  { 

allowedValues  =  {  "effBase",  "effChemBase"  }; 

description  =  "PulseDetonationCombustor  adiabatic  efficiency"; 

socketType  =  "BURN_EFFICIENCY"; 

} 

Socket  SQhx  { 
allowedValues  =  {  "Qhx"  }; 
description  =  "Thennal  storage  socket"; 
socketType  =  "HEATTRANSFER"; 

} 


//  TABLES 

// - 

//  *******  INTERNAL  solver  setup  ******* 

// - 

// - 

//  ******  ADD  SOLVER  INDEPENDENTS  &  DEPENDENTS  ****** 

// - 

// - 

//  *******  varjable  CHANGED  METHODOLOGY  ******* 
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// - 

void  variableChanged(  string  name,  any  oldVal  )  { 

//  Check  to  see  what  variables  were  changed.... 

//  Change  input/output  status  as  necessary  -  IA-  18  Dec  06 

if(  name  ==  "switchBurn"  )  { 
if  (  switchBurn  ==  "FAR"  )  { 

FAR.IOstatus  =  "input"; 

Wfuel.IOstatus  =  "output"; 

TtCombOut.IOstatus  =  "output"; 
eqRatio.  10  status  =  "output"; 

} 

//  else  if  (  switchBurn  ==  "FUEL"  )  { 

//  FAR.IOstatus  =  "output"; 

//  Wfuel.IOstatus  =  "input"; 

//  TtCombOut.IOstatus  =  "output"; 

//  } 

//  else  if  (  switchBurn  ==  "WFUEL"  )  { 

//  FAR.IOstatus  =  "output"; 

//  Wfuel.IOstatus  =  "output"; 

//  TtCombOut.IOstatus  =  "output"; 

//  } 

// _ added  5  Feb  2007  -IA- 

else  if  (  switchBurn  ==  "EQRATIO"  )  { 

FAR.IOstatus  =  "output"; 

Wfuel.IOstatus  =  "output"; 

TtCombOut.IOstatus  =  "output"; 
eqRatio. IOstatus  =  "input"; 

} 

// _ end  of  additions  -IA- 

}  ' 

else  if(  name  ==  "switchHotLoss"  )  { 
if  (  switchHotLoss  ==  "INPUT"  )  { 
dPqPRayleigh.IOstatus  =  "input"; 

} 

else  if  (  switchHotLoss  ==  "input"  ){  switchHotLoss  =  "INPUT";  } 
else  { 

dPqPRayleigh.IOstatus  =  "output"; 

} 

} 

else  if(  name  ==  "switchAud"  )  { 
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adPqPAud.IOstatus  =  "inactive"; 
sdPqPAud.IOstatus  =  "inactive"; 
if(  switchAud  ==  "AUDIT"  )  { 
adPqPAud.IOstatus  =  "input"; 
sdPqPAud.IOstatus  =  "input"; 

} 

} 

} 


// - 

//  *******  PERFORM  ENGINEERING  CALCULATIONS  ******* 

// - 

void  calcPreLossQ  { 


// - 

//  Check  to  see  if  the  pressure  sockets  are  empty,  if  not  thenexecute 

// - 

if  (  !S_dPqP.isEmpty()  )  { 

S_dPqP.execute(); 

}  ” 

dPqP  =  dPqPBase  *  s  dPqP  +  a  dPqP;  //  calculate  pressure  losses  (dry  duct  and 
Valve) 

if(  switchDes  ==  "OFFDESIGN"  )  { 
if(  switchAud  ==  "AUDIT"  )  { 

dPqP  =  dPqP  *  s  dPqPAud  +  a  dPqPAud; 

} 

} 

//comment  -IA-  Collect  total  enthalpy  at  inlet 
real  hin  =  Fl  l.ht; 

real  Pin  =  (  1  -  dPqP  )  *  Fl_I.Pt;  //coment  -IA-  apply  pressure  losses  as  calculated 
above 

//comment  -IA-  copy  flow  to  combustor  flow 
Fl_Icomb.copyFlowStatic(  "FI  I"  ); 

Fl_Icomb.setTotal_hP(  hin,  Pin  ); 


} 

void  calcBurn()  { 
real  TtCombOutTemp; 
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real  htStoich; 
real  WFuelLimit; 
real  WFuelHeat; 

Fl_Ocomb.copyFlow(  "Fllcomb"  ); 

// - 

//  Efficiency 

// - 

if  (  !S_eff.isEmpty()  )  { 

S_eff.execute(); 

}  ' 

eff  =  effBase  *  s_eff  +  a_eff; 

effChem  =  effChemBase  *  seffChem  +  aeffChem; 

// - 

//  Burn 

// - 

Fl_Ocomb.bum(  "Fu  I",  eff); 

// - 

//  if  inputting  a  PW  type  of  efficiency  adjust  the  temperature 

// - 

if  (  effChem  <  1.0)  { 

TtCombOutTemp  =  effChem  *(  Fl  Ocomb.Tt  -  Fl  Icomb.Tt  )  + 
FlIcomb.Tt; 

Fl_Ocomb.setTotalTP(  TtCombOutTemp,  Fl_Icomb.Pt );  //  use  Pin 

} 


void  calcRayleighLoss()  { 

flagRayleighChoked  =  0; 
flagRayleighLossTooMuch  =  0; 

PqPRayleigh  =  1.0; 
PqPRayleighError  =  0.0; 


// - 

//  self-convergent  iteration  loop  for  internal  momentum  pressure  drop  calc 

// - 
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for(  countRayleigh=0;  countRayleigh<=countRayleighMax; 
countRayleigh++)  { 


// - 

//  input  or  output  dPqPRayleigh 

// - 

if(  switchHotLoss  ==  "INPUT"  )  { 
PqPRayleigh  =1.0-  dPqPRayleigh; 

: 

else  if(  switchHotLoss  ==  "CALCULATE"  )  { 
dPqPRayleigh  =1.0-  PqPRayleigh; 

} 


// - 

//  calculate  momentum  pressure  drop 

// - 

real  PtCombOut  =  PqPRayleigh  *  Fl_Icomb.Pt; 

Fl_Ocomb.setTotal_hP(  FlOcomb.ht,  PtCombOut ); 


// - 

//  Check  momentum  pressure  drop 

// - 

PqPRayleighNew  =  PqPRayleigh; 

if  (  switchHotFoss  ==  "CAFCUFATE"  )  { 


// - 

//  make  this  thing  a  constant  area  burner 

// - 

FlOcomb.A  =  FlIcomb.A; 
flagRayleighChoked  =  0; 
if(  Fl_Ocomb.MN  >  1.0  )  { 

//  when  MN  >1.0  FlowStation  static  calc  is 
//  not  consistent  with  Area 

//  Fl_Ocomb.MN=  1.0; 

//  do  not  do  this  -  creates  major  iteration  problems 
flagRayleighChoked  =  1 ; 
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// - 

//  Calculate  the  exit  static  pressure  from  the  momentum  equation 
//  assume  the  fuel  has  the  same  velocity  as  the  entrance  flow 

// - 

real  PsMomMethl; 

PsMomMethl  =  Fl_Icomb.W*Fl_Icomb.V  -  Fl_Ocomb.W*Fl_Ocomb.V; 
PsMomMethl  =  PsMomMethl /CGRAVITY; 

PsMomMethl  =  PsMomMethl  +  Fllcomb.Ps  *  FI  Icomb.A; 

PsMomMethl  =  PsMomMethl/Fl_Ocomb.A; 
real  PsMomMeth2; 

//PsMomMeth2  =  Fl_Ocomb.W*Fl  Icomb.V; 

PsMomMeth2  =  Fl_Icomb.W*Fl_Icomb.V  +  Wfuel*Fl_Icomb.V*fuelFractV; 
PsMomMeth2  =  PsMomMeth2/C_GRAVITY; 

PsMomMeth2  =  PsMomMeth2  +  Fl_Icomb.Ps  *  FlIcomb.A; 

PsMomMeth2  =  P s M o m M c t h 2/ F I  O c om b A ; 

PsMomMeth2  = 

PsMomMeth2/(1.0+Fl_Ocomb.gams*Fl_Ocomb.MN*Fl_Ocomb.MN); 
//PsMomMethl  =  PsMonMeth2; 

// - 

//  Note  Methl  =  Meth2  when  MN  <=  1 .0 

//  Use  Meth2  -  seems  more  stable  the  Methl  when  MN  >1.0 

// - 

PqPRayleighNew  =  (PsMomMeth2/Fl_Ocomb.Ps)  *  PqPRayleigh; 


//  Check  against  tolerance 

PqPRayleighError  =  PqPRayleighNew  -  PqPRayleigh; 
if(  abs(PqPRayleighError)  <  tolRayleigh  )  {  break;  } 

//  Bounding  of  PqPRayleigh  movement  to  PqPRayleighStep 
real  sign; 

sign  =  PqPRayleighError/abs(PqPRayleighError); 

PqPRayleighDelta  =  sign  * 

min(abs(PqPRayleighError),  PqPRayleighStep); 

PqPRayleighNew  =  PqPRayleigh  +  PqPRayleighDelta; 

//  Lower  limit  of  PqPRayleigh  -  limit  too  much  loss  to  PqPRayleighMin 
if(  PqPRayleighNew  <  PqPRayleighMin  )  { 
if(  HagRayleighLossTooMuch  ==  1  )  { 

ESOreport(  1023 901, "Rayleigh  pressure  loss  limited,  too  much  loss",  FALSE  ); 
break; 
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} 

PqPRayleighNew  =  PqPRayleighMin; 
AagRayleighLossTooMuch  =  1; 

} 

else  { 

AagRayleighLossTooMuch  =  0; 

} 

/* 

//  debug  info 

cout «  Fl  Ocomb.A  «  " 
cout «  Fl  Ocomb.MN  «  " 
cout «  FlOcomb.Ps  «  " 
cout «  PsMomMethl  «  " 
cout «  PsMomMeth2  «  " 
cout «  PqPRayleigh  «  " 
cout «  PqPRayleighNew  «  " 
cout «  endl; 

*/ 


// - 

//  check  for  convergence 

// - 

if(  countRayleigh  >=  countRayleighMax  )  { 

ESOreport(  1023 901, "Rayleigh  iteration  failed  to  converge,  counter  exceed  max", 
FALSE  ); 

break; 

} 

PqPRayleigh  =  PqPRayleighNew; 


} 

if(  AagRayleighChoked  ==  1  )  { 

ESOreport(  1 023901, "Rayleigh  Fl  Ocomb.MN  exceed  choked 
condition",  FALSE  ); 

} 

} 


void  calculateQ  { 
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// - 

//  Prebuming  pressure  loss 

// - 

calcPreLoss();  //  creates  Fllcomb,  applies  pre-losses 

real  FARin  =  FlIcomb.FAR; 
real  WARin  =  FlIcomb.WAR; 

// - Added  6  Feb  2007  -  IA - 

if  (F1I.MN  ==  0.  &&  Fl_I.Aphy  ==  0.){ 

FIJcomb.MN  =  0.4; 

Fl_Icomb.setTotal_hP(Fl_Icomb.ht,  Fl_Icomb.Pt); 

} 

// - End  Additions  6  Feb  2007 

// - 

//  Pre-calculate  Burning  to  obtain  enthalpy,  burned  fuel  attrib. 

// - 

if  (  switchBum  ==  "FAR"  )  { 

// - 

//  detennine  the  fuel  weight  flow  from  the  input  FAR 

// - 

Wfuel  =  (  Fl  lcomb.W  /(  1 .  +  FARin  +  WARin))*(  FAR  -  FARin  ); 
Ful.Wfuel  =  Wfuel; 

eqRatio  =  FAR/Fu_I.FARst;  //  Added  5  Feb  2007  -  IA 

calcBum(); 

calcRayleighLoss(); 

TtCombOut  =  FlOcomb.Tt; 


} 

//mmmmmmmmmmmmmmmmmmmmmmmmmmmmmm# 

II  Added  5  February  2007  -  IA 
//  do  an  equivalence  ratio  calculation 
else  if  (switchBurn  ==  "EQRATIO")  { 

FAR  =  eqRatioH5Fu_I.FARst; 

Wfuel  =  (  Fl  lcomb.W  /(  1 .  +  FARin  +  WARin))*(  FAR  -  FARin  ); 
Fu_I. Wfuel  =  Wfuel; 

calcBum(); 

calcRayleighLoss(); 
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TtCombOut  =  FI  Ocomb.Tt; 

} 

// - 

//make  a  flow  station  that  has  properties  of  cold  vitiated  air... 

// - 

Fl_Vit.copyFlowStatic("Fl_Ocomb"); 

F1_V it.  setT otalTP(Fl  Icomb .  Tt,  Fl_Icomb.Pt); 

// - 

//  copy  inlet  flow  for  pure  air  reference  to  be  used  later 

// - 

//Take  a  snapshot  of  air  after  it  has  entered  the  detonation  tubes 
Fl_IcombAir.copyFlowStatic("Fl_Icomb"); 

//  Copy  input  flow  properties  for  internal  bypass  flow 
//  -  W  set  later 

Fl_02.copyFlow("Fl  Icomb  Air"); 

// - 

//  On-design  loop 

// - 

if  (switchDes  ==  "DESIGN")) 

// - 

//  Initialize  iterated  variables 

// - 

real  uCJ,  a_l,  rhoVit,  freq,  PcqPi,  errors; 

real  gamt,  Cpt,  beta,  MCJ2,  PcqPi2;  //  average  (static  gamma,  Cp) 
real  Atube,  Vtube;//,  mCycle,  Wtube; 
real  MFP,  Wvalve,  gma_I; 

real  mFillAir,  mPurgeAir,  mPureAir;  //tauVO,  WvalveOpen, 
real  tDetonation,  tDetProp,  tBlowdown,  tPurge,  tF ill,  iVel; 
real  gam_s,  gmm_fc;//,  a_inlet; 
real  WtotAir,  Wbypass; 

int  count; 

// —  initiated  but  not  iterated - 

//static  density  of  cool  vitiated  fluid 
rhoVit  =  FI  Vit.rhot;  //(lbm/ftA3) 
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//  speed  of  sound  in  pure  air,  stagnated  in  detonation  tube 

//  that  the  detonation  wave  propogates  in  to 

a_l  =  sqrt(Fl_Icomb.gamt*Fl_Icomb.Rt*Fl_Icomb.Tt*25037.); 


//================================================= 

//  Calculate  Chapman- Jouguet  Mach  number  for  wave  as  described 
//  in  Heiser  and  Pratt 

//================================================== 

//***  input  variables:  // 

//***  output  variables:  //MCJ,  deltaS,  qadd// 

//***  Flow  Stations:  //Fl_Ocomb,  Fl_Icomb  // 

//  local  variables:  //gamt,  Cpt,  qadd,  beta,  MCJ2  // 

// - Arithmetically  average  specific  heats - 

gamt  =  (Fl_Ocomb.gamt  +  Fl_Icomb.gamt)/2.0;  //  arithmetic  mean  of 
//  gamma  for  stopped 
//  fluid 

Cpt  =  (FlOcomb.Cpt  +  Fl_Icomb.Cpt)/2.0;  //  arithmetic  mean  of  Cp 

//  for  a  stopped  fluid 

// - Calculate  heat  addition  per  Heiser-Pratt  cycle - 

//  calculate  non-dimensional  heat  addition 

qadd  =  (Fl  Ocomb.ht  -  Fl_Icomb.ht)/(Cpt*Fl_Icomb.Tt); 

// - Calculate  Chapman- Jouget  Mach  number - 

beta  =  (gamt  +  1.0)*qadd+1.0; 

MCJ2  =  beta  +  sqrt(  beta*  *2  -  1.0  ); 

MCJ  =  sqrt(MCJ2); 

// - Calculate  Entropy  gain  based  on  CJ  detonation - 

deltas  =  Cpt*(-log(MCJ2*((gamt+l  .0)/ 

(1 .0+gamt*MCJ2))**((gamt+l  .0)/gamt))  ); 

// —  calculate  the  pressure  rise  using  the  H  &P  method - 

PcqPi  =  (1.0+  gamt*MCJ2)/(gamt+1.0); 
uCJ  =  a_l*MCJ; 


// - Calculate  tube  volume  and  Area 

Atube  =  (PI/4.)*dTube* *2/144.;  //  ftA2 
Vtube  =  Atube*(lTube/12);  //  ftA3 
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// - calculate  the  valve  inlet  mass  flow  rate - 

gma_I=Fl_Icomb  Air.  gamt; 

MFP  =  Mvalve*sqrt(  (gma_I*32.174)/(Fl_IcombAir.Rt*778.16)  ) 

*  ( 1 . +( gma_I- 1  ,)/2 .  *  M  valve  *  *  2)  *  *  ( (  gma_I+ 1 .  )/(2 .  *  ( 1 .  -gma_I))) ; 
Wvalve  =  (Fl_IcombAir.Pt/sqrt(Fl_IcombAir.Tt)) 

*(Atube*  144.*ARvalve)*MFP; 


// - 

//  On-Design:  Calculate  bypass  ratio 

// - 

//***  input  Variables:  //dTube,  lTube,  n_tubes,  fillFrac 
//  //  purgeFrac, 

//***  iterated  Variables  //  freq 
//***  output  Variables:  //  iBPR 

//***  local  variables:  //WfillAir,  WpurgeAir,  WpureAir,  WtotAir 
//  Wbypass,  WpurgeAir,  Wvit,  // 

//***  Flow  Stations:  //  Fl_IcombAir,  Fl_Icomb,  Fl_Iprg,  FI  Vit,  // 

// - Calculate  the  split  and  partition  of  flow - 

//  amount  of  air  that  will  be  mixed  with  fuel  -  one  tube 
mFillAir  =  Vtube*(rhoVit*fillFrac)/(l.+FAR); 

//  amount  of  air  that  will  purge  during  each  cycle  -  one  tube 
mPurgeAir  =  Vtube*(Fl  IcombAir.rhos*purgeFrac); 

//  total  air  per  cycle  flowing  though  one  tube 
mPureAir  =  mFillAir  +  mPurgeAir; 

// - 

//  Timing  -  calculate  frequency  -  Do  I  need  to  put  this  at  the  end? 

// - 

//***  input  Variables:  //  DDT,  tValve,  Ltube,  ff,  pf,  tCycle 

//***  iterated  Variables:  //  uCJ,  PcqPi 

//***  output  Variables  //  tCycle,  tauValveOpen,  freq 

//***  local  variables:  //  tDetonation,  tDetProp,  tBlowdown,  tPurge, 


//  //  tFill 

// - 

// - Detonation  time 


//  DetProp  time  is  relatively  independant  of  fill  fraction... 
//  tDetProp=  lTube/(uCJ*12); 
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//DDT  is  input,  tDetonationPropogatio  calcd  (may  need  to  iterate) 

//  tDetonation  =  DDT  +  tDetProp; 

// - Blowdown  time - 

//  assume  choked  flow  at  tube  exit  and  calculate  blowdown  based  on 
//  draw-down  time  of  a  pressurized  tank  calculated  on  pressure 
//  differential 

//  gam_s  =  Fl_IcombAir.gams;  //  larger  gamma  is  more  conservative 
//  gmm  fc  =  ((gam_s  +  l.)/2.)**(-(gam_s+l.)/(  2.*(gam_s-l.))  );  // 

//####  tBlowdown:  Use  —1/2  calcd  pressure  (to  match  experimental  data) 

//  we'll  use  CJ  det  wave  velocity  as  the  speed  of  sound  in  the  gas 
//  since  a  cannot  be  directly  calc'd 
//  note  tBlowdown  is  proportional  to  tube  length 
//  tauBlDn  is  proportional  to  tube  length... 

//  tBlowdown  =  (log(0.4*PcqPi)/gmm_fc)*(lTube/uCJ); 

// - Fill  an(i  Purge  time - 

//  Use  the  choked  flow  at  valve  inlet  and  the  mass  flow  rate  as 
//  calculated  outside  the  loop  to  calculate  fill  time  (m /  mdot) 

//  tPurge  =  tValve  +  mPurgeAir/Wvalve;  //(s) 

//  tFill  =  tValve  +  mFillAir/Wvalve;  //(s) 

//Improvement  could  be  made  by  calculating  vitiated  air  velocity... 

// - Cycle  Time  output  calculation - 

//  tCycle  =  tDetonation  +  tBlowdown  +  tPurge  +  tFill; 

//  tauValveOpen  =  (tPurge+tFill)/tCycle; 

fireq  =  1  ./tCycle;  //tCycle  is  user  input 

//cout  «  "\n  \n  tDetonation,  tBlowdown,  tPurge,  tFill  PcqPi"«"  "«  tDetonation  «" 
"«  tBlowdown«"  "«  tPurge«"  "«  tFill«"  fireq"  «  1/tCycle  «  "  "  «  PcqPi  « 
endl; 


// - Set  total  mass  flow  through  tubes - 

WtotAir  =  mPureAir*n_tubes*freq; 

//  steady-state  flow  rate  into  tubes 

//  conservation  of  mass  check 
if  (WtotAir  >  F1_I.W)  { 

fillFrac  =  fillFrac*(Fl  I.W/WtotAir); 
purgeFrac  =  purgeFracHs(Fl_I.W/WtotAir); 

mFillAir  =  Vtube*(rhoVit*fillFrac)/(l.+FAR); 
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//  amount  of  air  that  will  be  mixed  with  fuel  -  1  tube 


mPurgeAir  =  Vtube*(Fl_IcombAir.rhos*purgeFrac); 

//  amount  of  air  that  will  purge  during  each  cycle  - 1  tube 

mPureAir  =  mFillAir  +  mPurgeAir; 

//  total  air  per  cycle  flowing  though  one  tube 

WtotAir  =  F1I.W; 

//  cerr  «  "ATTENTION  !pf  &  ff  changed  to:  "  «  purgeFrac  «  "  "  « 

fillFrac  «  endl; 

ESOreport(  2222 100, "Purge  and  fill  fractions  changed  in  order  to  maintain 
conservation  of  mass  through  the  engine",  FALSE  ); 

//break; 


} 

// - Set  iBPR - 

Wbypass  =  (F1_I.W  -  WtotAir) *flowby; 

//  steady-state  flow  rate  sent  to  bypass 

iBPR  =  Wbypass/WtotAir; 

//  steady-state  internal  PDC  bypass  ratio 

iBPRdes  =  iBPR; 

// - Set  bypass  exit  flow  SPLIT - 

F1_02.W  =  Wbypass; 

F103.W  =  (F1I.W  -  W totAir)-((Fl_I. W -W totAir) * flowby); 

// - Set  purge  and  fill  stations  PARTITION - 

FlIprg.copyFlowStaticC'FlIcombAir"); 

//  copy  flow  for  purge  function 

// -  PURGE  AIR - 

Fl_Iprg.AphyDes  =  (Atube*144)*n_tubes;  //Set  physical  area 
Fllprg.W  =  m  Purge  A  i  r*  frcq  *  ntubes ;  //  set  mass  flow  rate 

// - FILL  AIR - 

Fl_Icomb.copyFlow("Fl_IcombAir"); 

FlIcomb.AphyDes  =  Atube*  1 44.*n_tubes*tauValveOpcn; 

//  Actual  area  is  multiplied  by  tauVO  to  get  equivalent 
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//  area.  -  Fluid  flows  steadily  through  this  area 

Fl_Icomb.W  =  mFillAir* ntubcs *  frcq ;  // 
Fl_Icomb.setTotal_hP(Fl_IcombAir.ht,  Fl_IcombAir.Pt); 

//sets  time-averaged  static  conditions 

// - 

//  Burning 

// - 

//  FAR  was  calculated  prior  to  enteringh  this  convergence  loop  - 
//  so  we  just  need  to  modify  Wfuel  based  on  changed  Fl_Icomb.W 
Wfuel  =  (  Fl  lcomb.W  /(  1 .  +  FARin  +  WARin))*(  FAR  -  FARin  ); 
Ful.Wfuel  =  Wfuel; 

calcBurn(); 

calcRayleighLossO; 

TtCombOut  =  FlOcomb.Tt; 

//================================================ 

//  Apply  Dyer-Kaemming  correction  to  obtain  tube  flow  at  exit 
//  (ignores  the  kinetic  energy  of  the  shock  wave.) 
//================================================ 

Fl  Ocomb .  setT otal_hS(Fl_Ocomb  .ht,  Fllcomb .  S+deltaS); 


} 

//  OFF-DESIGN  CODE  //added  Dec  09  -  CT 
if  (switchDes  ==  "OFFDESIGN"){ 

// - 

//  Initialize  iterated  variables 

// - 

//real  uCJ,  a_l,  rhoVit,  freq,  PcqPi,  errors; 

//real  gamt,  Cpt,  beta,  MCJ2,  PcqPi2;  //  average  (static  gamma,  Cp) 
//real  Atube,  Vtube;//,  mCycle,  Wtube; 

//real  MFP,  Wvalve,  gma_I; 

//real  mFillAir,  mPurgeAir,  mPureAir;  //tauVO,  WvalveOpen, 

//real  tDetonation,  tDetProp,  tBlowdown,  tPurge,  tFill,  iVel; 

//real  gam_s,  gmm_fc ;//,  a_inlet; 

//real  WtotAir,  Wbypass; 

//int  count; 
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// —  initiated  but  not  iterated - 

//static  density  of  cool  vitiated  fluid 
rhoVit  =  FIVit.rhot;  //(lbm/ftA3) 

//  speed  of  sound  in  pure  air,  stagnated  in  detonation  tube 

//  that  the  detonation  wave  propogates  in  to 

a_l  =  sqrt(Fl_Icomb.gamt*Fl_Icomb.Rt*Fl  Icomb.Tt*25037.); 


//================================================= 

//  Calculate  Chapman- Jouguet  Mach  number  for  wave  as  described 
//  in  Heiser  and  Pratt 

//================================================== 

//***  input  variables:  // 

//***  output  variables:  //MCJ,  deltaS,  qadd  // 

//***  Flow  Stations:  //Fl_Ocomb,  Fl_Icomb  // 

//  local  variables:  //gamt,  Cpt,  qadd,  beta,  MCJ2  // 

// - Arithmetically  average  specific  heats - 

gamt  =  (Fl_Ocomb.gamt  +  Fl_Icomb.gamt)/2.0;  //  arithmetic  mean  of 
//  gamma  for  stopped 
//  fluid 

Cpt  =  (FlOcomb.Cpt  +  Fl_Icomb.Cpt)/2.0;  //  arithmetic  mean  of  Cp 

//  for  a  stopped  fluid 

// - Calculate  heat  addition  per  Heiser-Pratt  cycle - 

//  calculate  non-dimensional  heat  addition 

qadd  =  (Fl  Ocomb.ht  -  Fl_Icomb.ht)/(Cpt*Fl_Icomb.Tt); 

// - Calculate  Chapman- Jouget  Mach  number - 

beta  =  (gamt  +  1.0)*qadd+1.0; 

MCJ2  =  beta  +  sqrt(  beta*  *2  -  1.0  ); 

MCJ  =  sqrt(MCJ2); 

// - Calculate  Entropy  gain  based  on  CJ  detonation - 

deltaS  =  Cpt*(-log(MCJ2*((gamt+1.0)/ 

(1 .0+gamt*MCJ2))**((gamt+l  .0)/gamt))  ); 

// —  calculate  the  pressure  rise  using  the  H  &P  method - 

PcqPi  =  (1.0+  gamt*MCJ2)/(gamt+1.0); 
uCJ  =  a_l*MCJ; 
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// - Calculate  tube  volume  and  Area - 

Atube  =  (PI/4.)*dTube**2/144.;  //  ftA2 
Vtube  =  Atube*(lTube/12);  //  ftA3 

// - calculate  the  valve  inlet  mass  flow  rate - 

gma_I=Fl_IcombAir.gamt; 

MFP  =  Mvalve*sqrt(  (gma_I*32.174)/(Fl_IcombAir.Rt*778.16)  ) 

*(1  A(gma_I-l.)/2.*Mvalve**2)**(  (gma_I+l.)/(2.*(l.-gma_I))); 
Wvalve  =  (Fl_IcombAir.Pt/sqrt(Fl_IcombAir.Tt)) 

*(Atube*  144.*ARvalve)*MFP; 


// - 

//  OFF-Design:  Calculate  bypass  ratio 

// - 

//***  input  Variables:  //dTube,  lTube,  n_tubes,  fdlFrac 
//  //  purgeFrac, 

//***  iterated  Variables  //  freq 
//***  output  Variables:  //  iBPR 

//***  local  variables:  //WfdlAir,  WpurgeAir,  WpureAir,  WtotAir 
//  Wbypass,  WpurgeAir,  Wvit,  // 

//***  Flow  Stations:  //  Fl_IcombAir,  Fl_Icomb,  Fl_Iprg,  Fl_Vit,  // 

// - Calculate  the  split  and  partition  of  flow - 

//  amount  of  air  that  will  be  mixed  with  fuel  -  one  tube 
mFillAir  =  Vtube*(rhoVit*fillFrac)/(l.+FAR); 

//  amount  of  air  that  will  purge  during  each  cycle  -  one  tube 
mPurgeAir  =  Vtube*(Fl_IcombAir.rhos*purgeFrac); 

//  total  air  per  cycle  flowing  though  one  tube 
mPureAir  =  mFillAir  +  mPurgeAir; 

// - 

//  Timing  -  calculate  frequency  -  Not  used  for  Thom  thesis,  frequency  is  input 

// - 

//***  input  Variables:  //  DDT,  tValve,  Ltube,  ff,  pf,  tCycle 

//***  iterated  Variables:  //  uCJ,  PcqPi 

//***  output  Variables  //  tCycle,  tauValveOpen,  freq 

//***  local  variables:  //  tDetonation,  tDetProp,  tBlowdown,  tPurge, 

//  //  tFill 

// - 
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// - Detonation  time - 

//  DetProp  time  is  relatively  independant  of  fill  fraction... 

//  tDetProp=  lTube/(uCJ*12); 

//DDT  is  input,  tDetonationPropogatio  calcd  (may  need  to  iterate) 

//  tDetonation  =  DDT  +  tDetProp; 

// - Blowdown  time - 

//  assume  choked  flow  at  tube  exit  and  calculate  blowdown  based  on 
//  draw-down  time  of  a  pressurized  tank  calculated  on  pressure 
//  differential 

//  gam_s  =  Fl_IcombAir.gams;  //  larger  gamma  is  more  conservative 
//  gmm_fc  =  ((gam_s  +  l.)/2.)**(-(gam_s+l.)/(  2.*(gam_s-l.))  );  // 

//####  tBlowdown:  Use  —1/2  calcd  pressure  (to  match  experimental  data) 

//  we'll  use  CJ  det  wave  velocity  as  the  speed  of  sound  in  the  gas 
//  since  a  cannot  be  directly  calc’d 
//  note  tBlowdown  is  proportional  to  tube  length 
//  tauBlDn  is  proportional  to  tube  length... 

//  tBlowdown  =  (log(0.4*PcqPi)/gmm_fc)*(lTube/uCJ); 

// - Fill  and  Purge  time - 

//  Use  the  choked  flow  at  valve  inlet  and  the  mass  flow  rate  as 
//  calculated  outside  the  loop  to  calculate  fill  time  (m/  mdot) 

//  tPurge  =  tValve  +  mPurgeAir/Wvalve;  //(s) 

//  tFill  =  tValve  +  mFillAir/Wvalve;  //(s) 

//Improvement  could  be  made  by  calculating  vitiated  air  velocity... 

// - Cycle  Time  output  calculation - 

//  tCycle  =  tDetonation  +  tBlowdown  +  tPurge  +  tFill; 

//  tauValveOpen  =  (tPurge+tFill)/tCycle; 

fireq  =  l./tCycle;  //frequency  is  input  (Thom  thesis) 


//cout  «  "\n  \n  tDetonation,  tBlowdown,  tPurge,  tFill  PcqPi"«"  "«  tDetonation  «" 
"«  tBlowdown«"  "«  tPurge«"  "«  tFill«"  fireq"  «  1/tCycle  «  "  "  «  PcqPi  « 
endl; 


// - Set  total  mass  flow  through  tubes 


125 


WtotAir  =  mPureAir*n_tubes*freq; 


//  steady-state  flow  rate  into  tubes 

//  conservation  of  mass  check 

if  (WtotAir  >  F1I.W)  { 

fillFrac  =  fillFrac*  (F1_I.  W/W  tot  A  i  r); 
purgeFrac  =  purgeFrac*  (Fl_I.W/Wtot  Air); 

mFillAir  =  Vtube*(rhoVit*fillFrac)/(l.+FAR); 

//  amount  of  air  that  will  be  mixed  with  fuel  -  1  tube 

mPurgeAir  =  Vtube*(Fl_IcombAir.rhos*purgeFrac); 

//  amount  of  air  that  will  purge  during  each  cycle  - 1  tube 

mPureAir  =  mFillAir  +  mPurgeAir; 

//  total  air  per  cycle  flowing  though  one  tube 

WtotAir  =  F1_I.W; 

//  cerr  «  "ATTENTION  !pf  &  ff  changed  to:  "  «  purgeFrac  «  "  "  « 

fillFrac  «  endl; 

ESOreport(  2222 100, "Purge  and  fill  fractions  changed  in  order  to  maintain 
conservation  of  mass  through  the  engine",  FALSE  ); 

//break; 


} 

// - Set  iBPR - 

Wbypass  =  (F1_I.W  -  WtotAir) *flowby; 
//  steady-state  flow  rate  sent  to  bypass 

iBPR  =  Wbypass/WtotAir; 

//  steady-state  internal  PDC  bypass  ratio 

//iBPRdes  =  iBPR; 

// - Set  bypass  exit  flow  SPLIT 

F1_02.W  =  Wbypass; 
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F1_03.W  =  (F1_I.W  -  WtotAir)-((Fl_I.W-WtotAir)*flowby);  //bleed  flow  needed 

//for  static  pressure 
//of  iBPR  entering 
//  Mixer39  to  equal 
//  static  pressure 
//entering  from  tubes 

// - Set  purge  and  fill  stations  PARTITION - 

FlIprg.copyFlowStaticC'FlIcombAir"); 

//  copy  flow  for  purge  function 

// -  PURGE  AIR - 

Fl  Iprg.AphyDes  =  (Atube*144)*n_tubes;  //Set  physical  area 
Fl  lprg.W  =  mPurgeAir*freq*n_tubes;  //  set  mass  flow  rate 

// - fill  air - 

Fl_Icomb.copyFlow("Fl  IcombAir"); 

FlIcomb.AphyDes  =  Atube*144.*n_tubes*tauValveOpen; 

//  Actual  area  is  multiplied  by  tauVO  to  get  equivalent 
//  area.  -  Fluid  flows  steadily  through  this  area 

Fl_Icomb.W  =  mFillAir* n_tubes *  frcq ;  // 

Fl_Icomb.setTotal_hP(Fl_IcombAir.ht,  Fl_IcombAir.Pt); 

//sets  time-averaged  static  conditions 

// - 

//  Burning 

// - 

//  FAR  was  calculated  prior  to  enteringh  this  convergence  loop  - 
//  so  we  just  need  to  modify  Wfuel  based  on  changed  Fl_Icomb.W 
Wfuel  =  (  Fl  lcomb.W  /(  1 .  +  FARin  +  WARin))*(  FAR  -  FARin  ); 

Fu_I. Wfuel  =  Wfuel; 

calcBum(); 

calcRayleighLoss(); 

TtCombOut  =  FlOcomb.Tt; 

//========================================================== 

//  Apply  Dyer-Kaemming  correction  to  obtain  tube  flow  at  exit 

//  (ignores  the  kinetic  energy  of  the  shock  wave.) 

//========================================================== 

Fl_Ocomb.setTotal_hS(Fl_Ocomb.ht,  Fl_Icomb.S+deltaS); 

} 
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//END  OFF-DESIGN  CODE 

// - 

//  Add  split  flows  back  to  combusted  flow 

// - 

Fl_Ocomb.add("Fl_Iprg");  //add  purge  flow  in  (uncorrected) 

//  Apply  corrections  to  the  flow  for  transition  to  steady  state... 

//  TTSS 

//—====—===——=====—=====— ========== 

//***  local  Variables:  //  Snew,  Pnew 

//***  Input  Variables:  //  deltaS,  TTSSeff,  TTSSdPqP 

//***  Flwo  stations:  //  FlOcomb,  Fl_Vit 

real  hnew,  Pnew;  // 

// - Calculate  new  Entropy  and  Pressure - 

//  eff  =  (dht)TTSF/(dht)comb  +  1 . 

//  current  h  -  ( h  gained)*(l.-eff) 

hnew  =  FlOcomb.ht  -  (FlOcomb.ht  -  FI  Icomb.ht)*(1.0-TTSSeff); 

Pnew  =  Fl_Ocomb.Pt*(  1 .0-TTSSdPqP); 

//End  of  12Jan2007  additinos  -  IA 

//mmmmmmmmmmmmmmmmmmmmmmmmmmmmm# 

Fl_01.copyFlow(  "Fl  Ocomb"  ); 

// - update  fluid  properties  based  on  new  Entropy  and  Pressure 

Fl_01.setTotal_hP(hnew,  Pnew);  //added  12Jan2007  -  IA 


// - 

//  Thennal  storage  calculations 

// - 

if  (  !S_Qhx.isEmpty()  )  { 
S_Qhx.execute(); 

} 

real  hout  =  F1_0 1  .ht  -  Qhx  /  F1_0 1 . W ; 
Fl_01.setTotal_hP(  hout,  Fl_01.Pt ); 


// - 

//  store  the  design  value  of  FAR  for  use  in  guessing 

// - 
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if  (  switchDes  ==  "DESIGN"  )  { 
FARDes  =  FAR; 

} 


// - 

//  register  the  appropriate  errors  at  build  time 

// - 

void  VCinit() 

{ 

ESOregCreate(  1023901,  8,  TRUE,  FALSE,  TRUE  );  //provisional 
ESOregCreate(  1093901,  8,  "",  TRUE,  FALSE,  TRUE  );  //provisional 
ESOregCreate(  2222100,  2,  "",  TRUE,  FALSE,  TRUE  );  //provisional 

} 


#endif 
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